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1 INTRODUCTION 

The dynamics of fluids is described by a classical field theory, whose origins lie in the nineteenth 
century, like Maxwell electrodynamics with which it shares some antecedents. The electromagnetic 
theory has enjoyed much development: it passed into quantum physics, and the resulting quantum 
electrodynamics served as a model for its non-Abelian generalization, Yang-Mills theory, which 
today is at the center of a quantum field theoretic description for fundamental physics. We believe 
that fluid dynamics can undergo a similar evolution and play a similar generative role in physics. 

Modern (quantum) field theory has expanded concepts and calculational possibilities beyond 
what was familiar to (classical) field theorists. One learned about higher and unexpected symme- 
tries, which also facilitate partial or complete integrability of the relevant differential equations. 
Topological and geometric concepts and structures, like solitons and instantons, were recognized 
as encoding crucial dynamical information about the models. New entities like Pontryagin den- 
sities and Chern-Simons terms entered into the description of kinematics and dynamics. Degrees 
of freedom were enlarged and unified by new organizing principles based on non-Abelian and su- 
persymmetries. Indeed application of field theory to particle physics has now evolved to a study 
of extended structures like strings and membranes, whose mathematical description bears some 
similarity to the theory of fluids. 

The novelties introduced in particle field theories can be also developed for the non-particle field 
theory of fluid dynamics. Correspondingly, fluid dynamics can illuminate some aspects of particle 
physics, especially as concerns the extended structures that these days are the focus of attention 
for particle physicists. 

We begin by addressing the question of why fluid mechanics would be interesting in its own right. 
Fluid mechanics, for most physical situations, is obtained from an underlying particle description by 
suitable averages of Boltzmann-type equations. Recall that the Boltzmann equation for the single 
particle distribution function f(X,P,t) is given by 

where X, P, refer to the coordinates and momenta, the phase space variables, of a single particle 
of mass m; F is the force acting on the particle. C(f) is the collision integral, which takes into 
account particle interactions. In the special case of the collisonless limit, i.e., with C = 0, the 
Boltzmann equation (1.0.1) is the equation for the distribution function for single particles obeying 
the standard classical equations of motion. Solving the Boltzmann equation is not very easy. 
The equilibrium distribution function is a solution of the equation; in particular, C(f) = for the 
equilibrium solution. The general strategy which has been used for solving equation (1.0.1) is to seek 
a perturbative solution of the form / = f(°'+f( 1 '+..., where f(°> = n p is the equilibrium distribution, 
appropriately chosen for bosons and fermions. Transport coefficients and fluid equations of motion 
can then be obtained from the perturbative corrections. 
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This approach has the virtue of simplicity and does capture many of the general features of the 
problem of deriving coarse-grained dynamics, fluid mechanics in particular, from the underlying 
particle dynamics. However, there are two limitations: First because we are only taking care of the 
single-particle distributions, and second because the treatment is essentially classical. 

Regarding the first point, one could indeed try to be more general by starting with the completely 
general iV-particle Liouville equation for the phase space distribution p(X n , P n ), n = 1, 2, N. The 
one-particle distribution function is then given by J dn^-i p(X„, P„), the two-particle distribution 
is given by j p(X n , P n ), etc., where d[iN denotes the phase space volume for N particles. The 

Liouville equation then leads to a hierarchy of kinetic equations, the so-called BBGKY (Bogolyubov- 
Born-Green-Kirkwood-Yvon) hierarchy, involving higher and higher correlated n-particle distribu- 
tion functions. (For the one-particle distribution function, we get the Boltzmann equation, but with 
the collision integral given in terms of the two-particle distribution function.) To be able to solve 
this infinite hierarchy of equations, one needs to truncate it, very often at just the single particle 
distribution function. Therefore, even though a more general formulation is possible, the feasibility 
of solving these equations limits the kinetic approach to dilute systems near equilibrium, where the 
truncation can be justified. 

Regarding the question of quantum corrections, the needed formalism is that of the Schwinger- 
Dyson equations with a time-contour, the so-called Schwinger-Bakshi-Mahanthappa-Keldysh ap- 
proach [1]. The Green's functions are defined by the generating functional 



where the time-integral goes from — oo to oo, folds back and goes from oo to — oo; Tq denotes 
ordering along this time-contour. M generically represents fields of interest, rf is a source function 
and po is the thermal density matrix. One can represent Z[rj] as a functional integral. 



The action is again defined on the time-contour. The Green's functions for which some of the 
fields are on the forward time line and some are on the reverse time line will represent the effect of 
statistical distributions. One has to solve the hierarchy of coupled Schwinger-Dyson equations which 
follow from (1.0.3), again truncating them at a certain level to get a description of nonequilibrium 
phenomena. In practice, one has to carry out semiclassical expansions to simplify these to the point 
where a solution can be found and again we have a formalism of limited validity. 

This discussion shows that we should expect that the regime of validity of the fluid dynamical 
equations derived within kinetic theory or within the time-contour approach is rather limited, 
basically a semiclassical regime for dilute systems not too far from equilibrium. However, fluid 
dynamical equations can also be derived from very general principles, showing that they have a 



Z[ V ] 



Trpo T c exp(il int + icf) ■ rj) 
Trpo 





(1.0.3) 
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much wider regime of validity, and, indeed in practice, we apply them over such a wider range. 
This is the 'universality' of fluid dynamics. It is this property which shows the value of a study of 
fluid mechanics in its own right, rather than its derivation from an underlying description in some 
approximation. It seems that a reconsideration of the whole setting and development of the ideas 
of fluid mechanics in the context of modern concepts in particle physics is entirely appropriate. The 
novelties introduced in particle field theories can be generalized to the non-particle field theory of 
fluid dynamics. Correspondingly, fluid dynamics can illuminate some aspects of particle physics, 
especially as concerns the extended structures that these days are the focus of attention for particle 
physicists. 

In this Introduction, we shall review Lagrange's and Euler's description of fluid kinematics and 
dynamics, and we shall describe the mapping between the two. We shall discuss the Hamiltonian 
(canonical) formulation, together with the associated Poisson brackets. We shall also present (con- 
figuration space) Lagrangians for fluid motion; in the Eulerian case a Clebsch parameterization 
is needed. Both nonrelativistic and relativistic systems in various spatial dimensions are treated. 
These are standard topics and they are well known [2]. However, in our presentation of this famil- 
iar material we shall approach the subject with an eye towards the various enhancements of fluid 
mechanics, which comprise our current research and which are reviewed in the remainder of this 
article. 



1.1 Lagrange and Euler descriptions of a fluid and the relationship 
between them. 

The Lagrange description of a fluid focuses on the coordinates of the individual fluid particles. 
These satisfy a Newtonian equation of motion (in the nonrelativistic case). On the other hand, 
in Euler's formulation, the fluid is described by a density p and velocity v, which are linked by a 
continuity equation, while Euler's equation describes the dynamics. Euler's method is akin to a 
classical field theory in physical space-time. 

In fact one may exemplify the two approaches, and the relation between them, already for 
a single particle, carrying mass m and located on the coordinate X(i), whose time evolution is 
governed by a force F. 

X(f) = -F(X(t)) (1.1.1) 

m 

(Over-dot denotes differentiation with respect to explicit time-dependence.) This is the "Lagrange" 
description (of course it is Newton's). Next, we introduce the Eulerian (single) particle density by 

p(t 7 r)=mS(X(t)-r). (1.1.2) 

The delta function follows the dimensionality of space so that the volume integral of p is m. Dif- 
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ferentiation (1.1.2) with respect to t leaves 

P(t,r) = ™^<W) - r)X*(f) 



■|-(x^(X(t)-r))=-A 



V%T) P (t,T) 



1.1.3) 



where the Eulerian velocity v is given by 

v(t, r) = X(t) with r = X(t). (1.1.4) 

Note that the velocity function v(t, r) is only defined at the point r = X(£); its value at other points 
being undetermined and irrevelent. Evidently the continuity equation is satisfied as a consequence 
of the above definitions. 

p(t,r) + V-j(t,r) = (1.1.5) 

j(t, r) = v(t, r)p(t, r) = X(t)mS(X(t) - r) (1.1.6) 

To arrive at the dynamical Euler equation, we differentiate j with respect to time. From (1.1.6) it 
follows that 

p(t, r)v(t, r) + v(t, r)p(t, r) = 

d 

X(t)mS(X(t) - r) + X(t)m— ,5(X(t) - r)X^t). (1.1.7a) 
Use of the continuity equation on the left side and Newton's equation (1.1.1) on the right leaves 
p(t, r)v(t, r) - v(t, r)V ■ (v(t, r)p(t, r)) = 

F(X(t))5(X(t) - r) - A(i(i)#(f)m<5(X(f) - r)Y (1.1.7b) 



Cancelling common terms gives Euler's equation (for a single-particle "fluid"!). 

vUt, r) + v j (t, r)d j v i (t, r) = —Fir) (1.1.7c) 

m 

[dj denotes a derivative with respect to the components of a spatial vector, which can be X or r, or 
below, x. If context does not determine unambiguously which vector is involved, it will be specified 
explicitly, as in (1.1.7a,b).] Strictly speaking, (1.17c) only holds at the point r = X. Suitable 
continuations of the function v(t, r) away from the point r = X can be found that make (1.1.7c) 
hold everywhere. Such continuations represent the velocity field of a fictitious accompanying fluid, 
of which X is one particle. For a complete dynamical description an appropriate expression for the 
force F still needs to be given. 

The same development holds for a collection of iV particles: (1.1.1) becomes replaced by 

X n (0 = —F(x 1 (t),...,X n (t)). (1.1.8a) 
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The particle label n ranges from 1 to N. We shall take the particles to be identical. Therefore the 
mass does not carry the n label, and the force F n has a functional dependence on X n independent 
of n, and is symmetric under exchange of the remaining N — 1 particle coordinates. 

X n (t) = -f(x„(0; {X k (t), k + n}) (1.1.8b) 



m 

The Eulerian mass density, velocity and current are defined as 

N 



p(t,T)=m^26(X n (t)-T), (1.1.9) 



n=l 

N 



j(t, r) = v(t, r)p(t, r) = m^ X n (t)«5(X n (t) - r). (1.1.10) 

n=l 

Similarly to the single-particle case, the function v(t,p) is defined only at the points r = X n (t). 
Evidently the continuity and Euler equations continue to hold. 

For the true fluid formulation, we promote the discrete particle label n to a continuous label 
x and the Lagrange coordinate X n (t) becomes X(t, x). Frequently x is specified by the statement 
that it describes the fluid coordinate X at initial time t — 0, i.e. 

X(0,x) = x. (1.1.11) 

Thus x is the comoving coordinate. Dynamics is again Newtonian. 

X(t, x) = — F(X(t, x)) (1.1.12) 

m 

The density and velocity are now defined by 

p(t,r) = p Jdx6(X(t,x)-r), (1.1.13) 

j(t,r) = v{t,r)p(t,r) = p JdxX{t,x.)6{X{t,x.)-r). (1.1.14) 

The integration is over the entire relevant volume, be it 1-dimensional, 2-dimensional or 3-dimensional. 
(The dimensionality of the measure will be specified only when formulas are dimension specific.) po 
is a background mass density, so that the volume integral of p is the total mass. In the above, we 
assume that the particles that are nearby in space have similar velocities and thus the function of 
discrete points v(t, X n ) goes over to a smooth continuous function v(£, r). This is the fundamental 
assumption in classical fluids. The more general case calls for a phase space density (Boltzmann) 
description of the collection of particles, which takes us away from the realm of classical Lagrange 
fluids. 
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The continuity and Euler equations follow as before, with the latter reading 

v% r) + v j {t, r)djV l {t, r) = — F*(r). (1.1.15) 

The determination of the force term in the above equation requires some discussion. For external 
forces, F(r) is simply a known function of space. For internal particle forces, however, the force could 
be a nonlocal function depending on the full distribution of particles in space. We shall postulate 
that interparticle forces are short-range and thus depend only on the distribution of nearby particles. 
For such forces, F would depend only on the density of particles and its derivatives at the point r. 
Under additional assumptions of isotropicity, F will be proportional to the gradient of the density, 
Vp. In this case, a standard argument shows that the right hand side of (1.1.15) should be set 
to — i VP, where P is the pressure [2]. Thus once an equation of state is given, i.e. once the 
dependence of P on p is known, we have a self-contained system of equations: the continuity and 
Euler equations. 

p(t,r) + V- (v(t,r)p(t,r)J =0 (1.1.16) 

v(t, r) + (v(t, r) ■ v(t, r) = --VP{p) (1.1.17) 

A fluid obeying these equations is called "perfect". Later this formalism will be generalized to 
account for an internal symmetry. For that purpose, it will be useful to delineate explicitly the 
effect of the 6 functions in (1.1.13) and ((1.1.14). In the course of the x integral, x becomes 
evaluated at a function x{t, r ) ; which is inverse to X(t, r). 



X(t, X (t,r)J = r (1.1.18a) 
x(a(t,x)) = x (1.1.18b) 



(Uniqueness is assumed.) Thus the x integration sets X equal to r and also there is a Jacobian, 
det^j- for the X — > r transformation. 

Consequently from (1.1.13) and (1.1.14) it follows that 

9 = p0 J^dW\ — ' (L1 ' 19) 

l uet 8x3 |x=x 

v = X| x=x . (1.1.20) 

In other words X effects a diffeomorphism x — > X = r, while x ac ts similarity for r — > x. The 
interchange between a dependent variable X and an independent variable r is called a hodographic 
transformation. 
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1.2 Lagrangian, Hamiltonian formulations and symmetries of dynamics 

(i) Lagrangian and Hamiltonian functions for Lagrange fluid mechanics 

Since the Lagrange method is essentially Newtonian, Lagrangian and Hamiltonian formulations for 
(1.1.12) are readily constructed provided the force is derived from a potential V(X). 



1 



dx\ -mX 2 - V(X) I (1.2.1a) 



V(X)) 



2 v 7 

F(X) = -VV(X) (1.2.1b) 

H L = J dx(J^P 2 + V(X)J (1.2.2a) 

P = mX (1.2.2b) 

A canonical formulation for the Lagrange description of fluid motion follows in the usual way, so 
that bracketing with Hi and using conventional Poisson brackets for X and P reproduces the equa- 
tion of motion (1.1.12). 



(ii) Diffeomorphism symmetry of Lagrange fluid mechanics 

In the discrete antecedent to the continuum formulation there is the obvious freedom of renaming 
the n label. The continuum version of this freedom manifests itself in that the Lagrange formulation 
of fluid dynamics enjoys invariance against volume-preserving diffeomorphisms of the continuous 
label x. 

An infinitesimal diffeomorphism of x, generated by a infinitesimal function f(x) reads 

<S/x = -f(x), (1.2.3a) 
and this is volume-preserving when f is transverse. 

V-f = (1.2.3b) 
Provided the Lagrange coordinate transforms as a scalar 

5 f X(t, x) = f (x) ■ VX(t, x), (1.2.4) 

Ll is invariant. 

5 f L L = JdxfmX'fdjX' - -!L-V(X).r<),X^ 

dxfdJ^X 2 - V(X) ) (1.2.5) 
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The last expression vanishes after partial integration by virtue of (1.2.3b). (Surface terms are 
ignored.) Noether's theorem then gives the constant of motion associated with the relabeling 
symmetry (volume-preserving diffeomorphism) . 

C f = J dxX'f'djX* (1.2.6) 

Since f is an arbitrary transverse function it can be stripped away from (1.2.6). Explicit formulas 
reflect the spatial dimensionality of the system. In three dimensions we can present the transverse 
f 1 as e^ k djf k , f k arbitrary, leading to conserved local vector quantities. 

q 3) = e^djX^kX 6 (1.2.7) 

In two dimensions f l involves a scalar function, j % = e^djf, and for planar systems the conserved 
quantity is a local scalar. 

C (2 ) = e lj diX k djX k (1.2.8) 

Finally, in one dimension, a "transverse" function is constant, so the constant of motion for lineal 
systems remains integrated. 

C (1) = f dx l Xd x iX = [ dXX (1.2.9) 



(iii) Hamiltonian function for Euler fluid mechanics 

A Hamiltonian for the continuity and Euler equations (1.1.16), (1.1.17) is obtained from Hl (1.2.2a) 
by transforming to Eulerian variables. First, however, to ensure proper dimensionality the kinetic 
term is divided by the spatial volume so that m is replaced by p . Also to reproduce the pressure 
form of the forces, we write V(X) as a function of the Jacobian, det^-. Thus Hl reads 

// 1 av i \ 

dx(-p X 2 + V(\det— \/po)Y (1.2.10) 

The transformation to Eulerian variables is effected by multiplying (1.2.10) by unity, in the form 
J dr5(X. — r), and interchanging orders of integration. In this way Hl becomes He- 

H e= [ ^pv 2 + ^V(-)^ (1.2.11) 



.2 Po P 

Agreement with the pressure expression for the force, as in (1.1.17), is achieved when we identify 

P(p) = -- V(l/p) (1.2.12) 
Po 

(Dash denotes derivative with respect to argument.) In the subsequent we drop the subscript E on 
H and rename ■^V(-) as V(p). Thus the Euler Hamiltonian reads 

H = J rfrQpv 2 + F(p)^), (1.2.13) 
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and the pressure is Legendre transform of V. 

P(p) = pV'(p)-V(p) (1.2.14) 



Conventional nomenclature for V'(p) is enthalpy and ^P'(p) = a/ pV"(p) is the sound speed s. 

Bracketing the Hamiltonian with p and v should generate equations (1.1.16) and (1.1.17). To 
verify this, we need to know the brackets of p, v with each other. These may be obtained from the 
canonical brackets in the Lagrange formulation, the only non-vanishing one being 

{X l (x), X J (x')} = — 5 ij 5(x - x'). (1.2.15) 
Po 

Using the definitions of p and j in terms of X and X, eqs. (1.1.13) and (1.1.14), as well as the 
canonical brackets (1.2.15), determines the brackets of p and j. 

{p(v),p(r')} = (1.2.16) 
{/(r),p(r')} = p(r)c^(r-r') (1.2.17) 
{f(r),f(r')} = j*(r)d i 6(T-T f )+j i (T f )d j 5(r-T f ) (1.2.18) 

Since j = vp these in turn imply that the brackets for p and v take the form [3] 

K(r),p(r')} = ^(r-r / ), (1.2.19) 

Mr), vi(v>)} = -^-6(r - r'), (1.2.20) 
p(r) 



where 



di Vj{r) - d^r) (1.2.21) 



is called the fluid vorticity. [In the above expressions, (1.2.15)-(1.2.21), all quantities are at equal 
times, so the time argument is omitted.] Of course the Jacobi identity is satisfied by the brackets. 

It is now straight forward to verify from (1.2.16), (1.2.19)-(1.2.21) that bracketing with H 
reproduces the equations of motion (1.1.16) and (1.1.17). 

p = {H,p} = -V -iyp) (1.2.22) 
v = {H, v} = -(v • V)v - W(p) (1.2.23) 

These equations may also be presented as continuity equations for a (nonrelativistic) energy 
momentum tensor. The energy density 

£ = Ipv 2 + V = T 00 , (1.2.24) 

together with the energy flux 

T^ = p^(iv 2 + y'), (1.2.25) 
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obey 

f°° + djT j0 = 0. (1.2.26) 
Similarly the momentum density, V, which in the nonrelativistic theory coincides with the current, 

V { = pv i = T°\ (1.2.27a) 

and the stress tensor T u 

T ij = 5 ij (pV -V)+ pv i v j = 5 ij P + /wV, (1.2.27b) 

satisfy 

T oi + djT^ — 0. (1.2.28) 

Note that T m ^ T l ° because the theory is not Lorentz invariant, but T*- 7 = T Jl because it is invari- 
ant against spatial rotations. {T^ u is not, properly speaking, a "tensor", but an energy-momentum 
"complex".) Because V = j, the current algebra (1.2.18) is also the momentum density algebra. 

(iv) Symmetries of Euler fluid mechanics 

The above continuity equations and other specific properties of the energy-momentum tensor 
allow constructing constants of motion, which reflect symmetries of theory. The Hamiltonian = 
energy, 

B = f*S (tune-t—,), (1.2.29) 
is constant as a consequence of time-translation invariance, while the constancy of the momentum, 

P J dr V = Jdrj (space-translation), (1.2.30) 

follows from space-translation invariance. The index symmetry of T % i is a consequence of rotational 
invariance and ensures that the angular momentum, 

M ij = j dr {r l V j - r j V l ) (spatial rotation), (1.2.31) 

is constant. The identity of the momentum density and the current density allows construction of 
the Galileo boost constant of motion, 

B = tP-Jdrrp (velocity boost), (1.2.32) 

whose time-independence signals invariance against velocity boosts. Finally, the total number 

N = J drp (number), (1.2.33) 
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is also conserved, as a consequence of the continuity equation (1.1.16). Upon bracketing with each 
other, the constants (1.2.29)-(1.2.33) form the extended Galileo Lie algebra, and they generate 
Galileo transformations on p and v. Brackets with N vanish; N provides the central extension for 
the Galileo algebra in the {B, P} bracket. 

Further constants of motion are present for specific forms of V. In particular if 

2T 00 = Sij T ij , (1.2.34a) 
which in d spatial dimensions requires that 

V(p) = Xp 1+1 d , (1.2.34b) 
two more constants exist. They are the dilation, 

D = 2tH - j dr r V (dilation), (1.2.35) 
and the special conformal generator. 

K = t 2 H — tD — ^ J dr r 2 p (conformal transformation) (1.2.36) 

The latter two, together with H, form the 50(2, 1) Lie algebra of the non-relativistic conformal 
group, which they generate by bracketing. Together with remaining Galileo group elements this is 
called the Schrodinger group [4]. 
Finally, in any dimension when 

V(p) = A/p, (1.2.37) 

we are dealing with a Chaplygin gas and, as we shall discuss below, this model supports remarkable 
higher symmetries related to relativistic extended objects [5]. 

What is the Eulerian image for the volume-preserving diffeomorphism symmetries of the La- 
grange formulation discussed in Section 1.2 (ii)? First we note that the Eulerian variables p and 
v do not respond to the volume-preserving transformations on the Lagrange variable X. This 
is seen from (1.1.13), (1.1.14), (1.2.4). Also the associated constants of motion, summarized by 
(1.2.6)-(1.2.9) do not possess in general simple Eulerian counterparts. But some do. 

Let us begin with the 3-dimensional quantity in (1.2.7) and integrate it over a surface S in the 
X parameter space bounded by a closed curve dS. This gives the constant 

A = j dS i e iik d j (d k X e X e ) = I dX <9; X e X e . (1.2.38a) 

Js JdS 

Upon performing the diffeomorphism x — > Xi (1 -2.38a) becomes expressed in terms of Euler variables 

as 

A= f dr-v= dS-uj, (1.2.38b) 

JdS JS 
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where dS is the (time-dependent) image of <9S under the diffeomorphism. The above quantity, 
known as the velocity circulation or vorticity flux, is therefore constant; this is Kelvin's theorem. 
[One can establish the result directly for (1.2.38b) from the equations of motion for Euler variables, 
provided one takes into account the time-dependence of the contour dS.} Similarly conserved 
velocity circulation exists also in 2-dimensional fluid mechanics; of course the integration surface or 
contour then lie in the plane. 

Additional important constants arise from the volume integrals of the C%, d = 1,2, 3, in (1.2.7)- 
(1.2.9). In three dimensions, we begin with (1.2.7) contracted with X m diX m and integrated over 
volume. 

(7 ( 3) = J d 3 x e ijk (X m d i X m )d j (X n d k X n ) 

d 3 x X m djX n e i]k diX m d k X n (1.2.39a) 

We use the identity e % i k diX m d k X n = e min ^j det and multiply (7(3) by unity in the form 
J d 3 r S(K — r) so that (7( 3 ) becomes 

/f)X n rM c)X p 

d 3 rv • (V x v) = / d 3 rv ■ u, (1.2.39b) 



(apart from an irrelevant factor of ±po)- Conservation of Ctg\, called the vortex helicity, also follows 
when the Euler equation (1.1.17) is applied to (1.2.39b). 

An expression like Crg\, with v an arbitrary 3- vector, is a well known mathematical entity, 
called the Chern-Simons term. In the last twenty years, Chern-Simons terms have come to play a 
significant role in physics and mathematics [6]. We shall have more to say about them, incarnated 
in various contexts, but always in odd- dimensional space. Indeed the one-dimensional constant 
(1.2.9), when written in terms of Euler variables reads 

C (1) = J dP-rv. (1.2.40) 

(Here "r" is not a radial coordinate, but lies on the real line.) This is just a 1-dimensional Chern- 
Simons term. (We shall further discuss Chern-Simons terms in Sidebar B, below.) 

Finally we turn to the planar constant (7(2) (1.2.8). Taking it to the M th power, multiplying by 
<5(X — r) and integrating over x and r gives 



Cg, = J d 2 x d 2 r(e ij d i X k d j X k ) M 5(X - r). (1.2.41a) 
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The 2x2 matrix identity £^^r = £ mk -§£^det converts the above to 

C& = [d 2 xd 2 r(^e mk ^-det^-) M 5(X-T) 
{2 > J \ dx l dX m dxi J v ; 



f> \ M r / \ m 

]2„„l ° „,k ^mk / „ I / j2„ i U 



d rp {dr^ V £ /P ) = J d rp (j J • (L2 - 41b) 

apart from irrelevant factors. Here uo is the planar vorticity e % ^divK Thus in the plane there is a 
denumarbly infinite set of particle relabeling constants. Again, time-independence of CK can be 
established directly from (1.2.41b) with the help of the continuity and Euler equations (1.1.16), 

(1.1.17) . 

The brackets of the above volume-integrated relabeling constants with the Euler Hamiltonian 
vanish, because they are time-independent. But this does not depend on the specific form of the 
Hamiltonian, since the fundamental brackets (1.2.16)-(1.2.21) give vanishing brackets for Cr^ with 
p and v. This is as it should be, because we have already remarked that the Eulerian variables 
do not respond to the relabeling transformations. Quantities whose bracket vanishes with all the 
elements of a (bracket) algebra, here p and v, are called Casimir invariants. So we see that the 
relabeling symmetry gives rise to Casimir invariants in the Euler formulation for fluids. This has a 
profound impact on the possibility of constructing a Lagrangian for Eulerian fluids. 

Finally, let us remark that even though the volume-preserving diffeomorphism transformations 
do not act on the Eulerian p and v, there remains in the formalism a related structure: The brackets 

(1.2.18) of currents (equivalently, momentum densities) present a local realization of the full (not 
merely volume-preserving) diffeomorphism algebra. For if we define for arbitrary f 

je = J drf(r).j(r), (1.2.42) 

then (1.2.18) implies 

{h,h} = h 2 , (1-2.43) 

where fi2 is the Lie bracket of fi and f2 

In = fidjfi - (1.2.44) 
(v) Lagrange function for Eulerian fluid mechanics 



While constructing the Euler Hamiltonian is straightforward, for example by transforming the 
Lagrange Hamiltonian, as in (1.2.10)-(1.2.13), an analogous construction for the Euler Lagrangian is 
problematic. First, the equations for the variables p and v are first order in time, so the Lagrangian 
should reproduce this. Second, time derivatives in the Lagrangian determine the canonical, bracket 
structure, which ultimately should reproduce (1.2.16)-(1.2.21). However, direct transcription of 
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Ll (1.2.1a) to Le, analogous to the passage from Hl in (1.2.2a) to He (1.2.13), would yield 
Le = f c/r(|pv 2 — V(p)), which contains no time derivatives. So something else must be done. 
Moreover, as we shall now explain, the presence of the Casimirs Cu) an d N poses obstructions to 
the construction of a Lagrangian, which must be overcome. 

Before proceeding, we present a Sidebar on the relation between Lagrangians and the canonical 
bracket structure. 
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A. Sidebar on canonical formalism determined by a Lagrangian 



(a) Easy case 

We begin with a Lagrangian that is first order in time. This entails no loss of generality because 
all second order Lagrangians can be converted to first order by the familiar Legendre transforma- 
tion that produces a Hamiltonian: H(p,q) = pq — L(q,q), where p = dL/dq. The equations of 
motion gotten by taking the Euler-Lagrange derivative with respect to p and q of the Lagrangian 
L(p,p; q, q) = pq — H(p, q) coincide with the "usual" equations of motion obtained by taking the q 
Euler-Lagrange derivative of L(q, q). [In fact L(p,p; q, q) does not depend on p.] Moreover, some La- 
grangians possess only a first-order formulation (for example, Lagrangians for Schrodinger or Dirac 
fields; also the Klein-Gordon Lagrangian in light-cone coordinates is first order in the light-cone 
"time" derivative). 

Denoting all variables by the generic symbol £ l , the most general first order Lagrangian is 

L = a,(Oe-#(0- (A.l) 

Note that although we shall ultimately be interested in fields defined on space-time, for present 
didactic purposes it suffices to consider variables that are functions only of time. The Euler- 
Lagrange equation that is implied by (A.l) reads 

h(oe = (a.2) 

where 

Ms) - -Qg Qg-- y A - 6 ) 

The first term in (A.l) determines the canonical 1-form: <2j(£)£ J dt = <2j(£)d£\ while gives the 
symplectic 2-form: da;(£) df = l/-^) dfdf . 

To set up a canonical formalism, we proceed directly. We do not make the frequently heard state- 
ment that "the canonical momenta dL/d^ 1 = a;(£) are constrained to depend on the coordinates 

and we do not embark on Dirac's method for constrained systems [7]. 

In fact, if the matrix possesses the inverse there are no constraints. Then (A.2) implies 

e = r'(0^ff^ (a.4) 

When one wants to express this equation of motion by bracketing with the Hamiltonian 
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one is led to postulating the fundamental bracket as 

{e,e'} = -r'(0- (a.6) 

The bracket between functions of £ is then defined by 

W(0,^tt)} = -^r^. (a.7) 

One verifies that (A.6), (A.7) satisfy the Jacobi identity by virtue of the Bianchi identity when / 
is given by (6.3.41). 



d d d 



(b) Difficult case 



When fij is singular, we may still proceed in the following manner [8]. Let us suppose that fij 
possesses N zero modes 

P\n)fij = 

n=l,...,A. (A.9) 
If we use a rank N projection operator P i 3 that satifies 

p i i p j h = p i k , lyj,^ c. (a. io) 

it is possible to find an inverse for on the projected subspace. Namely, the "inverse" f lJ is 
uniquely determined by 



Once P 3 is constructed, we define the Poisson bracket for functions of £ l by (A.7). 

It still remains to verify the Jacobi identity. An easy computation shows that (A.7) continues 
to satisfy that identity provided 

f,^ = 0. (A.12) 

Hence we use the brackets (A.7) only between functions that satisfy the admissability criterion 
(A.12). 
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(c) Obstructions to a canonical formalism 

Our problem in connection with Eulerian fluid mechanics is in fact the inverse of what has been 
summarized above. From (1.2.16), (1.2.19) and (1.2.20), we know the form of and that the Jacobi 
identity holds. We then wish to determine the inverse fy, and then a* from (6.3.41). Since we know 
the Hamiltonian from (1.2.13), construction of the Lagrangian (A.l) should follow immediately. 

However, an obstacle arises: Since there exist a Casimir invariants C(£) whose brackets with all 
the vanish, then 

= {C,C(O} = -f l ^C(O . (A.13) 

That is, f l i has zero modes g|jC(£), and the inverse to namely the symplectic 2-form /y, does 
not exist. In that case, something has to be done to neutralize the Casimirs. 

(d) Canonical transformations 

It is interesting, and will be later useful, to develop the theory further into a discussion of canonical 
transformations [9]. A canonical trasfromation is a transformation on the phase space coordinates 
£ l , given infinitesimally by a vector function, 

5C = -v\0, (A.14) 
which leaves the symplectic 2-form invariant. 

Sfij = v m d m f VJ + 0, v m f mj + dj v m f im = (A. 15a) 

[The exrpession in (A. 15a) is the Lie derivative of with respect to the vector field v m .] Use of 
the Bianchi identity (A. 8) allows casting (A. 15a) into 

di(v m f mj ) - <9> m /™) = 0. (A.15b) 

This shows that the quantity v m f mi may be presented as 

G(£) is called a "generator" for the vector field v m . Conversely, when f m i possesses an inverse, 
given any function G(£) on phase space we can define a vector field v m . Eq. (A. 15b) is the general 
requirement for v % to define a canonical transformation. Eq. (A. 16) is a necessary and sufficient 
condition, locally on phase space. If the phase space has non-trivial topology, one can have more 
general solutions to the condition (A. 15). 

The change of a function F(£) on phase space, due to a canonical transformation, is given by 

dF 

f)C f)F 
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where (A. 6) and (A. 7) and (A. 16) have been used. This shows that indeed G generates the in- 
finitesimal canonical transformation by Poisson bracketing. 

Note that the bracket (A. 17) has been established without the inverse 2-form f l K This enjoys an 
advantage over explicit evaluations relying on the methods presented above in (a), and especially in 
(b), where a projected inverse is employed. Moreover we see that the admissability criterion (A. 12) 
is automatically satisfied by any generator G that solves (A. 16). This follows immediately from 
(A.9). 

To construct a Lagrangian for Euler fluids, which leads to the correct equations of motion 
(1.1.16), (1.1.17) and brackets (1.2.16)-(1.2.21), we must neutralize the Casimirs. In three and one 
spatial dimensions, we must neutralize the velocity Chern-Simons terms (1.2.59) and (1.2.40), and 
also the total number N. In two dimensions the Casimirs which must be neutralized comprise the 
infinite tower (1.2.41), as well as N. 

In three dimensions this is achieved in the following manner, based on an idea of C.C. Lin [10]. 
We use the Clebsch parameterization for the vector field v [11]. Any three-dimensional vector, 
which involves three functions, can be presented as 

v = V0 + «VA (1.2.45) 

with three suitably chosen scalar functions 9, a, and (3. This is called the Clebsch parameterization, 
and (a, (3) are called Gaussian potentials. In this parameterization, the vorticity reads 



and the Lagrangian is taken as 



w = Vax V/3, (1.2.46) 



J d 3 rp(8 + a(5) - H\ v=ve+a ^ (1.2.47) 



with v in H expressed as in (1.2.45). Thus the canonically conjugate pairs are (p, 9) and (pa, f3) re- 
placing p and v. The phase space (p, 9, a, (3) is 4- dimensional, corresponding to the four observables 
p and v, and a straightforward calculation shows that the Poisson brackets (1.2.16), (1.2.19)-(1.2.21) 
are reproduced with v constructed by (1.2.45). 

But how has the obstacle presented by the Casimirs been overcome? Let us observe that in the 
Clebsch parameterization C% is given by 

C 3 = J d 3 re i]k 8i9 d j( y d k (3, (1.2.48) 

which is just a surface integral 

C 3 = I dS-(0w). (1.2.49) 
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In this form, C3 has no bulk contribution, and presents no obstacle to constructing a symplectic 
2-form and a canonical 1-form in terms of (p, 9, a, (3), which are defined in the bulk, that is, for all 
finite r. Moreover, the brackets with N are no longer universally vanishing. Specifically with the 
new dynamical variable 9 we find 

{N t 6} = -1. (1.2.50) 
Note that the response of 9 to a finite boost with velocity u 

u 2 

9(t, r) -> 9(t, r - ut) + u • r - —t, (1.2.51) 

contains the 1-cocycle u ■ r — ^-t, as is familiar from representations of the Galileo group. 

A "derivation" of (1.2.47) can be given, based on ideas of Lin and earlier ones of Eckart [12]. We 
begin with Le, the Euler transcription of the Lagrange Lagrangian, mentioned earlier and containing 
no time derivatives. This is supplemented with constraints, enforced by Lagrange multipliers that 
ensure various continuity equations. 

d 3 r Qpv 2 - V{p) +9(p + V ■ (vp)) - pa{/3 + v • V/3) J (1.2.52) 

The first two terms in the integrand reproduce Le] the first multiplier, 9, enforces the mat- 
ter/current continuity equation. The second continuity equation for (3, enforced by the multiplier 
pa, is physically obscure. (Lin argues that it reflects the conservation of "initial data" for fluid 
motion.) Varying v and eliminating it gives (1.2.47). 

Let us carry out all the variations and re-derive the two Eulerian equations. For greater gener- 
ality, useful for relativistic kinematics, we shall take an arbitrary kinetic energy pT(v), rather than 
the above non relativistic case T(v) = |v 2 , and define the momentum p be the derivative of T. 

Varying 9 gives the continuity equation (1.1.16); varying a gives 

/3 + v-V/? = 0, (1.2.54a) 
and varying (3 gives a similar continuity equation for a. 

d + v-Va = (1.2.54b) 

Next we vary v to find 

pp - p(V9 + aVf3) = 0. (1.2.55) 

Thus in the general case it is p (rather than v) that is given by the Clebsch parameterization. With 
(1.2.55) the Lagrangian (1.2.52) is rewritten, apart from a total time derivative, as 

L = J d 3 r(^pT(v) -V{p) -p(9 + a$) -pv-pj. (1.2.56) 
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Let h(p) be the Lagrange transform of T(v) 



h(p) = p-v-T(v), (1.2.57a) 
^(p) 



dp 

and L becomes 



1.2. 57b) 



L = J d 3 r(-p(9 + ap)-ph(p)-V(p)\ (1.2.58) 
The remaining variable to vary is p. 

6 + a(3 + h(p) + V'(p) = (1.2.59a) 
Differentiating with respect to r l converts this to 

9 .0 + 9 . a/ j + a ^ = " ^'(P)^P- (1.2.59b) 

But according to (1.2.55) ^p- 5 = djp 1 + diOtdjf3 — djOtdi(3. Thus the first term on the right in (1.2.59b) 



is 



Op 1 

—v l djp' 1 — diav j djP + v j dj<ydiP = —v j —rdjV k + diaf3 — adiP, 



where we have used (1.2.74) and (1.2.55). Rearranging (1.2.59b) leads to 

3$ + adi(3 + adi/3 + tfdrf = -V"(p)d iP , (1.2.60a) 

f + v'djp* = -V"(p)d iP . (1.2.60b) 

With Newtonian kinematics p l = v l , and Euler's equation is regained. With arbitrary kinematics 
p l = dT d ^ ■ Also, since 5p l = ^-Sv^ = t^Sv^ where r*- 7 = ; the above becomes 



or 



r ik (v k 4 



id jV k ) = -V"{p)d iP , (1.2.61a) 



or 



v* + vidjV* = -(T~ 1 Y ] V"(p)d j p, (1.2.61b) 

whenever the inverse to r*- 7 exists. For Newtonian kinematics r 1 - 7 = 6 13 . 

Note that the free Euler equation, (1.2.61b) with V" = 0, is insensitive to the form of the kinetic 
term T(v) (provided t 4 - 7 possess an inverse), and can be solved together with the continuity equation 
(1.1.16). The general solution in the non- interacting case is presented in terms of posited initial 
data. 



p(t = 0,r) = p (r) 
v(t = 0,r) = v (r) 



(1.2.62) 
(1.2.63) 
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Define the quantity x(t, r) by the equation 

r = *v (x(t,r))+ X (t,r). (1.2.64) 
Then (1.1.16) and (1.1.17) (with VP = 0) are solved by 

p(t,r) = Po ( x )\det^-\, (1.2.65) 
v(t,r) = v (x). (1-2.66) 

This result is verified by differentiation; alternatively it may be derived from (1.1.13), (1.1.14), with 
X(£, x) taken in the absense of forces to be a linear function of t. 

In one dimension, we parameterize the velocity as a derivative of a potential 9, 

v = 9', (1.2.67) 

and the phase space consist of (p,9). The Casimir (1.2.40) again becomes a surface term (only 
endpoints contribute) and is neutralized in the bulk. The two variables are the conjugate pair (p, 9), 
which capture the two degrees of freedom (p, v). Of course lineal fluids possesses no vorticity, so the 
velocity bracket (1.2.20) vanishes, while (1.2.19) is verified. [In an alternative approach to lineal 
fluids, we replace 9(r) by | f dr'e(r — r')v(r') and the canonical 1-form f dr9p is — ~ J drdr' p(r)e(r — 
r ; )v(r'), where e is the ±1 step function. Evidently this leads to a spatially non-local, but otherwise 
completely satisfactory canonical formulation for fluids on a line.] 

Two dimensions presents the additional problem that the number of physical variables is three: 
(p, v), but an odd number cannot form a symplectic structure. At the same time there is an 
infinite number of Casimirs, (1.2.41). One may then conclude, heuristically, that it should be 
possible to neutralize an infinite number of non-local Casimirs by suppressing one local degree of 
freedom, thereby decreasing the effective variables from three to two, an even number with which 
one can build a symplectic structure. But it is not known how to effect this suppression; rather the 
Lin/Clebsch method is adopted, which increases the degrees of freedom to four and the Lagrangian 
(1.2.47) is used in two dimensions as well. 

Finally note that the Lagrangian in (1.2.47), apart from a total time derivative, can also be 
written as 

L = - J dr[p{9 + a/3) + pv • ( V0 + aV/3)] + L E 

= - J drfid^ + ad^) + L E . (1.2.68) 

Although we are dealing with nonrelativistic dynamics, we have used covariant notation for the 
canonical 1-form, with = (cp, pv) and = (^, V). This formulation becomes our starting 
point for the relativistic generalization, discussed in Section 1.4 below. (That is why we have intro- 
duced the velocity of light c in the above definitions; of course it disappears in the nonrelativistic 
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theory, where it has no role.) 



B. Sidebar on Clebsch parameterization and the Chern-Simons term 

We elaborate on the Clebsch parameterization for a vector field [11], which was presented for 
the velocity vector in (1.2.45). Here we shall use the notation of electromagnetism and discuss the 
Clebsch parameterization of a vector potential A, which also leads to the magnetic field B = V x A. 
(Of course the same observations apply when the vector in question is the velocity field v, with 

V x v giving the vorticity.) 

The familiar parameterization of a three-component vector employs a scalar function 9 (the 
"gauge" or "longitudinal" part) and a two-component transverse vector A^: A = V# + At, 

V • At = 0. This decomposition is unique and invertible (on a space with simple topology). In 
contrast, the Clebsch parameterization uses three scalar functions, 6, a, and /3, 

A = V6 + aV(3, (B.l) 

which are not uniquely determined by A (see below). The associated magnetic field reads 

B = V x A = Va x V/3. (B.2) 

Repeating the above in form notation, the 1-form A = A4 dr % is presented as 

A = d6+adf3, (B.3) 

and the 2-form is 

dA = dad(3. (B.4) 

Darboux's theorem [13] ensures that the Clebsch parameterization is attainable locally in space [in 
the form (B.3)]. Additionally, an explicit construction of a, (3, and 9 can be given by the following 
procedure [14]. 

Solve the equations 

dx dy dz 

TT = TT = Tp (B.5a) 

£>x &y &Z 

which may also be presented as 

s ijk dr j B k = 0. (B.5b) 

Solutions of these relations define two surfaces, called "magnetic surfaces" , that are given by equa- 
tions of the form 

S n (r) = constant (n = 1,2). (B.6) 
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It follows from (B.5) that these also satisfy 

B-VS n = (n = l,2), (B.7) 

that is, the normals to S n are orthogonal to B, or B is parallel to the tangent of S n . The intersection 
of the two surfaces forms the so-called "magnetic lines" , that is, loci that solve the dynamical system 

^ = B(r(r)), (B.8) 

where r is an evolution parameter. Finally, the Gaussian potentials a and (3 are constructed as 
functions of r only through a dependence on the magnetic surfaces, 

a(r) = a(S n (r)), 

P(r) = P(S n (r)), (B.9) 

so that 

V« x V/3 = (VSi x VS 2 )e mn -^^-. (B.10) 

Evidently as a consequence of (B.7), Vax V/3 in (B.10) is parallel to B, and because B is divergence- 
free a and (3 can be adjusted so that the norm of V« x V/3 coincides with |B|. Once a and (3 have 
been fixed in this way, 9 can easily be computed from A — «V/3. 

Neither the individual magnetic surfaces nor the Gauss potentials are unique. [By viewing 
A as a canonical 1-form, it is clear that the expression (B.3) retains its form after a canonical 
transformation of a,/3.} One may therefore require that the Gaussian potentials a and (3 simply 
coincide with the two magnetic surfaces: a — S\, (3 = 5*2 . Nevertheless, for a given A and B it may 
not be possible to solve (B.5) explicitly. 

The Chern-Simons integrand A ■ B becomes in the Clebsch parameterization 

A • B = VO ■ (V« x V/3) = V • (0 B) = B • V0. (B.ll) 

Thus having identified the Gauss potentials a and f3 with the two magnetic surfaces, we deduce 
from (B.7) and (B.ll) three equations for the three functions (8, a, (3) that comprise the Clebsch 
parameterization. 

B • V« = B • V/3 = 

B • V# = Chern-Simons density A ■ B (B.12) 

Eq. (B.ll) also shows that in the Clebsch parameterization the Chern-Simons density becomes 
a total derivative. 

A-B = V-(0B) (B.13) 
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This does not mean that the Clebsch parameterization is unavailable when the Chern-Simons inte- 
gral over all space is nonzero. Rather for a nonvanishing integral and well-behaved B field, one must 
conclude that the Clebsch function 9 is singular either in the finite volume of the integration region 
or on the surface at infinity bounding the integration domain. Then the Chern-Simons volume 
integral over (Q) becomes a surface integral on the surfaces (dQ) bounding the singularities. 



d d rA-B= / dS-(0B) (B.14) 
n Jon 

Eq. (B.14) shows that the Chern-Simons integral measures the magnetic flux, modulated by 9 and 
passing through the surfaces that surround the singularities of 9. 

The following explicit example illustrates the above points. 

Consider the vector potential whose spherical components are given by 

A r = (cosO)a'(r), 

A® = — (sin Q)- sin a(r), 
r 

A$ = -(sin 9)-(l -cosa(r)). (B.15) 

(r, and 0, $ denote the conventional radial coordinate and the polar, azimuthal angles.) The 
function a(r) is taken to vanish at the origin, and to behave as 2itv at infinity [y integer or half- 
integer). The corresponding magnetic field reads 

B r = — 2(cos6) — (l — cos ct(r)), 

Bq = (sin 6) -a'(r) sin a(r), 
r 

B$ = (sine)-a'(r)(l - cosa(r)), (B.16) 

and the Chern-Simons integral - also called the "magnetic helicity" in the electrodynamical context 
- is quantized (in multiples of 167T 2 ) by the behavior of a(r) at infinity. 

f f°° d 

/ d 3 rA ■ B = — 8ir / dr — (a(r) — sin a(r)) 
J Jo dr v 

= -167rV (B.17) 

In spite of the nonvanishing magnetic helicity, a Clebsch parameterization for (B.15) is readily 
constructed. In form notation, it reads 

A = d(-2$) +2(l - (sin 2 ^) sin 2 e) d($ + tan _1 (tan|)cose). (B.18) 

The magnetic surfaces can be taken from formula (B.18) to coincide with the Gauss potentials. 

51 = 2(l — (sin 2 — ) sin 2 = constant 

52 = $ + tan -1 (tan — ) cos G = constant (B.19) 

2 
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The magnetic lines are determined by the intersection of Si and S2. 

cos — = e cos($ — 0o ) 



sin0 — \ \ 2 1 2 f. , , (B.20) 

' 1 - e 2 cos 2 ($ - O ) 

where e and 0o are constants. The potential 9 = —2$ is multivalued. Consequently the "surface" 
integral determining the Chern-Simons term reads 



00 

d 3 rA • B = I d a rV ■ (-2&B) = -4vr / r dr I dB B, 



Jo *= 27r 



(B.2i; 



That is, the magnetic helicity is the flux of the toroidal magnetic field through the positive-x (x, z) 
half-plane. 

Finally we remark on a subtle property of the Clebsch decomposition when used in variational 
calculations [15]. Consider an "action" of the form 

/ = J (B) + I J d 3 rA ■ B. (B.22) 



Variation of A gives 

SI = J d 3 r(V x B + fiB) ■ 5 A, (B.23) 

h 

where £>(r) = ^3 ^ ■ Demanding that I be stationary against variations of A requires the vanishing 
of the term in parentheses, which is transverse, since the transverse part of the variation 5 A is 
arbitrary. 

V x B + fiB = 0, (B.24) 

Now let us examine the same problem in the Clebsch parameterization. The Chern-Simons contri- 
bution to (B.22) reads, according to (B.ll) 

^ / d 3 rV • (0B) = ^ [dS-ffB. (B.25) 



2 7 2 „ 



In the gauge = (B.25) vanishes, and in any gauge it has no bulk contribution, so its variation 
will never produce the second left-hand term in (B.24). So how is (B.24) regained? 

Returning to (B.22), we accept the fact that the variation of the last term vanishes, while the 
variation of the first leaves 

61 = f d 3 r(V x B) ■ (V/3 5a - Va 5(3). (B.26) 
h 
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Since 6a and 6(3 are arbitrary, I is stationary provided 



(V x B) ■ V/3 = 0, 
(V x B) ■ V« = 0. 



(B.27a) 
(B.27b) 



We obtain two equations, which imply by (B.12), 



V x B + p(r)B = 0. 



(B.28) 



Transversality of V x B and B further requires that B ■ V/x(r) = 0, i.e., a non-vanishing V/i(r) 
(= non-constant /i) must lie in the (Vq, V/5) plane. 

Since (B.28) is weaker than the parameteriztion-independent (B.24), we conclude that the 
Clebsch paramterization is somewhat incomplete, when used in variational calculations that ignore 
surface terms. This is similar to the fact that in the Clebsch parameterization gauge potentials, 
which carry a non- vanishing Chern-Simons term (= velocities with non- vanishing vortex helicity), 
encode the non- vanishing value in a surface term. 

In Section 7.4, we shall present a different method, based on group theory, for obtaining the 
Clebsch parameterization. This approach is then generalized to non-Abelian vector fields. 



1.3 Irrotational fluids 

The simplification found for 1-dimensional fluids, (1.2.67), presenting the velocity as a derivative 
of a velocity potential 9, can be extended to two-and three-dimensional fluids that are irrotational: 
the vorticity vanishes. 



The Clebsch parameterization (1.2.45) holds trivially; the potentials a,f3, the Casimirs Cu) van- 
ish, and N obeys (1.2.50). This removes the obstruction to a canonical formalism. The Lagrangian 
(1.2.47) becomes 



and this can be derived by the Eckart procedure as in (1.2.52) - (1.2.58), where now only the 
continuity equation is enforced by the Lagrange multiplier 9, and the Gaussian potentials a and (3 
are omitted. 

The Euler equation (1.1.17), with the pressure P expressed in terms of the enthalpy V'(p) as in 
(1.2.14) may be integrated once to give the Bernoulli equation. 



u) = V x v = 



v = V6 



(1.3.1) 
(1.3.2) 




(1.3.3) 



V\p) 



(1.3.4) 
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This also follows by bracketing 9 with H in (1.3.3), since now there is only one non- vanishing and 
canonical bracket. 

{e(v),p(v')} = 8(v-v') (1.3.5) 

Observe that the Bernoulli equation (1.3.3) allows for the possibility of expressing p in terms of 
9 + \{yef, through the inverse to V . One can then eliminate p in the continuity equation, leaving 
a single (non-linear) equation for 9. 



1.4 Relativistic fluids 

We discuss the formalism for describing relativistic fluids; only the Euler approach is treated. Usu- 
ally the dynamics of a relativistic fluid is presented in terms of the energy- momentum tensor, 9^, 
and the equations of motion are just the conservation equations 8^9^ = 0. This is analogous 
to the nonrelativistic situation mentioned previously, where the nonrelativistic energy momentum 
complex T^ u encapsulates the equations of motion for a nonrelativistic fluid. [We denote the rela- 
tivistic energy-momentum tensor by 9^ u , to distinguish it from the nonrelativistic T^ u introduced 
in (1.2.24)-(1.2.28). The limiting relation between the two is given below.] But we shall begin with 
a Lagrange density. 

Inspired by the suggestive formula (1.2.68), we consider 

C = -fa,-f( v ^%). (1.4.1) 

Here j M is the current Lorentz vector j M = (cp, j) [16]. The a M comprise a set of auxiliary variables; 
in the relativistic analog of irrotational fluids we take a M = 8^9, more generally 

a M = 8^0 + ad^, (1.4.2) 

so that the Chern-Simons density of cij is a total derivative [compare (1.2.48), (1.2.49)]. The 
function / depends on the Lorentz invariant j^j^ = c 2 p 2 — j 2 and encodes the specific dynamics 
(equation of state). 

The energy momentum tensor for L is 

V = -g^C + _^ = /'( > /^). (1.4.3) 

V J Ja 

[One way to derive (1.4.3) from (1.4.1) is to embed that expression in an external metric tensor g^, 
which is then varied; in the variation and a M are taken to be metric-independent and j M = g^ u j u -] 
Furthermore, varying j M in (1.4.1) shows that 



so that (1.4.3) becomes 



VJ 3a 

d^v = -g»v[nf'(n) - f(n)} + u^nf '(n). (1.4.5) 
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We have introduced the proper velocity by factoring j M , as suggested by Eckart. 

jfi = nu^ u^u^ = 1 (1.4.6) 

One sees that n is proportional to the proper density and 1/n is proportional to the specific volume. 

Eq. (1.4.5) identifies the proper energy density e and the pressure P (which coincides with C) 
through the conventional formula [17] 

V = -9^P + u^u v (P + e). (1.4.7) 

Therefore, in our case 

e = f(n) (1.4.8) 

p = nf{n) - f(n) . (1.4.9) 

The thermodynamic relation involving entropy S reads 

Pd(-) +d(-) ocd£, (1.4.10) 



where the proportionality constant is determined by the temperature. With (1.4.8) and (1.4.9) the 
left side of (1.4.10) vanishes and we verify that entropy is constant, that is, we are dealing with an 
isentropic system, as has been stated in the very beginning. 

For the free system, the pressure vanishes, so we choose f(n) = cn. 



Co = - C y/jVjp (1.4.11) 

Other forms for / give rise to relativistic fluid mechanics with other equations of state. 
Taking the divergence or 9^ u in (1.4.5) leaves 

d^e^ = -nf"{n)d v n + nu^{u v f{n)) + d„(nu>")u v f'(n). (1.4.12) 

The first two terms on the right are orthogonal to u u , the last one is parallel. So the vanishing of 
the divergence of 6^ u implies the continuity equation. 

= 9,(nu") = d, L f (1.4.13) 

The vanishing of the remaining components is equivalent to 

uHd^ - d^u v )f\n) + {g Vil - u v u^nf"{n) = 0. (1.4.14) 

This is the relativistic Euler equation. 

The same two equations follow from the Lagrange density (1.4.1). Variation of 9 and the Gauss 
potentials a and (3 gives the current continuity equation (1.4.13) and equations satisfied by a and 

P 

u n d„a = = d„p (1.4.15) 
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Variation of j M gives (1.4.4), whose curl reads 

(d v - u y )f{n) + (up d u n-u u n)f"(n) = <9 M a u - d u = d^ad u [3 - d v ad^. (1.4.16) 

Contracting this with w M makes the right side vanish by virtue of (1.4.15) and the left side coincides 
with (1.4.14). 

It is especially intriguing to notice that 9^ u is symmetric but T^ v is not. To make the connection 
we recall that u M = 1/ y/l — v 2 /c 2 (l, v/c), we observe that n = a/ p 2 c 2 — j 2 , set j = vp and conclude 
that n = pcyl — v 2 /c 2 ~ pc — (pv 2 /2c). Also f(n) is chosen to be cn + V(n/c), and thus P = 
nf'{n) - f{n) = (n/c)V'(n/c) - V{n/c). It follows that 

Q00 ^ nc- (vW* + y ^ pc 2 - pv 2 /2 + 
1 — v^/r 1 — v z /cr 

2 

« pc 2 + ^- + \/(p) = pc 2 + T°° . (1.4.17) 



Thus, apart from the relativistic "rest energy" pc , # 00 passes to T°°. The relativistic energy flux is 
c#° (because ^f^ = + 6^'°) 



(!■' 



; pc 2 -pv 2 /2 + pF'(p) 



v 2 /c 2 V c / 1 — v 2 /c 2 

« j j c 2 + p^'(v 2 /2 + 1/'(p)) = fc + T jo (1.4.18) 

Again, apart from the 0(c 2 ) current, associated with the 0(c 2 ) rest energy in T JO is obtained 
in the limit. The momentum density is 6 m /c (because has dimension of energy density). 

9 oi / c = - - (nc + -V') « pv l = V 1 (1.4.19) 

1 — \ l jc z V c / 

Finally, the momentum flux is obtained directly from 9 % K 

n T r . \ u V ( n , 



= s» (-v'-v) + - — -[nc + -v 

« & j {pV\p) - V(p)) + uV> = T ij ' (1.4.20) 
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2 SPECIFIC MODELS (d ^ 1) 

We now examine irrotational models, both relativistic and nonrelativistic, for which we shall specify 
an explicit force law and discuss further properties. 

2.1 Two models, mostly in spatial dimensions d ^ 1. 
(i) Galileo-invariant nonrelativistic model: Chaplygin gas 

Recall that the nonrelativistic Lagrangian for irrotational motion reads 

£ Galilc ° = / dr(6p - - V(pj), (2.1.1) 

where V0 = v. The Hamiltonian density 7i is composed of the last two terms beyond the canonical 
1-form J dr9p 

H = J dri K p ( —^- + V{p)^ = J drH. (2.1.2) 

Varying (2.1.1) with respect to p produces the Bernoulli equation (1.3.4 ). Various expressions for 
V appear in the literature. V(p) oc p n is a popular choice, appropriate for the adiabatic equation 
of state. We shall be specifically interested in the Chaplygin gas [5] . 

V(p) = \/p, A>0 (2.1.3) 

According to what we said before, the Chaplygin gas has enthalpy V = —X/p 2 , negative pressure 
P = — 2A/p, and speed of sound s = V2X/p (hence A > 0). 

Chaplygin introduced his equation of state as a mathematical approximation to the physically 
relevant adiabatic expressions with n > 0. (Constants are arranged so that the Chaplygin formula 
is tangent at one point to the adiabatic profile.) Also it was realized that certain deformable solids 
can be described by the Chaplygin equation of state. These days negative pressure is recognized as 
a possible physical effect: exchange forces in atoms give rise to negative pressure; stripe states in the 
quantum Hall effect may be a consequence of negative pressure; the recently discovered cosmological 
constant may be exerting negative pressure on the cosmos, thereby accelerating expansion. 

For any form of V, the model possesses the Galileo symmetry, discussed previously as appropriate 
to nonrelativistic dynamics. There are a total of |(d + l)(d + 2) Galileo generators in d space plus 
one time dimensions. Together with the central term, N, we have a total of \{d + l)(d + 2) + 1 
generators. 

A useful consequence of symmetry transformations is that they map solutions of the equations 
of motion into new solutions. Of course, "new" solutions produced by Galileo transformations are 
trivially related to the old ones: they are simply shifted, boosted or rotated. 
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But we shall now turn to the specific Chaplygin gas model, with V(p) = A/p, which possesses 
additional and unexpected symmetries. 

The Chaplygin gas action and consequent Bernoulli equation for the Chaplygin gas in (d, 1) 
space-time read 

^Chaplygin = j dt j dr (Bp - - -) (2.1.4) 

This model possesses further space-time symmetries beyond those of the Galileo group [18]. First, 
there is a one-parameter (u, dimensionless) time rescaling transformation 

t^T = e UJ t, (2.1.6) 

under which the fields transform as 

6(t,T)^6 u {t,T)=e?6(T,T), (2.1.7) 
p(t,r)^p u (t,r)=e- u p(T t T). (2.1.8) 

Second, in d spatial dimensions, there is a vectorial, <i-parameter (w, dimension inverse velocity) 
space-time mixing transformation. 

t -> T(t,r) = t + u • r + §u; 2 0(T,R) (2.1.9) 
r -> R(t,r) = r + w#(T,R) (2.1.10) 

Note that the transformation law for the coordinates involves the 9 field itself. Under this transfor- 
mation, the fields transform according to 

0(t,r)->0 w (t,r) = 0(T,R), (2.1.11) 
p(t, r) - p^t, r) = p(T, R)-L (2.1.12) 

with J the Jacobian of the transformation linking (T, R) — > (£, r). 

(l - w • V0(T, R) - y ^ T > R )) 1 ^ 2 - 1 - 13 ) 



J = det 



/<9T 






dri 






\ dt 


dri / 



(The time and space derivatives in the last element are with respect to t and r.) One can tell the 
complete story for these transformations: The action (2.1.4) is invariant; Noether's theorem gives 
the conserved quantities, which for the time rescaling is 

S = tH- J dr P 9 (time rescaling), (2.1.14) 
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while for the space-time mixing one finds 

G = / dr (r7Y — 9T*) (space-time mixing). (2.1.15) 



The time independence of S and G can be verified with the help of the equations of motion 
(continuity and Bernoulli) [19]. Poisson bracketing the fields 9 and p with S and G generates the 
appropriate infinitesimal transformation on the fields. Note that unlike the Galileo constants of 
motion, the new constants of motion cannot be locally expressed in terms of v: their integrands 
depends locally on 9 = J? • v. 

So now the total number of generators is the sum of the previous + + 2) + 1 with 1 + d 
additional ones. 

\{d + + 2) + 1 + 1 + d = \(d + 2){d + 3) (2.1.16) 

When one computes the Poisson brackets of all these with each other one finds the Poincare Lie 
algebra in one higher spatial dimension, that is, in (d + 1, l)-dimensional space-time, where the 
Poincare group possesses |(d + 2)(d + 3) generators. Moreover, one verifies that (t,8,r) transform 
linearly as a (d + 2) Lorentz vector in light-cone components, with t being the "+" component and 
9 the "-" component. [20] 

Presently, we shall use these additional symmetries to generate new solutions from old ones, 
but, in contrast to the Galileo transformations, the new solutions will be nontrivially linked to the 
former ones. Note that the additional symmetry holds even in the free theory. 

Before proceeding, let us observe that p may be eliminated by using the Bernoulli equation to 
express it in terms of 9. In this way, one is led to the following p-independent action for 9 in the 
Chaplygin gas problem: 



jChaplygin = ^ ^ J ^ Q + ^_ZL (2.1.17) 

Although this operation is possible only in the interacting case, the interaction strength disappears 
from the equations of motion. 

!_^ + V.^^ = (2.1,8) 



A merely serves as an overall factor in the action. 

The action (2.1.17) looks unfamiliar; yet it is Galileo invariant. [The combination 9 + |(V6*) 2 
responds to Galileo transformations without a 1-cocycle; see (1.2.71).] Also (2.1.17) possesses the 
additional symmetries described above, with 9 transforming according to the previously recorded 
equations. 

Let us discuss some solutions. For example, the free theory is solved by 

9(t,r) = ^ (2.1.19) 
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which corresponds to the velocity 

v(t,r) = £. (2.1.20) 

Galileo transformations generalize this in an obvious manner into a set of displaced, rotated and 
boosted solutions. Performing on the above formula for 9 the new transformations of time-rescaling 
and space-time mixing (2.1.6), (2.1.7), (2.1.9)-(2.1.11), we find that the solution is invariant. 

We can find a solution similar to (2.1.19) in the interacting case, for d > 1, which we henceforth 
assume. (The d — 1 case will be separately discussed in Section 3.) One verifies that a solution is 

r 2 l0\ \i\ 1 9 X 1 

0{t, r) = iVf r ) = V -j( d - !)- = V 4 ~ t 2 ' 1 - 21 ) 

2 (a — l)t V a r V a v 

r /o \ 



Note that the speed of sound, s = \/2X/p = \fdv, exceeds v. Again this solution can be translated, 
rotated, and boosted. Moreover, the solution is time-rescaling-invariant. However, the space-time 
mixing transformation (2.1.9)-(2.1.13) produces a wholly different solution. This is best shown 
graphically, where the d = 2 case is exhibited (see Figure) [21]. 

Another interesting solution, which is essentially one- dimensional (lineal), even though it exists 
in arbitrary spatial dimension, is given by 

9{t,r) = 0(n-r) + u-r-it(u 2 -(n-u) 2 ). 

v(t,r) = u + n9'(n-r) (2.1.23) 

Here n is a spatial unit vector, u is an arbitrary vector with dimension of velocity, while G is an 
arbitrary function with static argument. The corresponding charge density is time-independent. 

/ s 2A 

P(t, r = —— rr = r 2.1.24 

n ■ u + G'fn • r) n • v v ' 



The current is static and divergenceless. 



j(t,r) = v / 2A(— —A (2.1.25) 



n ■ v 



The sound speed s — v2A/p = |n • v| is smaller than v. 

Finally, we record a planar static solution to (2.1.18), which depends on two orthogonal unit 
vectors ni, and n2 [22]. 

9(t,r) = e(fii ■ r/n 2 ■ r). 
(n^-n^T)) 
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10 6 

The original density p(t, r) oc — 

r 



(in two spatial dimensions, r = (x,y)) 




10 10 

The transformed density p(t,r). 
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This gives the density 

p = - 2.1.27 

v 

Now the sound speed coincides with v. 

(ii) Lorentz-invariant relativistic model: Born-Infeld model 

We now turn to a Lorentz-invariant generalization of our Galileo-invariant Chaplygin model in 
(d, l)-dimensional space-time. We already know from (1.2.73)-(1.2.78) how to construct the free 
Lagrangian with a relativistic kinetic energy. 



T(v) = -cVl-v 2 /c 2 (2.1.28) 

Mass has been scaled to unity, and we retain the velocity of light c to keep track of the nonrelativistic 
c —> oo limit. Evidently, the momentum is 

dT(v) v 
p = = . 2.1.29 

<9v - v 2 /c 2 

Thus the free relativistic Lagrangian, with current conservation enforced by the Lagrange multi- 
plier 9, reads [compare (1.2.72), with Gaussian potentials, a and j3 omitted] 

£Lorcntz = J py /\ - v 2 / (? + B{f> +V ■ (vp))} . (2.1.30) 

This may be presented in a Lorentz-covariant form in terms of a current four- vector = (cp,vp). 
^Lorentz Q f e q Ua ti on (2.1.30) is thus equivalent to [compare (1.4.11)] 

^Lorentz = J fr^fQ^ _ C y/ji%) . (2.1.31) 

Eliminating v in (2.1.30), we find that p is irrotational, 

9T v V9 , 

p = — = - — = V0, v = =, 2.1.32 

and the free Lorentz-invariant Lagrangian reads 

L^ ovcntz = j 'dr(ep- pc 2 ^lTJyW~R)- (2.1.33) 

To find L^ alilco [(2.1.1) with V = 0] as the nonrelativistic limit of L^ orcntz in (2.1.33), a nonrel- 
ativistic 9 variable must be extracted from its relativistic counterpart. Calling the former #nr and 
the latter, which occurs in (2.1.33), 6*r, we define 

9 R = -c 2 t + 9 NR . (2.1.34) 
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It then follows that apart from a total time derivative 

^Lorcntz > ^Galileo _ (2.1.35) 

c— >oo 

Next, one wants to include interactions. While there are many ways to build Lorentz-invariant 
interactions, we seek an expression that reduces to the Chaplygin gas in the nonrelativistic limit. 
Thus, we choose 

L~ = J d r (0 P - VT^V^NW) , (2-1.36) 

where a is the interaction strength [23]. (The reason for the nomenclature will emerge presently.) 
We see from (2.1.4) that, as c — > oo, 

r Born-Infeld k r Chaplygin /r> i orA 

L >Ll X=a 2 /2 ■ {Z.L.6() 

C-^OO I 

[Again #nr is extracted from 6*r as in (2.1.34) and a total time derivative is ignored.] 

Although it perhaps is not obvious, (2.1.36) defines a Poincare-invariant theory, and this will 
be explicitly demonstrated below. Therefore, /Bam-infeid p OS s esses Poincare symmetries in (d, 1) 
space-time, with a total of \{d-\- l)(d + 2) + 1 generators, where the last " + 1" refers to the total 
number N = J dr p. 

When a = 0, the model is free and elementary. It was demonstrated previously [eqs. (1.2.73)- 
(1.2.81)] that the free equations of motion are precisely the same as in the nonrelativistic free 
model, so the complete solution (1.2.84)-(1.2.86) works here as well. For a ^ 0, in the presence of 
interactions, one can eliminate p as before, and one is left with a Lagrangian just for the 9 field. It 
reads 



L Born-Infeld = _ fl / dr I # _ (^2, (2.1.38) 



This is a Born-Infeld-type theory for a scalar field 9; its Poincare invariance is manifest, and 
again, the elimination of p is only possible with nonvanishing a, which however disappears from the 
dynamics, serving merely to normalize the Lagrangian. 

Although manifestly Lorentz covariant, the Lagrangian (2.1.38) is not of the form (1.4.1). To 
achieve that expression, we choose f(n) = c\/ a 2 + n 2 , so that the pressure P becomes 

P= , aC . (2.1.39) 

V a 2 + n 2 

This leads to the Chaplygin pressure as c — > oo with A = a 2 /2. Now the Lagrange density C = P 
coincides with (1.4.1 ) when is written as 
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The equations of motion that follow from (2.1.36) read 



P+V-lV9 t / ff±|l |=0, (2.1.41) 



The density p can be evaluated in terms of 9 from (2.1.42); then (2.1.41) becomes 

d a ( ,„ * gag)=0, (2-1.43) 



V^ 2 " (P? 

which also follows from (2.1.38). After #nr is extracted from Or as in (2.1.34), we see that in the 
nonrelativistic limit L* OIn - lnield (2.1.36) or (2.1.38) becomes L^ haplygin of (2.1.4) or (2.1.17), 

L Born-Infcld > L^jf n , (2 . 1 .44) 

and the equations of motion (2.1.41)-(2.1.43) reduce to (2.1.5) and (2.1.18). 

In view of all the similarities to the nonrelativistic Chaplygin gas, it comes as no surprise that 
the relativistic Born-Infeld theory possesses additional symmetries. These additional symmetry 
transformations, which leave (2.1.36) or (2.1.38) invariant, involve a one-parameter (u, dimesionless) 
reparameterization of time, and a c/-parameter (w, dimension velocity) vectorial reparameterization 
of space. Both transformations are field dependent [24]. 

The time transformation is given by an implicit formula involving also the field 6, 



t->T(t,r) = — \ ^ T '^ tanhcj, (2.1.45) 

cosh oj c l 



while the field transforms according to 



6(t,r) -> 0Jt,r) = v ' ' + c 2 t tanhcu . (2.1.46) 
cosh a; 

[We record here only the transformation on 6; how p transforms can be determined from the 
(relativistic) Bernoulli equation, (2.1.42), which expresses p in terms of 9. Moreover, (2.1.46) 
is sufficient for discussing the invariance of (2.1.38).] The infinitesimal generator, which is time 
independent by virtue of the equations of motion, is [25] 

S = Jdr (cHp + 9^p 2 c 2 + oVc 2 + (V#) 2 ) , 
= J dr(c 4 tp + 9TC) (time reparameterization). (2.1.47) 



42 



R.Jackiw 



A second class of transformations involving a reparameterization of the spatial variables is 
implicitly defined by 



R{t, r) = r + ^9{t, R) + w(w ■ r) ( \j 1 + ^ - 1 ) , (2.1.48) 



a; 2 



0(t, r) -> ^(t, r) = y 1 + — 0(i, R) + w • r. (2.1.49) 
The time-independent generator of the infinitesimal transformation reads [25] 

G = J dr(c 2 rp + 9pV9), 

= J dr(c 2 vp + OV) (space reparameterization). (2.1.50) 

Of course the Born-Infeld action (2.1.36) or (2.1.38) is invariant against these transformations, 
whose infinitesimal form is generated by the constants. 

With the addition of S and G to the previous generators, the Poincare algebra in (d + 1, 1) 
dimension is reconstructed, and (t, r, 9) transforms linearly as a (d + 2)-dimensional Lorentz vector 
(in Cartesian components). Note that this symmetry also holds in the free, a = 0, theory. 

It is easy to exhibit solutions of the relativistic equation (2.1.43), which reduce to solutions of 
the nonrelativistic, Chaplygin gas equation (2.1.18) [after — c 2 t has been removed, as in (2.1.34)]. 
For example 



r 2 



9{t,r) = -c\lcH 2 + — - (2.1.51) 
solves (2.1.43) and reduces to (2.1.21). The relativistic analog of the lineal solution (2.1.23) is 

6(t, r) = e(n ■ r) + u • r - ct^/c 2 + u 2 - (n ■ u) 2 . (2.1.52) 

[Note that the above profiles continue to solve (2.1.43) even when the sign of the square root is 
reversed, but then they no longer possess a nonrelativistic limit.] 

Additionally there exists an essentially relativistic solution, describing massless propagation in 
one direction: according to (2.1.43), 9 can satisfy the wave equation 09 = 0, provided (d^9) 2 = con- 
stant, as for example with plane waves, 

9{t,r) = /(n-r±rf), (2.1.53) 

where {d^9) 2 vanishes. Then p reads, from (2.1.42), 

P = T- 2 f. (2.1.54) 
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2.2 Common ancestry: the Nambu-Goto action 

The "hidden" symmetries and the associated transformation laws for the Chaplygin and Born-Infeld 
models may be given a coherent setting by considering the Nambu-Goto action for a p-brane in 
(p + 1) spatial dimensions, moving on (p+1, l)-dimensional space-time. In our context, a p-brane 
is simply a p-dimensional extended object: a 1-brane is a string, a 2-brane is a membrane and so 
on. A p-brane in (p + 1) space divides that space in two. 
The Nambu-Goto action reads 

/ng = J d0°d<K NG = J d^d^-'-d^VG, (2.2.1) 
G= (-l) p det^^. (2.2.2) 



Here X M is a (p + 1, 1) "target space-time" (p-brane) variable, with p extending over the range 
p — 0,1, ... ,p,p + 1. The Q are "world-volume" variables describing the extended object with a 
ranging a = 0,1, ... ,p; <f) r , r = 1, . . . ,p, parameterizes the p-dimensional p-brane that evolves in 
0°. 

The Nambu-Goto action is parameterization invariant, and we shall show that two different 
choices of parameterization ("light-cone" and "Cartesian") lead to the Chaplygin gas and Born- 
Infeld actions, respectively. For both parameterizations we choose (X 1 , . . . ,X P ) to coincide with 
{<p l , . . . , <jf), renaming them as r (a p-dimensional vector). This is usually called the "static param- 
eterization" . (The ability to carry out this parameterization globally presupposes that the extended 
object is topologically trivial; in the contrary situation, singularities will appear, which are spuri- 
ous in the sense that they disappear in different parameterizations, and parameterization-invariant 
quantities are singularity-free.) 

(i) Light-cone parameterization 

For the light-cone parameterization we define X as -^{X° ± X p+1 ). X + is renamed t and 
identified with y/2X(f> . This completes the fixing of the parameterization and the remaining variable 
is X~ , which is a function of 0° and </>, or after redefinitions, of t and r. X~ is renamed as 6(t, r) 
and then the Nambu-Goto action in this parameterization coincides with the Chaplygin gas action 

7 ChaplygiB in (2 L17) [2g] _ 

(ii) Cartesian parameterization 

For the second, Cartesian parameterization X° is renamed ct and identified with c<j)°. The re- 
maining target space variable X p+l , a function of cjP and (p, equivalently of t and r, is renamed 
9{t,r)/c. Then the Nambu-Goto action reduces to the Born-Infeld action / dt L^ orn - Infcld , (2.1.38) 
[26]. 
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(iii) Hodographic transformation 



There is another derivation of the Chaplygin gas from the Nambu-Goto action that makes use 
of a hodographic transformation, in which independent and dependent variables are interchanged. 
Although the derivation is more involved than the light-cone/static parameterization used in Section 
2.2(i) above, the hodographic approach is instructive in that it gives a natural definition for the 
density p, which in the above static parameterization approach is determined from 9 by the Bernoulli 
equation (2.1.5). 

We again use light-cone combinations: -^(X° + X p+l ) is called r and is identified with <fr°, while 
-^(X° — X p+1 ) is renamed 9. At this stage the dependent, target-space variables are 9 and the 
transverse coordinates X: X 1 , (i = 1, . . . ,p), and all are functions of the world- volume parameters 
0° = r and (p: <p r , (r = l,...,p); d T indicates differentiation with respect to r = 0°, while d r 
denotes derivatives with respect to (jf . The induced metric G a p = ^^^njr takes the form 

G = (G 00 G os \f 2d T 9 - (d T Xf d s 9 - d T X ■ d s x\ 2 
aP \G ro -g rs ) \d r 9 - <9 r X • <9 r X -<9 r X ■ <9 S X J ' 1 ' ' ' 

The Nambu-Goto Lagrangian now leads to the canonical momenta 

9£ng 



p, (2.2.4) 
IT, (2.2.5) 



dd T X 

dC NG 
dd T 9 

and can be presented in first-order form as 

£ NG = p • d T X + Iid T 9 + ^ (p 2 + g) + u r (p ■ d r X + Ud r 9), (2.2.6) 
where g = det g rs and 

u r = d T X ■ d r X - d r 9 (2.2.7) 
acts as a Lagrange multiplier. Evidently the equations of motion are 

d T X = --p - u r d r X, (2.2.8a) 

d T 9 = ^(p 2 + g)-u r d r 9, (2.2.8b) 

d rP = -d r (^gg rs d s X) - d r (u r p), (2.2.8c) 
d T U = -d r (u r Tl). (2.2.8d) 



Also there is the constraint 



p • d r X + Ud r 9 = 0. 



(2.2.9) 
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[That u r is still given by (2.2.7) is a consequence of (2.2.8a) and (2.2.9).] Here g rs is inverse to g rs , and 
the two metrics are used to move the (r, s) indices. The theory still possesses an invariance against 
redefining the spatial parameters with a r-dependent function of the parameters; infinitesimally: 
5cf) r = — / r (r, </>), 59 = f r d r 8, 5X l = f r d r X l . This freedom may be used to set u r to zero and II 
to -1. 

Next the hodographic transformation is performed: Rather than viewing the dependent variables 
p, 9, and X as functions of r and </>, X(r, (f>) is inverted so that (f) becomes a function of r and X 
(renamed t and r, respectively), and p and 9 also become functions of t and r. It then follows from 
the chain rule that the constraint (2.2.9) (at II = —1) becomes 

= ^-4^ (2-2.10) 



dX 

and is solved by 

p = V0. (2.2.11) 

Moreover, according to the chain rule and the implicit function theorem, the partial derivative with 
respect to r at fixed <fi [this derivative is present in (2.2.6)] is related to the partial derivative with 
respect to r at fixed X = r by 

d T = d t + V9 -V, (2.2.12) 

where we have used the new name "t" on the right. Thus the Nambu-Goto Lagrangian - the 
integral of the Lagrange density (2.2.6) (at u r = 0, II = —1) - reads 

L NG = J d</>(p -V9-9-V9-V9- |(p 2 + (2.2.13a) 
But use of (2.2.11) and of the Jacobian relation d0 = dr det = shows that 

(W 

With the definition 

V^=V2A/p, (2.2.13c) 
Lng becomes, apart from a total time derivative 



L NG = Ur(--^-9-^—{V9f-\^g). (2.2.13b) 



7 NG 



^j ir (er- P ?f-^. (2.2.i3d) 



Up to an overall factor, this is just the Chaplygin gas Lagrangian in (2.1.4). 

The present derivation has the advantage of relating the density p to the Jacobian of the X — ► 
transformation: p = \/2A det J^-. (This in turn shows that the hodographic transformation is just 
exactly the passage from Lagrangian to Eulerian fluid variables.) 
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2.3 Interrelations 

The relation to the Nambu-Goto action explains the origin of the hidden (p+1, 1) Poincare group in 
our two nonlinear models on (p, 1) space-time: the (p+ 1, l)Poincare invariance is what remains of 
the reparameterization invariance of the Nambu-Goto action after choosing either the light-cone or 
Cartesian parameterizations. (In this context, recall that the light-cone subgroup of the Poincare 
group is isomorphic to the Galileo group in one lower dimension [27]) Also the nonlinear, field 
dependent form of the transformation laws leading to the additional symmetries is understood: it 
arises from the identification of some of the dependent variables (A^ M ) with the independent variables 

(<n 

The complete integrability of the d = 1 Chaplygin gas and Born-Infeld model is a consequence 
of the fact that both descend from a string in 2-space; the associated Nambu-Goto theory being 
completely integrable. We shall discuss this in Section 3. 

We observe that in addition to the nonrelativistic descent from the Born-Infeld theory to the 
Chaplygin gas, there exists a mapping of one system on another, and between solutions of one 
system and the other, because both have the same p-brane ancestor. The mapping is achieved by 
passing from the light-cone parameterization to the Cartesian, or vice-versa. Specifically this is 
accomplished as follows. 

(i) Chaplygin gas — > Born-Infeld: 

Given 9jsf R (t,r), a nonrelativistic solution, determine T(t,r) from the equation 

T+^9 NR (T,r) = V2t. (2.3.1) 

Then the relativistic solution is 



c 2 



9 R (t, r) = i=c 2 T - A=0 NR (T, r) = c\^2T - t). (2.3.2) 



(ii) Born-Infeld — ► Chaplygin gas: 

Given 9 R (t,r), a relativistic solution, find T(t,r) from 

T + \e R (T,r) = V2t. (2.3.3) 

cr 

Then the nonrelativistic solution is 

9 NR (t, r) = i=c 2 T - A=0 R {T, r) = c 2 (V2T - t). (2.3.4) 

The relation between the different models is depicted in the Figure below. 
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Dualities and other relations between nonlinear equations. 



A final comment: Recall that the elimination of p, both in the nonrelativistic (Chaplygin) 
and relativistic (Born-Infeld) models is possible only in the presence of interactions. Nevertheless, 
the ^-dependent (p-independent) resultant Lagrangians contain the interaction strengths only as 
overall factors; see (2.1.17) and (2.1.38). It is these ^-valued Lagrangians that correspond to the 
Nambu-Goto action in various parameterizations. Let us further recall the the Nambu-Goto action 
also carries an overall multiplicative factor: the p-brane "tension", which has been suppressed 
in (2.2.1). Correspondingly, for a "tensionless" p-brane, the Nambu-Goto expression vanishes, 
and cannot generate dynamics. This suggests that an action for "tensionless" p-branes could be 
the noninteracting fluid mechanical expressions (2.1.4), (2.1.36), with vanishing coupling strengths 
A and a, respectively. Furthermore, we recall that the noninteracting models retain the higher, 
dynamical symmetries, appropriate to a p-brane in one higher dimension. 
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3 SPECIFIC MODELS (d = 1) 



In this Section, we shall discuss nonrelativistic/relativistic models in one spatial dimension. Com- 
plete integrability has been established for both the Chaplygin gas [28] and the Born-Infeld theory 
[29] . We can now understand this to be a consequence of the complete integrability of the Nambu- 
Goto 1-brane (string) moving on 2-space (plane), which is the antecedent of both models. [Therefore, 
it suffices to discuss only the Chaplygin gas, since solutions of the Born-Infeld model can then be 
obtained by the mapping (2.3.1)-(2.3.2).] 

As remarked previously, in one dimension there is no vorticity, and the nonrelativistic velocity v 
can be presented as a derivative with respect to the single spatial variable of a potential 9. Simi- 
larly, the relativistic momentum p = v/ y/l — v 2 /c 2 is a derivative of a potential 9. In both cases 
the potential is canonically conjugate to the density p governed by the canonical 1-form J dx9p. 
Moreover, it is evident that at the expense of a spatial nonlocality, one may replace 9 by its an- 
tiderivative, which is p both nonrelativistically and relativistically (nonrelativistically p = v), so 
that in both cases the Lagrangian reads 





For the Chaplygin gas and the Born-Infeld models, H is given respectively by 




(3.0.2) 



(3.0.3) 



The equations of motion are, respectively 



Chaplygin gas: 




(3.0.4) 



(3.0.5) 







1 







P 



(3.0.6) 



or 




+ 







Born-Infeld model: 



. 8_ f / p 2 c 2 + a 2 \ = Q 
^ dx \ y c 2 + p 2 ) 



(3.0.7) 
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In the above, eqs. (3.0.6) and (3.0.9) result by determining p in terms of 9 (p = 9'; dash indicates 
different at ion with respect to spatial argument) from (3.0.5) and (3.0.8), and using that expression 
for p in (3.0.4) and (3.0.7). 

3.1 Chaplygin gas on a line 

(i) Specific solutions 

Classes of solutions for a Chaplygin gas in one dimension can be given in closed form. For 
example, to obtain general, time-rescaling-invariant solutions, we make the Ansatz that 9 oc 1/t. 
Then (2.1.18) or (3.0.6) leads to a second-order nonlinear differential equation for the x-dependence 
of 9. Therefore solutions involve two arbitrary constants, one of which fixes the origin of x (we 
suppress it); the other we call k, and take it to be real. The solutions then read 

0(t, x) = -— 2- cosh 2 kx . (3.1.1) 
2/u t 

[Other solutions can be obtained by relaxing the reality condition on k and/or shifting the argument 
kx by a complex number. In this way one finds that 9 can also be ^27 sinh 2 kx, ^27 sin 2 kx, 
cos 2 kx; but these lead to singular or unphysical forms for p.] The density corresponding to 
(3.1.1) is found from (2.1.5) or (3.0.5) to be 

k \t 



p(t,x) = V2A . (3.1.2) 



cosh kx 



The velocity/momentum v = p = 9' is 



v(t,x) = p(t,x) = — — sinh kx cosh kx, (3.1.3) 

Kb 

while the sound speed 

cosh 2 kx 

s(t,x) = — — p (3.1.4) 

is always larger than \v\. Finally, the current j = p|| exhibits a kink profile, 

j(t,x) = -e(t)V2Xtanhkx, (3.1.5) 

which is suggestive of complete integrability. 

Another particular solution is the Galileo boost of the static profiles (2.1.24), (2.1.25): 



p(t, x) = p(x — ut), 



(3.1.6) 
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Here u is the boosting velocity and p(x — ut) is an arbitrary function of its argument (provided 
p 7^ u). Clearly this is a constant profile solution, in linear motion with velocity u [30]. 

Further evidence for complete integrability is found by identifying an infinite number of constants 
of motion. One verifies that the following quantities 



I 



n 



dxp[p±^— J , n = 0,±l, ... (3.1.8) 



are conserved. 



The combinations p ± are just the velocity (±) the sound speed, and they are known as 



/2A 

V =c - 

Riemann coordinates. 

\/2\ 

R±=p±- — (3.1.9) 
P 

The equations of motion for this system [continuity (3.0.4) and Euler (3.0.5)] can be succinctly 
presented in terms of R± : 

R± = —R T R'± . (3.1.10) 



(ii) General solution for the Chaplygin gas on a line 

The general solution to the Chaplygin gas can be found by linearizing the governing equations 
(continuity and Euler) with the help of a Legendre transform, which also effects a hodographic 
transformation that exchanges the independent variables (t,x) with the dependent ones (p,9); 
actually instead of p we use the sound speed s = V2X/p and instead of 6 we use the momentum 
p = 0'. 

Define 

V>(p,s) = 9{t,x)-t9{t,x)-x9'(t,x) . (3.1.11) 
From the Bernoulli equation we know that 

6 = -\p 2 + \s 2 . (3.1.12) 

Thus 

^(p, s) = 9(t, x) + \{P 2 ~ s 2 ) - xp, (3.1.13) 
and the usual Legendre transform rules govern the derivatives. 

dip 



tp-x (3.1.14a) 
op 

= -ts (3.1.14b) 

OS 
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It remains to incorporate the continuity equation (3.0.4) whose content must be recast by the 
hodographic transformation. This is achieved by rewriting equation (3.0.4) in terms of s = v2A/ p 

dt ^ dx dx ^ 
Next (3.1.15) is presented as a relation between Jacobians 

which is true because here dx/dt = dt/dx = 0. Eq. (3.1.16a) implies, after multiplication by 

d(t,x)/d(s,p) 

Q _ d(M) ^ <9(t,p) 

d(s,p) d(s,p) 9(s,p) 

= T-*T-ir- < 3 - hl6b > 

ap as 

The second equality holds because now we take ds/dp = dp/ds = 0. Finally, from (3.1.13), (3.1.14) 
it follows that (3.1.16b) is equivalent to 

dp 2 ds 2 s ds 

This linear equation is solved by two arbitrary functions of p ± s (p ± s being just the Riemann 
coordinates) 

V>(p, s) = F(p + s) - sF'(p + s) + G(p -s) + sG'{p - s) (3.1.17) 

In summary, to solve the Chaplygin gas equations, we choose two functions F and G, construct 
ifj as in (3.1.17), and regain s (= y/2X/p), p (= 6'), and 6 from (3.1.13), (3.1.14). In particular, the 
solution (3.1.1), (3.1.2) corresponds to 

F(z) =G{-z) = ±^-\nz, (3.1.18) 

where the sign is correlated with the sign of t. 



C. Sidebar on the integrability of the cubic potential in d=l 

Although it does not belong to the models that we have discussed, the cubic potential for 1- 
dimensional motion, V(p) = £p 3 /S, is especially interesting because it is secretly free [see (1.2.34b)] 
- a fact that is exposed when Riemann coordinates are employed. For this problem these read 
R± = p ± \r2£p and again they are just the velocity (±) the sound speed. In contrast to (3.1.10) 
the Euler and continuity equations for this system decouple: R± = —R±R'±. Indeed, it is seen that 
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R± satisfy essentially the free Euler equation [compare with (1.2.61b) at V" = and identify R± 
with v]. Consequently, the solution (1.2.62)-(1.2.66) works here as well. 

Recall the previous remark in Section 1.2 (iv) on the Schrodinger group [Galileo and SO(2,l)]: 
in one dimension the cubic potential is invariant against this group of transformations, and in 
all dimensions the free theory is invariant. Therefore a natural speculation is that the secretly 
noninteracting nature of the cubic potential in one dimension is a consequence of Schrodinger group 
invariance. 

Another interesting fact about a one- dimensional nonrelativistic fluid with cubic potential is 
that it also arises in a collective, semiclassical description of nonrelativistic free fermions in one 
dimension, where the cubic potential reproduces fermion repulsion [31]. In spite of the nonlinearity 
of the fluid model's equations of motion, there is no interaction in the underlying fermion dynamics. 
Thus, the presence of the Schrodinger group and the equivalence to free equations for this fluid 
system is an understandable consequence. 



3.2 Born-Infeld model on a line 

Since the Born-Infeld system is related to the Chaplygin gas by the transformation described in 
Section 2.3, there is no need to discuss separately Born-Infeld solutions. Nevertheless, the for- 
mulation in terms of Riemann coordinates is especially succinct and gives another view on the 
Chaplygin/Born-Infeld relation. 

The Riemann coordinates R± for the Born-Infeld model are contructed by first defining 

— 6' = p/c = tan0 p , 

a/pc = tan0 p , (3.2.1) 

and 

R± = <P P ± P . (3.2.2) 

The 1-dimensional version of the equations of motion (2.1.41), (2.1.42), that is, (3.0.7), (3.0.8) can 
be presented as 

R ± = - c (smR T )R' ± . (3.2.3) 

The relation to the Riemann description of the Chaplygin gas can now be seen in two ways: a 
nonrelativistic limit and an exact transformation. For the former, we note that at large c, <f) p ~ p/c, 
4>p ~ a/pc so that 

1 / n\ 1 .„, . . 

(3.2.4) 



^.Born-Infeld ^ = - ^'Chaplygin 

± C\ pJ C ± 



X=a 2 /2 
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Moreover, the equation (3.2.3) becomes, in view of (3.2.4), 

-i? Chaplygin = j^Chaplygin \ & jgChaplygin (3 2 5) 

c ± T c dx ± 

so that (3.1.10) is regained. On the other hand, for the exact transformation we define new Riemann 
coordinates in the relativistic, Born-Infeld case by 

K± = csinR± . (3.2.6) 

Evidently (3.2.3) implies that TZ± satisfies the nonrelativistic equations (3.1.10), (3.2.5) when R± 
solves the relativistic equation (3.2.3). Expressing TZ± and R± in terms of the corresponding nonrel- 
ativistic and relativistic variables produces a mapping between the two sets. Calling p NR , p NR and 
Pr, Pr the momentum and density of the nonrelativistic and of the relativistic theory, respectively, 
the mapping implied by (3.2.6) is 

c 2 PrPr 

Pnr - 



V(p'r + c 2 )(pW + * 2 ) 

PNR = ^ V /(p R + c 2 )(p| C 2 + a2) . (3.2.7) 

As can be checked, this maps the Chaplygin equations into the Born-Infeld equations. But the 
mapping is not canonical. 

We record the infinite number of constants of motion, which put into evidence the (by now 
obvious) complete integrability of the Born-Infeld equations on a line. The following quantities are 
time-independent: 

I± = ac n - 1 [ dx {<f)p± , n = 0,±l,... (3.2.8) 
J sin <fip cos (f>p 

The nonrelativistic limit takes the above into (3.1.8), while expressing 1^ in terms of 1Z± according 
to (3.2.6) shows that the integrals in (3.2.8) gives rise to a series of the integrals in (3.1.8). 

In the relativistic model p need not be constrained to be positive (negative p could be interpreted 
as antiparticle density). The transformation p — > —p, p — > —p is a symmetry and can be interpreted 
as charge conjuguation. 

Further, p and p appear in an equivalent way. As a result, this theory enjoys a duality transfor- 
mation. 

P^±4p P^±-P (3.2.9) 

c z a 

Under the above, both the canonical structure and the Hamiltonian remain invariant. Solutions 
are mapped in general to new solutions. Note that the nonrelativistic limit is mapped to the 
ultra-relativistic one under the above duality. Self-dual solutions, with p = ± Jrp, satisfy 

P = Tc-^p', (3.2.10) 



54 



R.Jackiw 



and are, therefore, the chiral relativistic solutions that were presented at the end of Section 2.1(h). 
In the self-dual case, when p is eliminated from the canonical 1-form and from the Hamiltonian with 
the help of (3.2.9), one arrives at an action for p, which coincides (apart from irrelevant constants) 
with the self-dual action, constructed some time ago [32]. 

\ J d£dxdy/j(x)e(x — y) P (y) - f it d* ^TO^T^d* 
i j AtAxAyp(x)t{x - y)p{y) -^J dtdx(p 2 (x) + ^) } (3.2.11) 



c 2 



2c^ 
a 



3.3 General solution of the Nambu-Goto theory for a (p=l)-brane 
(string) in two spatial dimensions (on a plane) 

The complete integrability of the 1-dimensional Chaplygin gas and Born-Infeld theory, as well as 
the relationships between the two, derives from the fact that the different models descend by fixing 
in different ways the parameterization invariance of the Nambu-Goto theory a for string on a plane. 
At the same time, the equations governing the planar motion of a string can be solved completely. 
Therefore it is instructive to see how the string solution produces the Chaplygin solution [20]. 
We follow the development in Section 2.2(iii). The Nambu-Goto action reads 



(3.3.1a) 
(3.3.1b) 
(3.3.1c) 



Here X^, fi = 0, 1,2, are string variables and ((p ,^ 1 ) are its parameters. As in Section 2.2(iii), we 
define light-cone combinations X ± = -^(X° ± X 2 ), rename X~ as 9, and choose the parameteriza- 
tion X + = (f>° = r. After suppressing the superscripts on (p 1 and X 1 , we construct the Nambu-Goto 
Lagrange density as 





f d0°L NG , 


Lng = j 


[ d^ 1 £ NG , 


£ng — 


— det — — -5- 



, ,„ 1 2d T e - (d T xf u 



u = dj,9 — d T Xdd,X. 



(3.3.2) 
(3.3.3) 



Equations of motion are presented in Hamiltonian form: 

<9£ng tt <9£ng 



V 



n 



dd T 9 ' 



(3.3.4) 
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and there is the constraint 



d T X = ~Y[P — ud^X, (3.3.5a) 

d T 9 = (p 2 + {d^Xf) - ud^O, (3.3.5b) 

d T p = -d^—d^X^j - d^up), (3.3.5c) 

d T U = -d^uU), (3.3.5d) 



pd^X + Ud^O = . (3.3.6) 



There still remains the reparameterization freedom of replacing by an arbitrary function of r 
and 0; this freedom may be used to set u = 0, II = — 1. Consequently, in the fully parameterized 
equations of motion eq. (3.3.5d) disappears; instead of (3.3.5a) and (3.3.5c), we have d T X = p, 
d T p = d^X, which imply 

(d T 2 - %)X = 0, (3.3.7a) 

(3.3.5b) reduces to 

<9 T = ±[(<9 T X) 2 + OVO 2 L (3.3.7b) 

and the constraint (3.3.6) requires 

0^6 = drXd^X . (3.3.7c) 
Solution to (3.3.7a) is immediate in terms of two functions F±, 

x(r,0) = F+(r + 0) + F_(r-0) (3.3.8) 

and then (3.3.7b), (3.3.7c) fix 9. 

/r-\-(j) rr—(p 

dz[F^z)] 2 + j dz[FL(z)] 2 (3.3.9) 



This completes the description of a string moving on a plane. 

But we need to convert this information into a solution of the Chaplygin gas, and we know from 
Section 2.2(iii) that this can be accomplished by a hodographic transformation: instead of X and 
6 as a function of r and 0, we seek as a function of r and X, and this renders 9 to be a function 
of r and X as well. The density p is determined by the Jacobian \dX/d(f>\. 

Replace r by t and X by x and define to be f(t,x). Then from (3.3.8) it follows that 



x = F + (t + f{t,x))+F-(t-f{t,x)) 



(3.3.10) 
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This equation may be differentiated with respect to t and x, whereupon one finds 

FUt + f)+FL(t-f) 



f 

r 



F+(t + f) — F-(t — f)' 
1 



(3.3.11a) 
(3.3.11b) 



F' + {t + f) - F'_{t - f) 

Thus the procedure for constructing a Chaplygin gas solution is to choose two functions F± 
solve the differential equations (3.3.11) for /, and then the fluid variables are 



9(t,x) 
V2X 



t+f(t,x) 



p 



[F' + (z)Ydz + 
F>(t + 4>)-FL{t-<]>, 



t-f{t,x) 



[FL{z)] 2 dz, 



(3.3.12) 
(3.3.13) 



One may verify directly that (3.3.12) and (3.3.13) satisfy the required equations: Upon differ- 
entiating (3.3.12) with respect to t and x, we find 

^ = (f;) 2 (i + /) + (F) 2 (i-/) 

= -2F' + FL (3.3.14a) 

o' = (K) V - (FLff 

= F' + + F'_. 

The second equalities follow with the help of (3.3.11). From (3.3.14) one sees that 

A 



(3.3.14b) 



?> \2 



P 1 



(3.3.15) 



the last equality being the definition (3.3.13). Thus the Bernoulli (Euler) equation holds. For the 
continuity equation, we first find from (3.3.13) and (3.3.14) 

. d V2X 



P 



' dtF' + — F'_ 



2\ 



=F 



{F' + - FL) 2 
2V2X 



F"(l + f) - F'l(l - f) 



(n 



(F»FL + F»F[).. 



(3.3.16a) 



d_ 

dx 



f 2\ 

(n - F ' 



T fj?> )?i \2 (F+FL + F'!_F' + ) f 



(3.3.16b) 
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The last equalities follow from (3.3.11); since (3.3.16a) and (3.3.16b) sum to zero, the continuity 
equation holds. 

We observe that the differentiated functions F± are just the Riemann coordinates: from (3.3.14b) 
and (3.3.13) [with the absolute value ignored] we have 



v2A 

p ± —- = R ± = 2F' ± . (3.3.17) 

Also it is seen with the help of (3.3.11) that the Riemann formulation (3.1.10) of the Chaplygin 
equations is satisfied by 2F±. 

The constants of motion (3.1.8) become proportional to 

=/d*g[Fi (t± /)r 

oc J &z[F' ± (z)] n . (3.3.18) 
Finally we remark that the solution (3.1.1), (3.1.2) corresponds to 

F + (z) = -F.(z) = ±^ . (3.3.19) 

There exists a relation between the two functions F and G in (3.1.17), which encode the Chaplygin 
gas solution in the linearization approach of Section 3.1(h), and the above two functions F±, which 
do the same job in the Nambu-Goto approach. The relation is that 2F\ is inverse to 2F" and 2F'_ 



is inverse to 2G" , that is, 



2F"[2F[(z)} = z, 

2G"[2F'_(z)\ = z. (3.3.20) 
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4 SUPERSYMMETRIC FLUID MECHANICS 

As explained in Section 1, classical fluids describe particles moving collectively. The fluid inherits its 
mechanical properties, such as energy, nomentum and angular momentum from the corresponding 
underlying particle properties. 

One consequence of this is that classical fluids cannot carry intrinsic spin. To be precise, the 
angular momentum with respect to the center of mass of a small volume V scales as its mass (which 
scales as V) times the residual velocity of the fluid about the center of mass (which scales like I- the 
linear dimension of V) times the distance to the center of mass (which also scales like £). Therefore, 
the self- angular momentum density scales lie £ 2 , and goes to zero with I. 

The inclusion of a spin density in fluids can be achieved in an essentially quantum mechani- 
cal formulation by introducing Grassmann (anticommuting) variables in the description, the spin 
density being represented as a bilinear in the Grassmann variables. This description reveals the pos- 
sibility of implementing within fluid mechanics supersymmetry transformations, which effectively 
mix spin and kinematical degrees of freedom. Particular forms of the Hamiltonian, generalizing the 
classical Chaplygin gas, admit these supersymmetry transformations as an invariance and generate 
conserved quantities. 

Supersymmetry poses severe restrictions in the dynamics. In general, supersymmetric extended 
objects cannot be formulated in arbitrary dimensions, and this holds true for supersymmetric fluids. 
It is natural, therefore, that the supersymmetric Chaplygin models are essentially related to, and 
derive from, higher-dimensional supersymmetric membrane actions in a way similar to the one 
already exposed in Section 2.2. As such, they also enjoy nontrivial higher-dimensional relativistic 
symmetries which are not apparent from their action. 

In the following we shall analyze the case of planar [33] and lineal [34] fluids, which devolve from 
the motion of membranes or strings in (3 + 1) or (2 + 1) dimensional spacetimes, respectively. 

Lineal and planar supersymmetric fluid models seem to exhaust the possibilities for the super- 
symmetric Nambu-Goto/fluid connection. For a higher dimensional generalization, the reduction 
program would begin with a p-brane in p + 2 dimensional space-time, giving rise to a fluid in p + 1 
dimensional space-time. While there are no constraints on p in the purely bosonic case, supersym- 
metric extensions are greatly constrained: the list of possible "fundamental" super p-branes (i.e. 
with only scalar supermultiplets in the world-volume) contains only the above two cases: p = 2 in 
four dimensions and p = 1 in three dimensions [35] . 

4.1 Supersymmetric fluid in (2 + 1) dimensions 

We begin by positing the fluid model. The Chaplygin gas Lagrangian is supplemented by Grassmann 
variables ip a that are Majorana spinors [real, two-component: if)* = i/) a , a = 1,2, (^1^2)* = tpl^]- 
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The associated Lagrange density reads 

A \/2A 

£ = - p (9 - 1^0) - \p(V9 - i^V^) 2 ^— • ■ (4.1.1) 

P 2 

Here a 1 are two (« = 1,2), 2x2, real symmetric Dirac "alpha" matrices; in terms of Pauli matrices 
we can take a 1 = a 1 , a 2 = a 3 . Note that the matrices satisfy the following relations, which are 
needed to verify subsequent formulas. 

e ab a l bc = e lJ a J ac 

<b4c = $ ij Sac - e i3 e ac 

a lb a ld - Sacfad - SabScd + S ad 5 bc (4.1.2) 

e a b is the 2x2 antisymmetric matrix e = ia 1 . In equation (4.1.1) A is a coupling strength which is 
assumed to be positive. The Grassmann term enters with coupling -\/2A, which is correlated with 
the strength of the Chaplygin potential V(p) = A/p in order to ensure supersymmetry, as we shall 
show below. It is evident that the velocity should be defined as 

v = V6-\^Vip. (4.1.3) 

The Grassmann variables directly give rise to a Clebsch formula for v, and provide the Gauss 
potentials. The two-dimensional vorticity reads uj = e^d^ = — \e 1 ^ d^djip = — x V?/>. The 
variables {0,p} remain a canonical pair, while the canonical 1-form in (4.1.1) indicates that the 
canonically independent Grassmann variables are y/pip so that the antibracket of the ^'s is 

{^r),VUr')} = -4V(r-r') ■ ( 4L4 ) 
p(r) 

One verifies that the algebra (1.2.16) - (1.2.21) is satisfied, and further, one has 

{0(r), ^(r)} = -_L-V(r)*(r - r'), (4.1.5) 
2p(r) 

{v(r),^(r0} = -^^5(r-r'), (4.1.6) 
p(r) 

{V(r),ip(r')} = -V^(r)£(r - r') . (4.1.7) 

The momentum density "P is given by the bosonic formula *P = vp, but the Grassmann variables 
are hidden in v, by virtue of (4.1.3). 
The equations of motion read 

p + V-(vp) = 0, (4.1.8) 

+ v V9 = I v 2 + 4 + ^^«- W, (4.1.9) 
z p z 2p 

ip + v-V?p = ex-V^, (4.1.10) 

P 
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and together with (4.1.1) they imply 



v + v ■ Vv = v4r 



2A 



Vtp)cx ■ Vip . 



P P 

All these equations may be obtained by bracketing with the Hamiltonian 

d 2 rH 



H = J d 2 r(lpv 2 + - + " V^) = J 



(4.1.11) 



(4.1.12) 



P 2 

when (1.2.16), (1.2.20) as well as (4.1.4)-(4.1.6) are used. 

We record the components of the energy-momentum "tensor" , and the continuity equations they 
satisfy. The energy density S = T°°, given by 

A V2A 
P 



£ = > 2 



ipa ■ = T°°, 



(4.1.13) 



satisfies a continuity equation with the energy flux T^°. 



A 



T^ = p^'(iv 2 --) 



rpO 



— ipa j v ■ — —ipdjip + —e* k ibedkib 
2 p p 

djT jo = 



(4.1.14) 
(4.1.15) 



This ensures that the total energy, that is, the Hamiltonian, is time- independent. Conservation of 
the total momentum 

(4.1.16) 



d 2 rV 



follows from the continuity equation satisfied by the momentum density V 1 = T m and the momen- 
tum flux, that is, the stress tensor T % K 



rpjl 



pv l v j — 5 lj (— 



2A 



P 2 
f oi + diT j 



ipcx ■ Vt/> ) + 



ipa 3 diip 







(4.1.17) 
(4.1.18) 



But T u is not symmetric in its spatial indices, owing to the presence of spin in the problem. 
However, rotational symmetry makes it possible to effect an "improvement", which modifies the 
momentum density by a total derivative term, leaving the integrated total momentum unchanged 
(provided surface terms can be ignored) and rendering the stress tensor symmetric. The improved 
quantities are 



T\° = pv l v J 



Vj = Tf = pv l + ^d 3 (p^) 
_ 6ij( 2X V2X 



- | d k \(e kl v j + e^v^fnpeip 
Tf + djT? = 



ifja l djifj + ipa^ditf) 



(4.1.19) 

(4.1.20) 
(4.1.21) 
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It immediately follows from the symmetry of Tf that the angular momentum 

M = J d 2 rj j r i V 3 j = J d 2 r pe ij r V + \ J d 2 r pipeip (4.1.22) 

is conserved. The first term is clearly the orbital part (which still receives a Grassmann contribution 
through v), whereas the second, coming from the improvement, is the spin part. Indeed, since 
|e = |cr 2 = E, we recognize this as the spin matrix in (2+1) dimensions. The extra term in 
the improved momentum density (4.1.19), |e u <9j (/??/> e-?/>), can then be readily interpreted as an 
additional localized momentum density, generated by the nonhomogeneity of the spin density. This 
is analogous to the magnetostatics formula giving the localized current density j m in a magnet in 
terms of its magnetization m: j m = V x m. All in all, we are describing a fluid with spin. 
Also the total number 

N = d 2 rp (4.1.23) 



is conserved by virtue of the continuity equation (4.1.8) satisfied by p. Finally, the theory is Galileo 
invariant, as is seen from the conservation of the Galileo boost, 

B = tP - J d 2 rvp (4.1.24) 

which follows from (4.1.8) and (4.1.16). The generators H, P,M, B and iV close on the (extended) 
Galileo group. [The theory is not Lorentz invariant in (2 + l)-dimensional space-time, hence the 
energy flux does not coincide with the momentum density, improved or not.] 

We observe that p can be eliminated from (4.1.1) so that £ involves only 9 and i/j. From (4.1.9) 
and (4.1.10) it follows that 

p = V\(6- \^ + \v 2 y 1/2 . (4.1.25) 

Substituting into (4.1.1) produces the supersymmetric generalization of the Chaplygin gas Lagrange 
density in (2.1.17). 

£ = -2V\ | yj2d-i/)ijj + (V0- lipVip) 2 + i^cx • V^} (4.1.26) 

Note that the coupling strength has disappeared from the dynamical equations, remaining only as a 
normalization factor for the Lagrangian. Consequently the above elimination of p cannot be carried 
out in the free case, A = 0. 



4.2 Supersymmetry 



As we said earlier, this theory possesses supersymmetry. This can be established, first of all, by 
verifying that the following two-component supercharges are time-independent Grassmann quanti- 
ties. 



Qa 



d 2 r 



pv ■ (a ab ip b ) + V2\ip a 



(4.2.1) 



(32 



R.Jackiw 



Taking a time derivative and using the evolution equations (4.1.8)-(4.1.11) establishes that Q a = 0. 

Next, the supersymmetric transformation rule for the dynamical variables is found by con- 
structing a bosonic symmetry generator Q, obtained by contracting the Grassmann charge with a 
constant Grassmann parameter rj a , Q = r] a Q a , and commuting with the dynamical variables. Using 
the canonical brackets one verifies the following field transformation rules. 

6p = {Q,p} = -V-p(riaif>) (4.2.2) 



\/2A 

59 = {Q, 9} = -Uv<*i>) ■ Vfl - Hrjaip) ■ ipViP + ^— ^ (4.2.3) 

zp 

v^2A 

Sijj = {Q, ip\ = -(rictip) ■ Vip — v • Oct/ T) (4.2.4) 

P 



v 2A 

5v = {Q, v} = —(rjctip) ■ Vv H r]V^p (4.2.5) 

P 

Supersymmetry is reestablished by determining the variation of the action J dt d 2 r £ consequent to 
the above field variations: the action is invariant. One then reconstructs the supercharges (4.2.1) 
by Noether's theorem. Finally, upon computing the bracket of two supercharges, one finds 

KQa,%U} = 2M 2 )#, (4.2.6) 

which again confirms that the charges are time-independent. 

{H,Q a } = (4.2.7) 

Additionally a further, kinematical, supersymmetry can be identified. According to the equa- 
tions of motion the following two supercharges are also time- independent. 

Q a = j d 2 rp^ a (4.2.8) 

Q = fj a Q a effects a shift of the Grassmann field. 

Sp = {Q,p} = (4.2.9) 

60 = {Q,6} = -±(ifil>) (4.2.10) 

fy = {Q,il>} = -fj (4.2.11) 

Sv = {Q,v} = (4.2.12) 

This transformation leaves the Lagrangian invariant, and Noether's theorem reproduces (4.2.8). 
The algebra of these charges closes on the total number TV, 

{v a 1 Q a ,r ] b 2 Q b } = (r ]l T ] 2)N, (4.2.13) 
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while the algebra with the generators (4.2.1), closes on the total momentum, together with a central 
extension, proportional to volume of space = J d 2 r. 

{v a Q a , V b Qb} = {fjotrj) ■ P + V2X (fjeri)n (4.2.14) 

The supercharges Q a ,Qa, together with the Galileo generators (H, P,M, B), and with N form a 
sup erext ended Galileo algebra. The additional, nonvanishing brackets are 

{M,Q a } = le ab Q b , (4.2.15) 
{M,Q a } = le ab Q b , (4.2.16) 
{B, Qa} — CXabQb ■ (4.2.17) 



4.3 Supermembrane connection 

The equations for the supersymmetric Chaplygin fluid devolve from a supermembrane Lagrangian, 
Lm- We shall give two different derivations of this result, which make use of two different parame- 
terizations for the parameterization-invariant membrane action and give rise, respectively, to (4.1.1) 
and (4.1.26). The two derivations follow what has been done in the bosonic case in Sections 2.2 (i) 
and 2.2 (hi). 

We work in a light-cone gauge-fixed theory: The supermembrane in 4-dimensional space-time 
is described by coordinates X M (/i = 0, 1,2,3), which are decomposed into light-cone components 
X ± = -^(X° ± X 3 ) and transverse components X % {i = 1,2}. These depend on an evolution 
parameter 0° = r and two space-like parameters <p r {r = 1,2}. Additionally there are two- 
component, real Grassmann spinors ip, which also depend on r and cjf . In the light-cone gauge, X + 
is identified with r, X~ is renamed 9, and the supermembrane Lagrangian is [36] 



J d 2 0£ A/ = - J d 2 <j){VG - \e rs d r ^cxd s ^ -X}, (4.3.1) 



where G = det G a 




G Q ,= \— ^o S > = (2d T 6-(d T X) 2 -^ 

'ro 9rs J y U r 

G = gT 

T = 2d T 6 - (<9 r X) 2 - ifjd T ip + g rs u r u s 
g rs = d r X ■ d s X , g = det g rs 

u s = d s 6 - ^d s tfj - <9 r X • d s X (4.3.3) 
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Here d T signifies differentiation with respect to the evolution parameter r, while d r differentiates 
with respect to the space-like parameters r ; g rs is the inverse of g rs , and the two are used to move 
the (r, s) indices. Note that the dimensionality of the transverse coordinates X 1 is the same as of 
the parameters r , namely two. 

(i) Hodographic transformation 

To give our first derivation following the procedure in Section 2.2(iii), we rewrite the Lagrangian in 
canonical, first-order form, with the help of bosonic canonical momenta defined by 

0C M 



dd T X 
dC M 



-Tld T X - Uu r d r X, (4.3.4a) 



99t6 n =v ^r. (4.3.4b) 

(The Grassmann variables already enter with first-order derivatives.) The supersymmetric extension 
of (2.2.6) then takes the form 



C M = P • d T X + Ud T 9 - \li^d T ^ + ^(p 2 + g) + \e rs d r ^otd s iP ■ X 



+ u r (p ■ <9 r X + Ud r 6 - lU^d r ip) . (4.3.5) 

In (4.3.5) u r serves as a Lagrange multiplier enforcing a subsidiary condition on the canonical vari- 
ables, and g = detg rs . The equations that follow from (4.3.5) coincide with the Euler-Lagrange 
equations for (4.3.1). The theory still possesses an invariance against redefining the spatial param- 
eters with a r-dependent function of the parameters. This freedom may be used to set u T to zero 
and fix II at —1. Next we introduce the hodographic transformation, as in Section 2.2(iii), whereby 
independent- dependent variables are interchanged, namely we view the (ff to be functions of X 1 . It 
then follows that the constraint on (4.3.5), which with II = —1 reads 

p • d r X - d r 9 + lipd r ip = 0, (4.3.6) 

becomes 

<9 r X • (p - V6 + l<ipVip) = . (4.3.7) 

Here p, 9 and i/j are viewed as functions of X, renamed r, with respect to which acts the gradient 
V. Also we rename p as v, which according to (4.3.7) is 

v = Vfl - ItpVip . (4.3.8) 

As in Section 2.2(iii), from the chain rule and the implicit function theorem it follows that 



d T = d t + d T X ■ V, 



(4.3.9) 
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and according to (4.3.4a) (at II = — 1, u r = 0) <9 r X = p = v. Finally, the measure transforms 
according to d 2 — > d 2 r -^=. Thus the Lagrangian for (4.3.5) becomes, after setting u T to zero and 
n to -1, 

/A 2 
-j= (v 2 - 6 - v V0 + \^ + v ■ V^) - |(v 2 + (?) 

- \e rs i\)ot djip d s X j d r X l ^j . (4.3.10a) 
But e rs d s X j d r X l = e ij det drX* = e ij ^/g. After y/g is renamed V2X/p, (4.3.10a) finally reads 

L M = -)= [ d 2 r (-p{9 - \^) - \p{V6 - ^V^) 2 - - - ^^>a x Vxp) . (4.3.10b) 
v2A J ^ p 2, J 



Upon replacing ip by ^(1 — e)^, this is seen to reproduce the Lagrange density (4.1.1), apart from 
an overall factor. 

(ii) Light-cone parameterization 

For our second derivation, we return to (4.3.1)-(4.3.7) and use the remaining reparameterization 
freedom to equate the two X 1 variables with the two <p r variables, renaming both as r\ Also r is 
renamed as t. This parallels the method in Section 2.2(i) . Now in (4.3.1)-(4.3.3) g rs = 5 rs , and 
<9 T X = 0, so that (4.3.3) becomes simply 

G = T = 26 -ifjtp + u 2 (4.3.11) 
m = V6 -\ipV%l) . (4.3.12) 

Therefore the supermembrane Lagrangian (4.3.1) reads 

L M = - J d 2 r|y / 2^-^+ (V0- |^V^) 2 + \^ol x wj • (4.3.13) 

Again a replacement of by ^(1 — e)?/> demonstrates that the integrand coincides with the Lagrange 
density in (4.1.26) (apart from a normalization factor). 

(iii) Further consequences of the supermembrane connection 

Supermembrane dynamics is Poincare invariant in (3+l)-dimensional space-time. This invariance 
is hidden by the choice of light-cone parameterization: only the light-cone subgroup of the Poincare 
group is left as a manifest invariance. This is just the (2 + 1) Galileo group generated by H, P, M, 
B, and N. (The light-cone subgroup of the Poincare group is isomorphic to the Galileo group in one 
lower dimension. [27]) The Poincare generators not included in the above list correspond to Lorentz 
transformations in the "— " direction. We expect therefore that these generators are "dynamical" , 
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that is, hidden and unexpected conserved quantities of our supersymmetric Chaplygin gas, similar 
to the situation with the purely bosonic model. 
One verifies that the following quantities 



are time-independent by virtue of the equations of motion (4.1.8)-(4.1.11), and they supplement 
the Galileo generators to form the full (3 + 1) Poincare algebra, which becomes the super-Poincare 
algebra once the supersymmetry is taken into account. Evidently (4.3.14), (4.3.15) are the super- 
symmetric generalizations of (2.1.14), (2.1.15). 

We see that planar fluid dynamics can be extended to include Grassmann variables, which also 
enter in a supersymmetry-preserving interaction. Since our construction is based on a supermem- 
brane in (3+l)-dimensional space-time, the fluid model is necessarily a planar Chaplygin gas. In 
the next section we shall derive a lineal version of the supersymmetric Chaplygin gas starting from 
a superstring in (2+l)-dimensional space-time. 

4.4 Supersymmetric fluid in (1 + 1) dimensions 

The one-dimensional case is in principle simpler since, in one spatial dimension, the canonical struc- 
ture can be straightforwardly realized. The physical implications of adding Grassmann variables, 
however, are somewhat limited since there is no vorticity and no spin in one space dimension. 

Nevertheless, a supersymmetric version of the lineal Chaplygin gas can be constructed. This 
is achieved along the same lines as in 2+1 dimensions, by considering a superstring moving on 
a plane and again fixing the parametrization invariance. The construction is analogous to what 
has already been done in one higher dimension: the Nambu-Goto action for a supermembrane in 
(3+l)-dimensions gives rise, in a specific parametrization, to a supersymmetric planar Chaplygin 
gas. 

As we shall demonstrate, the supersymmetric extension enjoys the same integrability properties 
as the purely bosonic, lineal Chaplygin gas, as a consequence of the complete integrability for the 
dynamics of the superstring on the plane. 

(i) Superstring formulation 




(4.3.14) 



G 




(4.3.15) 



We begin with the Nambu-Goto superstring in 3-dimensional space-time, 
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I = - J drda (V? - ie ij diX^d^), (AAA) 

where 

g = -det{nfn^ M 4, (4.4.2) 

nf = fi^" - ijrfdoj) . (4.4.3) 

In these expressions u., v are spacetime indices running over 0, 1,2 and i,j are world-sheet indices 
denoting r and a. We now go to the light-cone gauge where we define X = (X° ± X 2 ). X + is 
identified with the timelike parameter r, X~ is renamed 8, and the remaining transverse component 
X 1 is renamed x. We can choose a 2-dimensional Majorana representation for the 7-matrices: 

7 = a 2 , 7 1 = -ia 3 , 7 2 = z'cr 1 , 

such that ip is a real, two-component spinor. A remaining fermionic gauge choice sets 

= 0, 

where 7 ± = ^75(7° ± 7 2 )- Thus ^ is further reduced to a real, one-component Grassmann field. 
Finally we define the complex conjugation of a product of Grassmann fields (^1^2)* — ^1^2 so as t° 
eliminate % from Grassmann bilinears in our final expression. The light-cone gauge-fixed Lagrange 
density becomes: 

C = + (4.4.4) 

where 

9 = (d a x) 2 , (4.4.5) 

u 2 

T = 2d T 9 - (d T x) 2 - 2V2ipd T ip + — , (4.4.6) 

9 

u = d a 9- d T xd a x - V2tpd T ip . (4.4.7) 

In the above equations, d a and d T denote partial derivatives with respect to the spacelike and 
timelike world-sheet coordinates. The canonical momenta 

d(d T x) V r 9 
IT = J¥L = - A [i (44 9) 

d(d T 8) Vr' 1 } 

satisfy the constraint equation 

pd a x + UdJ - V2U^d a ip = (4.4.10) 
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and can be used to recast £ into the form 

C = pd T x + Ud T 9 + ^{p 2 + g) + V2tpd a 4j - V^Utpd^ + u(pd a x + UdJ - V^Utpd^), (4.4.11) 

where u is now a Lagrange multiplier enforcing the constraint. We use the remaining parameteri- 
zation freedom to fix u = and II = —1 and perform a hodographic transformation, interchanging 
independent with dependent variables. The partial derivatives transform by the chain rule: 

d a = (d a x)d x = ^gd x , (4.4.12) 
d T = d t + (d r x)d x = d t + vd x , (4.4.13) 

and the measure transforms with a factor of 1/y 7 ^. Finally, after renaming ^fg as \p2~\f p, we obtain 
the Lagrangian for the Chaplygin "super" gas in (l+l)-dimensions, 

If 1 1 A \/2\ 

L = — J dx {-p(0 - -W>) - -pv 2 - - + , (4.4.14) 

where according to (4.4.8) and (4.4.10) (at u — and II = —1) 

v = p = 9' - ~W . (4.4.15) 

We have used p and v in anticipation of their role as the fluid density and velocity, and we demon- 
strate below that they indeed satisfy appropriate equations of motion. For convenience we have also 
rescaled ip everywhere by a factor of 2~ 3//4 . The Lagrangian (4.4.14) agrees with the limiting case of 
the planar fluid in (4.1.1). We note that as for the planar case, a more straightforward derivation 
leads to the fluid Lagrangian of (4.4.14) with p integrated out. Specifically, if the parameterization 
freedom is used directly to equate the spacelike and timelike coordinates a and r with x and t, we 
obtain 

U = - J dx (^J 29 - i/jip + v 2 - ^ip , (4.4.16) 

where v is defined as in (4.4.15). This form of the Lagrangian can be obtained from (4.4.14) after 
p is eliminated using the equations of motion for 9 and ip, shown below. 

(ii) Supersymmetric Chaplygin gas 
(a) Equations of motion 

The following equations of motion are obtained by variation of the Lagrangian (4.4.14). 

P + d x (pv) = (4.4.17) 



v 2A\ 

ip + (v + l-^)ip' = (4.4.18) 

9 + v9' = -v 2 + A;-—^' (4.4.19) 
2 p 2 2p 

v + vv' = d x [^j (4.4.20) 
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Naturally, there are only three independent equations of motion as (4.4.20) is obtained from (4.4.18), 
(4.4.19) and (4.4.15). Equations (4.4.17) and (4.4.20) are seen to be just the continuity and Euler 
equations for the Chaplygin gas. Note that these do not see the Grassmann variables directly. 

We now pass to the Riemann coordinates, which for this system are (velocity ± sound speed 
V2X/P): 



R± = iv±—) . (4.4.21) 



P 

In terms of the Riemann coordinates, the equations of motion obtain the form 

R± = —R T R'±, (4.4.22) 
ip = -R+tl)', (4.4.23) 
9 = - l -R + R_- l -R + ^' . (4.4.24) 

The equations in (4.4.22) contain the continuity and Euler equations and are known to be integrable. 
It is readily verified that equation (4.4.23) for ip is solved by any function of 

ip = V(R-), (4.4.25) 

and hence the fluid model is completely integrable. That this is the case should come as no surprise 
considering that we began with an integrable world-sheet theory. 

At this point it may seem curiously asymmetric that equation (4.4.23) for the Grassmann field 
should contain the R + Riemann coordinate and not the R- companion coordinate. In fact, the 
reverse would have been the case if the sign of the a/2A term in (4.4.14) had been opposite. The 
entire model is consistent with this substitution, which is just the choice of identifying yj~g with plus 
or minus the sound speed V2X/ p. 

The energy-momentum tensor is constructed from (4.4.14), and its components are 



T 00 = H = \pv 2 + --^^\ (4.4.26) 
T 01 = V = pv, (4.4.27) 



T 10 = ^R + R--^±R + W\ (4.4.28) 
T 11 = pR + R_ . (4.4.29) 



The expected conserved quantities of the system, the generators of the Galileo group, are verified 
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to be time-independent using the equations of motion. We have 

N = J dxp, (4.4.30) 

P = J dxpv, (4.4.31) 

dx(lpv 2 + ^-y^^), (4.4.32) 

B = I dxp(x-vt) = dxxp-tP, (4.4.33) 



H 



Although some generators look purely bosonic, there are still Grassmann fields hidden in v according 
to its definition (4.4.15). 

In going to Riemann coordinates, we can observe a ladder of conserved charges of the form 



J± = / dxpRl . (4.4.34) 



The first few values of n above give 

It = N, (4.4.35) 

if = P±V2\n, (4.4.36) 

1+ = 2H , (4.4.37) 

where Q is used to denote the length of space J dx. (Note that J7, would correspond to the 
Hamiltonian of the theory with v2A replaced by its negative). 

In Section (4.2) we identified two different supersymmetry generators, which correspond in one 
space dimension to the time independent quantities 



Q = j dxpip, (4.4.38) 

Q = J dx p(v - . (4.4.39) 

These are again special cases (n = and n = 1) of a ladder of conserved supercharges described by 

Q n = f dxpR n _^p . (4.4.40) 



We see that the supercharges evaluated on the solution (4.4.25) reproduce the form of the bosonic 
charges (3.1.8). 

Let us observe that there exist further bosonic and fermionic conserved charges. For example, 
one may verify that the bosonic charges 

f R '±\ m 



dxpRl^j (4.4.41) 
dxpRl(^j (4.4.42) 
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are conserved, as are the fermionic charges 

jdxpR n _{^-\ (4.4.43) 

Conserved expressions involving higher derivatives may also be constructed. The conservation of 
these quantities is easily understood when the string world-sheet variables are used. Then the above 
are written as J daR±(d cr R±) m , J dcrR™(i/)d a ip), and j daR^^d^), respectively. Furthermore when 
R± are evaluated on solutions, they become functions of r ± a [see Section 3.3], so that integration 
over u extinguishes the r dependence, leaving constant quantities. 



(b) Canonical structure 

The equations of motion (4.4.17)-(4.4.19) can also be obtained by Poisson bracketing with the 
Hamiltonian (4.4.26) if the following canonical brackets are postulated. 

{9(x),p(y)} = S(x-y) (4.4.44) 
{9{x)^{y)} = -t-8{x-y) (4.4.45) 

U>{x),il){y)} = --S(x-y) (4.4.46) 

where the last bracket, containing Grassmann arguments on both sides is understood to be the 
anti-bracket. With these one verifies that the conserved charges in (4.4.30)-(4.4.33) generate the 
appropriate Galileo symmetry transformations on the dynamical variables p, 9, and ip. Correspond- 
ingly the supercharges (4. 4. 38), (4. 4. 39) generate super transformations. 

5p = 5p = -r)d x {(nl)) (4.4.47) 

89 = ~r]%l) 59 = --r)R+ip (4.4.48) 
Sip = -7] 5ip = -rjipip' - r]R_ (4.4.49) 

which leave the Lagrangian (4.4.14) invariant. The algebra of the bosonic generators reproduces 
the algebra of the (extended) Galileo group. The algebra of the supercharges is 

{f}Q,VQ} = ZVVH, (4.4.50) 

{VQ,VQ} = VVN, (4.4.51) 

{fjQ,r)Q} = fjri(P-V2Xn), (4.4.52) 

{B,Q} = Q. (4.4.53) 
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(c) Additional symmetries of the fluid model 

As mentioned above, since the fluid model descends from the superstring, it should possess an 
enhanced symmetry beyond the Galileo symmetry in (l+l)-dimensions. In fact, the following 
conserved charges effecting time rescaling and space-time mixing are also verified: 



D = tH - J dxpO , (4.4.54) 
G = /d, (««-«»), (4.4.55) 

The Galileo generators supplemented by D and G together satisfy the Lie algebra of the (2+1)- 
dimensional Poincare group, with N, P, and H corresponding to the three translations and with 
B, D and G forming the (2+l)-dimensional Lorentz group 50(2, 1): 

{B,D} = B, {G,B}=D, {D,G} = G, (4.4.56) 

with Casimir 

C = BoG + GoB + DoD . (4.4.57) 

Adjoining the supercharges results in the super- Poincare algebra of (2+l)-dimensions. The Lorentz 
charges do not belong to the infinite towers of constants of motion mentioned earlier. Rather, they 
act as raising and lowering operators. One verifies for the Q n and I+: 

{ J B,/ n +} = -n/ n + _ 1 , {D,I+} = (n-l)I+, {O, /+} = (f - 
{B,Q n } = -nQ nr . 1 , {D,Q n } = {n-\)Q n , {G,Q n } = (f - \)Q n +i- 

(Note that the {B, J^} bracket coincides with {B,2H}, which should equal — 2P according to the 
Galileo algebra. But the above result, viz. —2/7+ gives — 2(P — \/2\Vt). This central addition arises 
from a term of the form J dxdyV2Xx5'(x — y), whose value is ambiguous, depending on the order 
of integration.) The brackets with the I~ do not close, but the I~ can be modified by the addition 
of another tower of constant quantities, namely those of (4.4.42): 

/-=/-- V2\n(n -I) J dxR^ipip' ■ (4.4.59) 

The modified constants obey the same algebra as 

{£,/-} = {Dj-} = (n-l)I- 1 {0,/-} = (f-l)/- +1 . (4.4.60) 

Evidently I+, I~ , and Q n provide irreducible, infinite dimensional representations for 50(2, 1), with 
the Casimir, in adjoint action, taking the form 1(1 + 1), and / = 1 for /+, I~, and / = 1/2 for Q n . 

We inquire about the algebra of the towers of extended charges J+, J~, and Q n . While some 
(bosonic) brackets vanish, others provide new constants of motion like those in (4.4.41)- (4.4.43) 
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and their generalizations with more derivatives. Thus it appears that one is dealing with an open 
(super) algebra. 

A final comment: We have presented supersymmetric versions of fluid dynamical models in two 
and one space dimension. These models enjoy supersymmetry as well as extra "dynamical" symme- 
tries tracing back to their origin in higher dimensional supermembrane models. Other investigations 
of supersymmetric fluids are reported in Ref. [37]. 

There remain some obvious open questions. One is what other fluid interactions can be obtained 
from the rich factory of (super)branes. For example, string theory D-branes have gauge fields living 
on them. Such gauge fields would presumably remain in passing to a fluid model and may thus 
provide a model of (super) magnetohydrodynamics from D-branes. 

Another question is the construction and properties of fluid models with Grassmann variables 
in arbitrary dimensions. Such models would not descend from supermembrane models and, conse- 
quently, would not enjoy supersymmetry or other extended symmetries. Note, nevertheless, that 
Grassmann Gauss potentials ip can be used even in the absence of supersymmetry. For example, 
the above models with the last explicitly fermion- dependent term omitted, possess a conventional, 
bosonic Hamiltonian without supersymmetry, while the Grassmann variables are hidden in v and 
occur only in the canonical 1-form. These models would describe fluids with fundamental spin de- 
grees of freedom and it would be worthwhile to explore their physical properties in this description. 
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5 NON-ABELIAN FLUIDS 



5.1 Introduction 



In this Section we generalize fluid dynamics to systems with non-Abelian charges. We begin with 
some comments on the physical contexts in which such a generalization might be useful and the 
scope and limitations of possible approaches to the problem. 

The quintessential example of a physical system with non-Abelian charges is the quark-gluon 
plasma. High energy collisions of heavy nuclei can produce a plasma state of quarks and gluons. 
This new state of matter has recently been of great interest both theoretically and in experiments at 
the RHIC facility and at CERN. In fact, there is growing evidence that such a state has already been 
achieved at the RHIC facility [38]. In attempting a theoretical description, there are basically two 
approaches that we can use. Since the plasma is at high temperatures, one can argue that the average 
energy per particle is high enough to justify the use of perturbative Quantum Chromodynamics 
by virtue of asymptotic freedom. However, it is known that because of the infrared divergences 
various resummations, such as summing hard thermal loop contributions, have to be done before 
a perturbative expansion with control of the infrared degres of freedom can be set up [39]. One 
has to address also the question of chromomagnetic screening, because unlike the Abelian plasma, 
there can be spatial screening of magnetic type interactions [40]. The expected end result is then 
a good description valid at high temperatures and for plasma states that are not too far from 
equilibrium, since one is perturbing around the equilibrium state. An alternative approach, which 
may be more suitable for nondilute plasmas or for situations far from equilibrium, would be to use 
a fluid mechanical description. 

We begin by observing that many of the general comments given in the introduction, on deriving 
fluid mechanics from an underlying particle theory by statistical averages, will apply in the non- 
Abelian context as well. Specifically for the quark-gluon plasma, some work along these lines was 
done many years ago using single particle kinetic equations [41]. The one-particle kinetic equation 
takes the form 







gQ a F t 



fit' 



+ gf ab cA b u Q 



dP, 



f(X,P,Q)=C(f). 



(5.1.1) 



Here f(X, P, Q) is the one-particle distribution function and C is the collision integral term taking 
account of scatterings of the particles. and are the potential and field for a non-Abelian 
theory based on a gauge group with structure constants j a bc- Here Q a represents the classical color 
charge of the particle. 

It may be interesting to note that, for a collisonless plasma, i.e., with C = 0, the Boltzmann 
equation (5.1.1) is the equation for the distribution function for single particles obeying the standard 
classical equations of motion for non-Abelian particles - the so-called Wong equations [42] - which 
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are 



dX» 



m 



dr 



m 



dr 

dQa 



gQaF-P", 



(5.1.2) 



m 



dr 



gf abc (P»Al)Q c . 



As we shall see shortly, the motion of the color degrees of freedom can also be described in a phase 
space way; the appropriate space is the Lie group modulo the maximal torus. 

The Boltzmann equation (5.1.1) is invariant under gauge transformations in the sense that if 
f(X,P,Q) solves (5.1.1), then so does / (X , P, U^ 1 QU) . As in the Abelian case, one may, in a 



equilibrium distribution; the perturbative corrections can then give the transport coefficients and 
fluid equations of motion [41]. Such an approach suffers from the same limitations mentioned in 
the introduction, namely, that it can only be justified for dilute systems near equilibrium and in 
a semiclassical approximation. However, the fact that the equations for the Abelian fluid have a 
fairly large regime of validity, despite the fact that it can be derived from an underlying particle 
theory in a limited context, then prompts us to ask for an a priori derivation of a non- Abelian 
fluid mechanics, which incorporates the non- Abelian degrees of freedom, coupling to a non- Abelian 
gauge field, etc. This theory may be valid for dense, nonperturbative and nondilute systems. This 
is the goal of this section. In proceeding with the development of such a theory, it is useful to keep 
in mind some guidelines or desirable features. First of all, a canonical or symplectic formulation (at 
least in the conservative limit) is important for quantization, so we should aim for this. At the same 
time, the analysis based on the kinetic equations still remains useful to us as a guide for arriving 
at the equations of interest. Our analysis is based on the paper in Ref [43]. 

5.2 Non- Abelian Euler variables 

A possible form for the non- Abelian Euler fluid variables may be inferred from the single-particle 
equations of motion (5.1.2) by a procedure analogous to what was one in Section 1.1 for the Abelian 
case. 

The single-particle non- Abelian current is defined in terms of the Lagrange variables X M and 



dilute system, seek a perturbative solution of the form / = + gf^ 1 ' + where 



n p is the 



Q, 



■a ■ 



With the parametrization X' 




(5.2.1) 



Pa(t,r) = Q a (t) 5(X(t)-r), (5.2.2) 
Ja(*,r) = Q a (t) X(t) 5(X(t)-r). (5.2.3) 
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This generalizes (1.1.2) and (1.1.6) by inclusion of the dynamical non-Abelian charge Q a which 
satisfies (5.1.2), so that J£ is covariantly conserved. 

(£>„ J^)a = <Va + fabcA"^ = (5.2.4) 

(In the present case, the current is defined without the mass factor and it is normalized by Q a .) 
Passage to the fluid description is achieved as in Section 1.1. In the many-particle case, dynamical 
quantities are decorated with the particle label n, as in (1.1.8)-(1.1.10), which is summed in the 
definition of the current. Then in the continuum fluid limit, n is replaced by the continuous variable 
x and the charge and current densities read 



p a (t,r) = J d A x Q (t,x) 5(X(t,x)-r), (5.2.5) 

J a (*,r) = J d 3 x Qa(t,x) X(t,x) 5(X(t,x) -r), (5.2.6) 

with <2 a (t,x) satisfying 

Q a (t, x) + f abc [cA b (t, X(t, x)) + X(t, x) • A\t, X(i, x))] Q c (t, x) = 0. (5.2.7) 

[Notice that replacing a discrete sum by an integral over x forces Q a (t,x) to be a charge density] 
Observe that, just as in the Abelian case discussed in Section 1.1, the x-integration evaluates x 
at x(t, r), the inverse of X, and the Jacobian factor | det dX 1 / dxi\x=-x is just the Abelian charge 
density p [see (1.1.18), (1.1.19)]. Thus the non-Abelian quantities factorize. 

p a {t,r) = Q a (t,r) p(t,r) (5.2.8) 
Ja(*,r) = Q a (t,r) p(t,r) v(*,r) (5.2.9) 

Equivalent ly 



where 



JZ(t,r) = Q a (t,r)f(t,r), (5.2.10) 

Q a {t,r) = Q Q (t,x)| x=x , (5.2.11) 
pit, r) Q a (t, r) = J dxQ a (t, x) <f(X(t, x) - r). (5.2.12) 

As a consequence of its definition, the Abelian current factor j M = (cp, vp) satisfies its own continuity 
equation (1.1.16). Moreover, differentiating (5.2.12) with respect to time and using (1.1.16) and 
(5.2.7) results in an equation for Q a , 

Q a (t, r) + v(t, r) • VQ a (t, r) = -f abc [cA b {t, r) + v(t, r) • A b (t, r)] Q c (t, r) (5.2.13) 

which can also be written as 

j^D^a = 0. (5.2.14) 
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This is analogous to the Abelian equations (1.2.74). Equations (5.2.13) and (5.2.14) can be un- 
derstood from the fact that the covariantly conserved current (5.2.1), (5.2.4) factorizes according 
to (5.2.10) into a group variable Q a and a conserved Abelian current j M . Consistency of (1.1.16), 
(5.2.4) and (5.2.10) then enforces (5.2.14). Evidently this is a generalization of the particle Wong 
equation (5.1.2); thereofre we shall refer to it as the fluid Wong equation. 

We recognize that the formulas (5.2.5) and (5.2.6) are the non-Abelian version of the Lagrange 
variable-Euler variable correspondence [see (1.1.3) and (1.1.4)]. Also, (5.2.7) is the field generaliza- 
tion of the particle Wong equation, presented in the Lagrange formalism, and (5.2.13) and (5.2.14) 
are the Euler version of the same. The decomposition of the non-Abelian current in (5.2.10) is the 
non-Abelian version of the Eckart decomposition (1.4.6). Indeed, (5.2.10) may be further factored 
as in (1.4.6). 

J£(t, r) = Q a (t, r) n(t, r) u»(t, r) (5.2.15) 

In the remainder of Section 5, we are guided in our construction of a dynamical model for non- 
Abelian fluid mechanics and "color" hydrodynamics by the above properties of the non-Abelian 
current, which follow from the very general arguments, based on a particle picture for the substratum 
of a fluid. In Sidebar G, at the end of this Section, we present a different model, based on a field 
theoretic fluid substratum. 



5.3 Constructing the action 

The equations of motion for the non-Abelian fluid in the Euler formulation include the kinematical 
equations: continuity (5.2.4) and Wong (5.2.14) that are satisfied by the non-Abelian current, which 
is factorized as in (5.2.10)-(5.2.15). Still needed is the Euler force equation, analogous to (1.4.14), 
which specifies the dynamics. We present this by constructing an action whose variation reproduces 
the kinematical equations and gives a model for the dynamical equation. 

The algebra underlying the non-Abelian theory is realized with ant i- Hermit ian generators T a 
satisfying 

[T\ T b ] = f abc T c , (5.3.1) 

and normalized by 



tr T a T b = - 




(5.3.2) 



In a canonical particle theory we expect that the algebra (5.3.1) is reproduced by Poisson brackets 
for corresponding symmetry generators. In a field theory, we expect to find a copy of (5.3.1) at each 
point in space, leading to the Poisson brackets 

{Pa(r), pb(r')} = f abc p c (r) 5(r - r), (5.3.3) 

which generalize (1.2.16). [A common time argument is suppressed.] Upon quantization the brackets 
become commutators and acquire the factor i/h. (Schwinger terms do not spoil the quantum algebra 
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unless there are anomalies in the gauge symmetry [44]. Of course, we assume that the theory is 
anomaly-free.) 

The action which leads to the commutation rules (5.3.3) is the field theoretic version of the 
Kirillov-Kostant form, which for a particle (not a field) reads 

'« = *•/*«*»-'*■ < 5 - 3 - 4 > 

where n is a normalization constant and we have taken the group £77(2) as a concrete example: 
g G £77(2), p = exp(T a ip a ),T a = a a /2i and a a are Pauli matrices. The action (5.3.4) is invariant 
under g —>■ g expT 3 ip, modulo surface terms. Therefore, the theory governed by (5.3.4) is defined on 
the 2-dimensional sphere SU(2)/U(1) = S 2 . The observables are given by q a = —2tr (gT 3 <7 _1 T a ), 
one can show that they obey commutation rules {q a , qb} = e a b c q c , which is the single-particle version 
of (5.3.3). Upon quantization, the quantum Hilbert space will consist of one unitary irreducible 
representation of the group ££7(2) with the highest weight or j-value given by j = |n. The 
symplectic 2-form associated with (5.3.4) is the field of a magnetic monopole on the two-sphere 
£ 2 = SU(2)/U(1). The Dirac quantization rule requires that n be an integer, consistent with j 
being half-integral. Interpreting the single irreducible representation of the group as representing 
the charge degrees of freedom, the action (5.3.4) describes a single particle with ££7(2) non-Abelian 
charges. One can in fact use (5.3.4) as part of an action for the Wong equations of motion [8], [45]. 

More generally, consider a Lie group G with H denoting its Cartan subgroup (or maximal torus). 
The required generalized action is given by 

Iq = ^ n s I dt trK {s) g~ l g (5.3.5) 

s=l 

where g G G and n s are the highest weights defining a unitary irreducible respresentation of G, 
are the diagonal generators of the Cartan subalgebra H of G. The summation in (5.3.5) extends 
up to the rank r of the algebra, but some of the n's could vanish. The action (5.3.5) is invariant 
under g — > gh, h G H and time independent, so that the phase space is G/H, which is known to 
be a Kahler (and symplectic) space. Quantization leads to a finite dimensional Hilbert space which 
carries a unitary irreducible representation of G labelled by the highest weights n s . 

Given this structure, we see that field theoretic generalization, which would give rise to (5.3.3), 
appears as j dt dr Y7 s =i n (s) tr {Kt s \ g~ x g), where now n s and g depend on r, while the remain 
as constant elements of the Cartan subalgebra of the group. Thus we take as the Lagrange density 
for our non-Abelian fluid dynamics the formula 

r 

£ = X^M 2 tlK {s)9~ lD v9 ~ /(n(i),ri(2), . . . , n {r) ) + C gauge . (5.3.6) 

s=l 

Here are a set of Abelian currents; they may be taken to be in the Eckart form j?, = n^u 1 ?^, 
where w^'s are four- velocity vectors, with u^u^^ = 1. As far as the variational problem of this 
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action is concerned, we regard ri( s ) as given by via = yj^j^s)- The space and time 
components of the currents are given by j?. = (cp( s ), p( s )V( s )), p( s ) = ri( a )U® a y In equation (5.3.6) 
7i( s ) are the invariant densities for the diagonal directions of the Lie algebra. 

Comparison with the usual form of the action for an Abelian fluid shows that what we have 
obtained is the non-Abelian analogue of the irrotational part of the flow. In the Abelian case, 
and without the gauge field coupling, equation (5.3.6) entails a single contribution, s = 1, and 
2tr(Kig~ 1 d /1 g) = —d^O, with vanishing vorticity. In the non-Abelian case, the vorticity is still 
nonvanishing. One can easily generalize (5.3.6) to include the other Gaussian components of the 
Clebsch parametrized vector which couples to j?y This gives the Lagrangian 

r 

£ = { 2trK (s)9~ lD fi9 + a M«} _ /( n (i)» n (2)» . . . ,n (r) ) + C gauge (5.3.7) 

8=1 

where a M ( s ) is given by 

<VW = a (s) d nP(s)- (5.3.8) 

For the rank - one group SU (2), with its single Cartan element, the s-sum in (5.3.6) is exhausted 
by a single element. We see that the SU(2) fluid has one component. For higher rank groups, the 
non-Abelian fluid can have up to r components, but fewer - indeed even just a single flow - are 
possible when some of the densities 7i( s ) vanish. Mathematically, single-component fluids are the 
simplest, but physically it is unclear what kinematical regimes of a quark-gluon plasma, for example, 
would admit or even require such a reduction in flows. 

The covariant derivative in (5.3.6), namely, 

Dtf = d^g + A^g, (5.3.9) 

involves a dynamical non-Abelian gauge potential = A^T a whose dynamics is provided by 
£gauge- The first term in £ contains the canonical 1-form for our theory and determines the sym- 
plectic structure ad the canonical brackets. We have added the Hamiltonian density part, the 
function f(nnyri(2),---,nr r )) which describes the fluid dynamics. The theory is invariant under 
gauge transformations with group element U 

9 U ~] 9 (5.3.10) 
A, - U-^A^ + djU. 

where U G G. In Side bar D we show that the canonical structure of our theory leads to the 
charge- density algebra (5.3.3). 

D. Sidebar on the charge density algebra 

The portion of the Lagrange density (5.4.1) that determines the Poisson bracket is 

^canonical = p2trKg' 1 g = p2tiQgg~ l . (D.l) 
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where Q = gKg l . The phase space is parameterized by the scalar density p and parameteris </? a 
specifying the elements g of the group G. Thus the phase space is identified as as 

{set of all maps p (r) = R 3 — > R + , g(r) = R 3 — >■ G} 

With a parametrization of the group element, e.g. g(f) = e TaiPa , one sees that gg~ l has the form 
ip a C a b (<p)T\ where the non-singular matrix C a b (<p) is defined by 

c a MT b = ^g-Hv) • (d.2) 

Thus 

^-canonical = ~P <$aP\Qh = —^aC\pb (D-3) 

We give a new name to the combination C a b pb, 

rr = -c> (d.4) 

The (D.3) reads 

£ canonical II 0a G^""-v 

Consequently, applying the formalism explained in Sidebar A(a), we conclude immediately that IP 
and ip a are canonically conjugate. Moreover, the charge density can be expressed in terms of IP 
and c b a , the inverse to C a b , 

p a = -c\U b . (D.6) 

The non-Abelian charge density p a is a function of (t, r) and for (5.3.3) we need the bracket 
with another density evaluated at (t, r'). (The common t-dependence is suppressed.) Since the 
dependence of c a b on </? involves no spatial derivatives of ip, it is clear that the brackets will be local 
in r — r', just as is the bracket between ip and II. 

{p a (r), Pb (r>)} = (c a a ,^U b '- a^&W-r') 



a dip a , 



„a „b ^ c" „c"-nb 



c c h ,U b - a b }5 (r - r' 



^ a' c c\ 

OPa' 

c a a'C b c i^^p c n - a b ) 5 (r - r') (D.7) 



To evaluate the derivative with respect to ip, return to (D.2) and observe 

dc i" = 9 2tr gg g -i T c" 

dip a i dcp a i dip c i 



-^iik:f l - c ^ c ^ ]T '' (D - 8) 
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The first term in the parentheses is symmetric in (a', c'); when inserted in (D.7) it produces a sym- 
metric contribution in (a, b) and does not contribute when antisymmetrization in (a, b) is effected. 
What is left establishes (5.3.3). 

{p a (r),p b (r')} = ( C > 6 c ,C c ;^„2trT rf yY c V»-a-6)Mr-r') 

= - (2trT a T b T c " p d , - a <-> b\ 6 (r - J) 

= -2trf abd T d T c " Pc „5(r-r') 

= f abc Pc(r)5(r-r') (D.9) 

For simplicity here we examined the single-channel case. The multi-channel case can be treated 
similarly. 

It is instructive to rederive the above result by using the general theory of canonical transfor- 
mations, described in Sidebar A (a), (d). 

From (D.3) we see that the canonical 1-form has components 

a p = 0, 

a Va = -pC\Q h = -C\p h (D.10) 
[These generalize a« in (A.l).] The symplectic 2-form [generalizing in (A. 3)] reads 

•/pp( r > r ') = 0, 
/ Wa (r,r') = -5(v-v')C a b Q b , 
/^(r,r') = 5(r-r')pC a c C b d Q e f cde . (D.ll) 

(The delta- function arises because the 1-form is a function of the coordinates and involves no 
coordinated derivatives.) It follows form (A. 12) that the generator of an infinitesimal canonical 
transformation 

5p = -v p 

5^ a = -v*» (D.12) 

obeys 

J dr'[v^UJ = i;** C\Q b = ^ (D.13a) 

dr'[v p f pipa + U b J = -v? C\ Q b - p C\ C\ Q e f cde = ^- (D.13b) 

Let us now consider the left translation of g: 5g = e a T a g, or equivalently 5ip a = € b c b a , 
fiQa = fabc e b Qc, & pa = fabc ^b Pc- Here € a is a function of r. Thus we have v p = 0, v Va = e b c b a and 
(D.13) are solved by 

G = - I dr P Q a e a = - [ dr e a p a . (D.14) 
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From the general theory (A. 13), we learn that by Poisson bracketing G generates the above trans- 
formation on any function of the phase space variables. Thus 

J dre a (r)p a (r),p fe (r')} = -Mr') 

= /a 6c e 6 (r')p c (r'). (D.15) 

Stripping away the arbitrary function e a (r) reproduces (D.9). 

Next we consider a right translation of g by the Lie algebra element K, present in the canonical 
1-form: 5g = g KX, or equivalently 5ip a = XQbC b a , leading to v p = 0, v Va = —\QbC b a . Here A is 
a function of r. Eqs. (D.13) are now solved by 



G = - J drXp, (D.16) 

when K is normalized to tr K 2 = — |. It then follows from (A. 13), that once A is stripped away, 

{p(r),g(r')} = -5(r-r')gK. (D.17) 

and 

{p(r), p(r')} = 0. (D.18) 



5.4 Equations of motion 

In working out the equations of motion and other consequences of this theory, it is instructive to 
consider first the simpler case of a single flow. This is obtained for G = SU(2), but could also occur 
in higher rank groups. 

(i) Single component flow 

The single flow Lagrangian takes the form 

£ = 3 ntiKg- l D^g - f(n) + £„ , (5.4.1) 
where K = a^/2i. Now we have a single Abelian current j M which can be decomposed as 

f = ( C p, vp) =nu p , n% = 1 . (5.4.2) 
The current J% to which couples is easily worked out from (5.4.1). Upon defining 

gKg- 1 = Q = Q a T a , (5.4.3) 
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we find 

J£ = {c Pa J a )=Q a f = Q a nu\ (5.4.4) 

Pa = PQa- 

Q a may be thought of as the charge of a single particle, with p a = pQ a as the non-Abelian charge 
density at the point r. Notice that this is of the Eckart form where the currents are given by charge 
densities multiplied by the velocity vector. This agrees with the general discussion in Section (5.2), 
where we showed that a particle based model for the fluid leads to the Eckart forms (5.2.8), (5.2.9). 

The gauge invariance of the Lagrangian (5.4.1) shows that the current J^ 1 = J£T a must be 
covariantly conserved, i.e., 

(D^) a = <Va + fabcAlJ? = (5.4.5) 
Further the current j M satisfies an ordinary conservation law, 

O^f = (5.4.6) 

In Sidebar E, we show that both conservation laws are a consequence of invariance of the action 
with respect to variations of the group element g: arbitrary left variations of g lead to covariant 
conservation of J£ (5.4.5) while the particular variation 8g = gKX, with A an arbitrary function 
of space-time, ensures that is conserved as in (5.4.6). We shall see later that these can be 
interpreted in terms of the Wong equations. Notice that these conservation laws lead to the fluid 
Wong equation 

f{D,Q) a = (5.4.7) 

as has already been noted in (5.2.14). 

E. Sidebar on varying the group element g 

We determine the variation of 



I = I dtdrJ^j^txK^g-'D^g, (E.l) 

J s=l 



when g is varied either arbitrarily or in the specific manner 

g- 1 5g = K {s , ) \, (E.2) 

where A is an arbitrary function on space-time. This will provide the needed results (5.4.5) and 
(5.4.6) for the single channel situation, as well as for many channels in (5.4.22) and (5.4.23), below. 

Recall the definitions Qt s \ = gK^g~ x and D^g = d^g + A^g, which implies D^g^ 1 = d^g^ 1 — 
g~ 1 A fJL . First, the variation of g~ 1 D fl g is established. 

5 {g-'D.g) = -g^Sgg^D^g + g^D,5g (E.3) 
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To evaluate the last term, note that D^5g = (gg 1 5g) = (D^g 1 5g + gD^ (g 1 5g). Thus 

5 {g-'D^g) = d, (g^Sg) + [g^D.g, g'Hg] . (E.4) 

Inserting (E.4) into the variation of Jo i n (E.l), integrating by parts, and rearanging the trace with 
K {s) gives 

/ r 
dt dr foifo ^K [s) g- l 5g + f {s) 2tx [g^Dtf, K (s) ] g- l 5g\ . (E.5) 

8=1 

Considering first arbitrary variations: the vanishing of 5Iq requires 

r 

E { d ^s)K(s) + jfa [g-'D.g, K (s) ] ) = (E.6) 

s=l 

or, after sandwiching the above between g . . . g^ 1 , 

r 

E (34)Gw [ D ^g-\Q( S) ]) = o . (e.7) 

s=l 

Finally we verify that 

[D^gg-\Q {S) ] =D^Q {S) , (E.8) 
so that the desired results (5.4.5) and (5.4.22) follow. 

r / r \ 

E ( d ^(s)Q(s) + j(s) D »Q(s)) = £#)GW = DpJ" = (E.9) 

8=1 \S=1 J 

Next we consider the specific variation (E.2) and separate the sum (E.5) into the term s = s' 
and s 7^ s'. After a rearrangement of the last term in (E.5), we get 

SIq = - J dt dr (d,jfc 2trK {sl) K {s/) \ + jfc 2trg- 1 D, l g[K {s '), a) 

+ E { d ^(s) 1trK {s) K {sl) \ + jfc 2tig- 1 D^g[K {s) , K^Xj . (E.10) 

The first commutator vanishes; so does the second when Kr s ) and Kr s ^ commute, i.e. when they 
belong to the Cartan subalgebra. Also 2tiK^K^) = —K?,K? s ,y for s' = s this is constant, while 
for s' s it vanishes when it is arranged that distinct elements of the Cartan algebra are selected. 
Thus for stationary variations j?s must be conserved, and (5.4.6) as well as (5.4.23) are validated. 

It remains to derive the Euler equation. This is accomplished by varying j^; stationary variation 
requires 

2tr [Q{D^g)g- 1 ] = 2tx{K g' 1 D ^g) = -£/'(n) , (5.4.8) 
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which we call the non-Abelian Bernoulli equation. The Euler equation then follows, as in the 
Abelian case, by taking the curl . 

d, {2trQ(D u g)g- 1 ) - d u {2trQ(D,g)g~ 1 ) = d, (%f(nj) - d u (^./») (5.4.9) 

In Sidebar F, we show that manipulating the left side allows rewriting (5.4.9) as 

2tr(D,Q)(D l/ g)g~ 1 + 2trQF Ml , = d, (^f/'(n)) - d v (^/'(n)) • (5.4.10) 

Finally, contracting with j M = nu^ and using (5.4.5) produces the relativistic, non-Abelian Euler 
equation. 

— (u v f(n)) - ndj'(n) = 2trJ^ (5.4.11) 

The left side is of the form of the usual Abelian Euler equation; the right side describes the non- 
Abelian Lorentz force acting on the charged fluid. 

F. Sidebar on manipulating equation (5.4.9) 

Observe that the first term in (5.4.9) equals 

d fl 2trQ(D u g)g~ 1 = 2tr(^D fl Q)(D u g)g- 1 + Q{D^D v g)g- x - Q{D u g)g~\D^g)g-^ . (F.l) 

The first term on the right side is rewritten with the help of (E.8) and combined with the last term, 
leaving 

2tr {Q{D, l D v g)g- 1 - Q{D tl g)g-\D v g)g- 1 ) . 
After antisymmetrization in (/i, is), the left side of (F.l) reads 

2trQ (([£>„, D v \g)g- X - [{D^g' 1 , (D^g' 1 )) = 2tr (QF, U - [Q, (D^g-^D^g' 1 ) . (F.2) 

When (E.8) is used again, (F.2) becomes the left side of (5.4.10). 

The curved space generalization of the action for the Lagrangian (5.3.10) is given by 

1 = J V=vd 4 x (j^trKg^D^g - f{n) J + I„, (5.4.12) 

where n = W ' i^j u n^ u . (We use rj for the metric even in curved space to avoid confusion with g, the 
group element.) The variation of this action with respect to the metric, 51 = — | J ■sJ—T) Sr]^ 
identifying the total energy-momentum tensor as 

gauge i 

= _^ n f^ _ f( n y + ^_ n f( n ) . (5.4.13) 
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We have used equation (5.4.7) to eliminate tiKg~ l D^g in the matter part of the energy-momentum 
tensor 9^ v '. This behaves as the corresponding Abelian expression (1.4.5), except that now there is 
an interaction with a non-Abelian gauge field. 

The divergence of 9^ u entails two independent parts: one proportional to u u and the other 
orthogonal to it, compare (1.4.12). 



u v f'(ri) 
8,9^ = d,{nu»)^±-L + n 
cr 



c 2 A* 



W(n)) - d v f{n) 



(5.4.14) 



The first vanishes by the virtue of (5.4.6) and the rest is evaluated from the Euler equation (5.4.11), 
leaving 

d^ u = 2trJ M F^ , (5.4.15) 
which is canceled by the divergence of the gauge-field energy-momentum tensor. 

deluge = -2trJ^ (5.4.16) 

Thus we have conservation of the total energy-momentum tensor. 

We record the nonrelativistic limit of the Euler equation (5.4.11). For small velocities, we may 
write, as in Section 1.4. 

1 2 

n P "2? PV 

u" « (c,v) (5.4.17) 

Further, we take f(n) to be of the form 

/ = nc 2 + V{n) . (5.4.18) 

With these simplifications, we find that the nonrelativistic limit for the spatial component of (5.4.11) 
gives the Euler equation with a non-Abelian Lorentz force 

v + v ■ Vv = force + Q a E a + - x Q a B a , (5.4.19) 

c 

where force is the pressure force coming from the potential V (and is therefore Abelian in nature), 
while non-Abelian force terms involve the chromoelectric and chromomagnetic fields. 



E\ = cF« , B\ = ~e ijk F^ k (5.4.20) 



It is seen that the non-Abelian fluid moves effectively in a single direction specified by j = pv. 
Nevertheless, it experiences a non-Abelian Lorentz force. 

Eventhough we wrote down the curved space action (5.4.12) primarily to obtain the energy- 
momentum tensor, we note that it may be useful for relativistic astrophysics with the quark-gluon 
plasma as in early universe or perhaps in the interior of neutron stars. 
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(ii) Multi-component flow 

We now return to the more general multi-component action (5.3.6). The current which couples to 
the gauge potential is now 

r 

J* = Y1 Gwifo ' with = aK^g- 1 ■ (5.4.21) 

s=l 

Arbitrary variation of g ensures that (5.4.21) is covariantly conserved, D^J^ = 0, but we also need 
the conservation of individual j?y This is achieved by considering special variations of g of the 
form 5( s )g = gK^X^ s y These variations of g lead to 

d»J(s) = P(s) + V • (v W p W ) = (5.4.22) 

The fluid Wong equation which follows from the conservation of the Abelian and non-Abelian 
currents now reads 

r 

J2j(s)D,Q(s) = 0. (5.4.23) 

s=l 

Varying the individual j?. in (5.3.7) produces the Bernoulli equations 

2trQ {s) (D,g)g~ l = ^f s \ (5.4.24) 



c 

where 

d 

/w = /( n(1) , n (2) , . . . , n (r) ). (5.4.25) 

As in the single channel case, the curl of equation (5.4.24) can be cast in the form 

\ 2 {uUf (s) ) ~ Q v (<)/ (s) ) } = 2tr({D»Q {s) ){D»g)g- 1 + Q {s) F^ . (5.4.26) 
When contracted with = n^u^, this leaves 



c 2 



d, - n {s) d v f^ = j, is) 2tr((D»Q is) )(D»g)g- 1 + Q (s) F^ . (5.4.27) 



However, unlike in the single channel case, the right side does not simplify since j/x(s)Q(s) cannot 
be replaced by J M because the latter requires summing over s. Also the first right-hand term in 
(5.4.27) does not vanish since (5.4.23) requires summation over s. Equations (5.4.27) are the Euler 
equations for the multicomponent non-Abelian fluid. 

The matter part of the energy-momentum tensor is now given as 



W = -<T £ ^)/ (s) - / + E Jf ¥ fl ^)/ (S) • (5-4-28) 



s=l 
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Its divergence is of the form (5.4.14). 



s=l 



U 



(s) 



(5.4.29) 



Using (5.4.22) and (5.4.27), the right side of this equation is evaluated as 

E W.)2tr ({D^Q [s) ){D v g)g- 1 + Q {s) fA. 

s=l ^ ' 

Since now we are summing over all channels, it follows from (5.4.21) and (5.4.23) that, as before, 



d^ v = 2tr JpFi™ . 



(5.4.30) 



Some simplifications which lead to a more transparent physical picture occur if the dynamical 
potential separates 



f(n {1) ,...,n {r) ) = ^2f(s)(n {s )] 

8=1 

f (s) = 4) 



(5.4.31) 



Then the left side of (5.4.27) refers only to variables labeled s, while the right side may be rewritten 
with the help of the generalized Wong equation (5.4.23) to give 



n (s) U (s) 



(5.4.32) 

Thus in the addition to the Lorentz force, there are forces arising from the other channels s' ^ s. 
Note that for separated dynamics (5.4.31), the energy-momentum tensor also separates, 



^ = E*W = E -9" V Mis) - / W (n W )) + 



n is)f(. 



is) 



(5.4.33) 



8=1 



S=l 



but the divergence of each individual does not vanish. This is just as expected since energy can 
now be exchanged between the different channels and with the gauge field; this is also evident from 
the equation of motion (5.4.32). It is clear that this fluid moves with r different velocities V( s ). 

The single- channel Euler equation (5.4.11) is expressed in terms of physically relevant quantities 
(currents, chromomagnetic fields); the many-channel equation (5.4.27) involves, additionally, the 
gauge group element g. One may simplify that equation by going to special gauge, for example 
g = I, so that the right side of (5.4.27) reduces to 

3 [ s) 2ti({D^Q {s) ){D u g)g' 1 + Q [s) F^ = 3 ^2tiK {s] (d,A u - cU M ) (5.4.34) 
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while the Wong equation (5.4.23) becomes 



r 



(5.4.35) 



8=1 



It is interesting that in this gauge the nonlinear terms in F^ v disappear. 

Observe that the inclusion of the a, /^-components of the current, or the use of the Lagrangian 
(5.3.7), does not change the form of the equations of motion for fluids when expressed in terms of 
the velocities and densities. The expressions for these quantities in terms of the group parameters 
and a, (5 will, of course, be altered. 

G. Sidebar on field-based fluid mechanics 

In the body of our review, we presented Eulerian variables, their connection to Lagrangian 
variables and their equations of motion in a picture derived from an underlying physical reality 
composed of point particles, whose discrete distribution is approximated by a continuum. For the 
Abelian case this was developed in Section 1.1. The hallmark of this approach is that the fluid 
current factorizes in an Eckart form. For the Abelian current we have 



which is no restriction at all, but for the non-Abelian current the further factorization is nontrivial, 



Here we present an alternative view of fluid mechanics, which is field-based as opposed to 
particle-based. In the Abelian case, it results in fluid equations that coincide with those of the 
particle-based derivation. This is in keeping with the fact that the Abelian Eckart form is not 
restrictive. However for the non-Abelian situation, the field-based picture results in equations 
which are different and much less elegant than the particle-based ones discussed in the body of the 
text. 

A field-based realization of the Euler equations for an Abelian fluid is provided by the Madelung 
"hydrodynamical" rewriting of the Schrodinger equation [46]. 



(G.l) 



J£ = Qaf = Qa nU». 



(G.2) 



(G.3) 



(We consider only the free equation.) Upon presenting the wave function as 



(G.4) 



we find that the imaginary part of (G.3) results in the continuity equation for j M = (cp, j), where 
the spatial current j is also the quantum current (h/m) Im.ip*'Vip. When j is written as vp, v is 
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identified as V#; the velocity is irrotational. The real part of (G.3) gives the Bernoulli equation, 
with a quantum force derived from 

fr 2 ft 2 CVn) 2 

V = f^Vv/p) 2 = f- 2 ^ (G.5) 
2m M 8m z p 

The Euler equation follows by taking the gradient of the Bernoulli equation. In this way, we arrive 
again at the conventional irrotational fluid. 

The story changes if we start from a non-Abelian Schrodinger equation. Again we consider 
the free case with nonrelativistic kinematics. Thus the equation involves a multi-component wave 
function with 

ft 2 

im = v 2 ^ (G.6) 

2m v ' 

The color degrees of freedom lead to the conserved non-Abelian current. 

Ja = (cPa,Ja) 

p a = i^T a ^, 3 a = —Re^T a W (G.7) 

m 

The singlet current j M = (cp, j) is also conserved. 

For definiteness and simplicity, we shall henceforth assume that the group is SU(2) and that the 
representation is the fundamental one: T a = a a /2i, {T a ,T b } = —5 ao /2. We shall also set the mass 
m and Planck's constant h to unity [47] . The non-Abelian analogue of the Madelung decomposition 
(G.4) is 

V = Jpgu, (G.8) 

where p is the scalar ^^r ) ^ is a group element, and u is a constant vector that points in a fixed 
direction [e.g., for SU{2) the two-component spinor u could be taken as u% = 1, u-i = 0, then 
iv)T a u = #j/2]. The singlet density is p, while the singlet current j = im\I/*V\I/ is 

j = vp, \ = -iu ] g^V g u . (G.9) 

With the decomposition (G.8), the color density (G.7) becomes 

Pa = QaP, Q a = m ] g- 1 T a gu = iR ab v?T b u= R ab t b /2, (G.10) 

where R ab is in the adjoint representation of the group and the unit vector t a is defined as t a /2 = 
iwT a u. On the other hand, the color current reads 

Ja = \p R ab ^ (T b g-'Vg + g- l Vg T b ) u , (G.ll) 



which with the introduction of 



g- l Vg = -2v a T a , (G.12) 
v = v a t a , (G.13) 
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may be presented as 



Ja = ^Rob • (G.14) 



Unlike the Abelian model, the vorticity is nonvanishing. 

V x v a = e abc w b x v c (G.15) 

A difference between the Madehmg approach and the previous particle based one is that the 
color current is not proportional to the singlet current. Equation (G.14) may be decomposed as 

Ja = QaPV + ^R ah v\ (G.16) 

where the "orthogonal" velocity v" is defined as 

= (5 ab - t a t b ) v b . (G.17) 

Equation (G.15) shows that color current possesses components that are orthogonal to the singlet 
current. 

In a Postscript at the end of this Sidebar, we derive for the SU(2) case the decomposition of 
the Schrodinger equation with the parametrization (G.8). Two equations emerge: one regains the 
conservation of the Abelian current and the other is the "Bernoulli" equation. 



2 r^i 1 

t° + -V- ( P e abc vH c ) (G.18) 



It is further verified that the covariant conservation of the color current is a consequence of the 
Abelian continuity equation and (G.18). However, there is no Wong equation because the color 
current is not proportional to the conserved singlet current. Finally, using the identity, which 
follows from the definition (G.12), 

v a = -\v (rt) a + e ab( y b (g- l g) c , (G.19) 

one can deduce an Euler equation for v a from (G.18). 
We record the energy and momentum density 

£ = -V^ ■ = -pv a ■ v a + Vp ' Vp (G.20) 
2 2 8p 

V = l - (V^% - ^ f W) = vp (G.21) 

Both parallel and orthogonal components of the velocitiy contribute to the energy density but only 
the parallel component v contributes to the momentum density. It is clear that within the present 
approach the fluid color flows in every direction in the group space, but the mass density is carried 
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by the unique velocity v. This is in contrast to our previous approach where all motion is in a 
single direction or at most in the directions of the Cartan elements of the Lie algebra (Section 5.4). 

The difference between the two approaches is best seen from a comparison of Lagrangians. For 
the color Schrodinger theory in the Madelung representation 

= \ - *t *) - i W> • V* (G.22) 

= vuVju- ipv".v°- (G.23) 

With u®u* = 1/2 — 2iK, the free part of the above reads 

^Schrodinger = P ^Kg^g - lpv a ■ V a . (G.24) 

On the other hand, the free part of the Lagrange density (5.4.1) in the nonrelativistic limit, with 
f(n) given by (5.4.18), is 

C° = p 2ti( [ Kg~ 1 g + Kw • g~ l Vg\ - \/ p\c 2 - v 2 ) , 

« p 2tr {Kg-'g + Kv ■ g^Vg^j - pc 2 + ^pv 2 , (G.25) 

= p 2tr Kg~ l g - ^pv 2 - pc 2 , 

where we have used v = — 2tr Kg~ l Vg, which follows upon the variation of v, in the next-to-last 
equality above. Thus the canonical 1-form is the same for both models while the difference resides 
in the velocity dependence of their respective Hamiltonians. Only the singlet v enters (G.25) while 
the Madelung construction uses the group vector v a . 

Finally, note that while the Euler equation, which emerges when (G.18) and (G.19) are combined, 
intricately couples all directions of the fluid velocity v a , it does admit the simple solution v a = vt a , 
with v obeying the Abelian equations that arise from (G.3)-(G.4). 

Postscript: 

When (G.8) is inserted into (G.6), and use is made of the definition (G.12), we find in the SU(2) 
case 

Upu + tp{g- 1 g) a T a u = ~^pV 2 ^pu + V(pv a )T a M + ipv a ■ v a M . (G.26) 

Next (G.26) is premultiplied by u\ where it implies 

ip + p(g^g) a t a = -^pV 2 yfp - iV (pv a t a ) + pv a • v a . (G.27) 

The imaginary part reproduces the continuity equation for the singlet current, while the real part 
gives 

{g~ 1 g) a T a = — ^V 2 v/p + v a - v a . (G.28) 
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To obtain further information, we premultiply (G.26) with u^T b . This gives 

pt'-ipig^gY + pig^gYe^r = iJpV 2 JpT b - V (pv fe ) - l V ■ {pv a ) e hac T c -ipv a -v a T h . (G.29) 

The imaginary part gives (G.18) while the real part is identically satisfied by virtue of the Abelian 
continuity equation and (G.18). 
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6 NONCOMMUTATIVE FLUIDS 



Most of the fluids examined so far arise out of underlying particle systems. [The sole exception is 
the Madelung-like construction based on an underlying field - the Schrodinger field (both Abelian 
and non- Abelian) - discussed in Sidebar G] . This becames manifest in the Lagrange description, in 
which the coordinates of the underlying particle substratum are explicitly involved. The transition 
to the Euler description, then, allows us to express the system in terms of purely fluid quantities, 
namely density and velocity. 

It is possible to consider fluids whose Euler description does not derive from an explicit under- 
lying Lagrangian description. Although such fluids would look conventional in terms of density and 
velocity, they would nevertheless be effective descriptions of exotic underlying theories. 

A concrete and interesting realization of this is noncommutative fluids. In these, the fundamental 
degrees of freedom represent "particles" on a noncommutative space, in which coordinates become 
non commuting operators and the notion of points breaks down. 

Noncommutative spaces were introduced by Heisenberg to ameliorate ultraviolet infinites in 
quantum field theory. They have been considered in mathematics and they arise as particular 
projections of quantum systems in an intense magnetic field (see Section 7.2) and also as special 
limits of string theory. Gauge theory on these spaces becomes particularly attractive, since it fuses 
spatial and internal degrees of freedom into one coherent formalism. As we shall demonstrate, such 
theories can be viewed as noncommutative fluids and the transitition to the Euler description will 
emerge as the transformation mapping them to equivalent commutative effective theories, known 
to the string literature as the Seiberg-Witten map. The presentation follows Ref. [48]. 

6.1 Review of noncommutative spaces 

Noncommutative spaces are described in terms of coordinates that are noncommuting operators. 
In the simplest realization, 'flat' noncommuting coordinates are characterized by a constant, anti- 
symmetric tensor Q l K 



[x\ x 



■?] = i6 lJ 



(6.1.1) 



We must specify whether possesses an inverse u>. 



(6.1.2) 



An inverse can exist in even dimensions, provided 6 is nonsingular, but oj will not exist in odd 
dimensions, where the antisymmetric 6 always possesses a zero mode. We shall assume the generic 
situation: nondegenerate with no zero modes in even dimensions, where we give a Euclidean 
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treatment: all coordinates are spatial; while in odd dimensions there is one zero mode, chosen in 
the time direction. Therefore, in all dimensions only spatial variables are noncommuting. 

Fields f are defined as arbitrary functions of the operator coordinates x l (and t in odd dimen- 
sions). They can be expressed in a Taylor expansion as sums of monomials involving (ordered) 
products of the coordinates. They are, effectively, themselves operators on the same footing as the 
x % (explicitly depending on the commuting ccordinate t in odd dimensions). The derivative of a 
function <9jf is defined, as usual, by eliminating x % once from each place in monomials in which it 
appears. Equivalently, this can be achieved through the adjoint action of the x J themselves: 

d i f=-iui j [x 3 ,f\, (6.1.3) 

which clearly has the desired effect on monomials. The volume integral over the full noncommutative 
space, on the other hand, can be expressed as the trace of the corresponding operator over a 
Heisenberg-like Hilbert space on which the x % act. 

J f = v / det(27r6>)Trf (6.1.4) 

Products of fields in the above space are defined as the usual ordered products of the corresponding 
operators. We may trade the above noncommutative description for one involving usual, commut- 
ing functions by introducing a new, nonlocal, noncommutative product between functions, called 
star-product. First, order all monomials in a fully symmetric way in the x % ("Weyl ordering"), 
which can always be achieved by using (6.1.1). This creates a one-to-one correspondence between 
noncommutative operators and commutative functions. The derivative and volume integral of the 
noncommutative function as defined in (6.1.3) and (6.1.4) map, in fact, to the usual derivative and 
integral of the commutative function f. The product of two operators f and q, on the other hand, 
upon Weyl-reordering, corresponds to a new function, called the star-product of the corresponding 
commutative functions f and q and denoted by f*g, is defined by 

(f*9)(x) = exp±0U^£ i f(x)g(y) \ x=y . (6.1.5) 

We may further define the ^-commutator of two functions as their antisymmetrized ^-product: 

[f,fl]* = f*0-fl*f (6.1.6) 

The noncommuting coordinates x l , in particular, map to the usual commuting coordinates and their 
^-commutator reproduces the noncommutative space relations (6.1.1). 

\x\x%=iQ ij (6.1.7) 

The two formulations, operator and star, are obviously equivalent and will be used interchangeably. 
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For future use let us record here the formulas of noncommuting, electrodynamics in 4-dimensional 
space. The vector potential is a function of the noncommuting coordinates x M , and it undergoes 
gauge transformations, which infinitesimally read 

SA^ = <9 M A — i [Ap, A]* = A. (6.1.8) 

Here A is the gauge function, also depending on the non commuting coordinates. A field strength 
F^v is defined so that it transform covariantly: under gauge transformations 

6F lu , = -i[F ia ,,\] m . (6.1.9) 

The formula for F^ v therefore reads 

F^ = d^K -dvAp- A v \. , (6.1.10) 
and the equations satisfied by the F^ are 

le^DpF^ = 0, (6.1.11) 
D^F^ = J v , (6.1.12) 

where J u is a noncommuting source current. Note that even though we are dealing with electrody- 
namics, the noncommutativity requires a non-Abelian, Yang-Mills-like structure for the relevanat 
expressions. 

The noncommutative space structure (6.1.1) remains invariant under a group of transformations. 
Subjecting the spatial coordinates to an infinitesimal coordinate transformation 

5x = -f(x), (6.1.13) 

for some operator functions f l , and requiring that (6.1.1) remain unchanged results in the condition 

-[f(x),x^]-[x 4 ,^(x)] = 0, (6.1.14) 
which in turn implies by (6.1.1) that 

-d k f(x)0^ - d k fi(x)6 ik = . (6.1.15) 
The left side is recognized as the Lie derivative of a contravariant tensor, 

L f0 ij = f k d k 6 ij - d k fe kj - d k f j 6 ik , (6.1.16) 

with the first term on the right vanishing since is constant. So the noncommutative algebra (6.1.1) 
is preserved by those coordinate transformations that leave 6 invariant: Lf6 = 0. 
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To solve for f in even dimensions, we define 

/' (> ;j !l, ('-J 1 -V). (6.1.17) 

This entails no loss of generality, because is nonsingular (by hypothesis). Then (6.1.15) becomes 

e u d kg ^ + e^d k9i e ik = o . (e.i.isa) 

Because is nonsingular and antisymmetric, this implies 

dk9e-d e g k = 0, (6.1.18b) 

or 

9i = . (6.1.19) 

Thus we have 

f = e ij d j( f) (6.1.20) 

for the coordinate transformations (in even dimensions) that leave 6 invariant. Since 

V-f = (6.1.21) 

the transformations are volume preserving; the Jacobian of the finite diffeomorphism is unity. How- 
ever, except in two dimensions, these are not the most general volume-preserving transformations. 
Nevertheless, they form a group: the Lie bracket of two transformations like (6.1.15), f\ = O^dj^i 
and fl = l Wj<f>2, takes the same form, 9' l Wj(9 dk^id^)- The group is the symplectic subgroup 
of volume-preserving diffeomorphisms that also preserve O 1 ^ . In two dimensions, where we can set 
8^ = 9e l \ the above transformations exhaust all the area-preserving transformations. 

In odd dimensions, where (by assumption) 6 possesses a single zero mode, for definiteness we 
orient the coordinates so that the zero mode lies in the first direction (labeled — > time) and 6, 
confined to the remaining (spatial) dimensions, is nonsingular. 

9^=(° °.) = l,...,2n) 

V0 l V 1 1 (6.1.22) 

9 ij oj jk = 4 

The infinitesimal diffeomorphisms that preserve 6 are 

fit) 

f lt ={ ■ d (6.1.23) 
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These still form a group. Two transformations, (/i,0i) and (f2,4>2), possess a Lie bracket of the 
same form (6.1.23), with (/2<9t/i — j\d t fi-, /2<9t0i — fidt<t>2 + Q kl dk4>idi(f)2) ■ But the space-time volume 
is not preserved: <9 M / M 7^ 0. (Of course, at fixed time, the spatial volume is preserved.) 

Unit-Jacobian diffeomorphisms also leave invariant the equations for an ideal fluid, in the La- 
grange formulation of fluid mechanics, and in particular a planar (two dimensional) fluid supports 
area-preserving diffeomorphisms.; see Section 1.2 (ii). This coincidence of invariance suggests that 
other aspects of noncommutativity possess analogs in the theory of fluids, whose familiar features 
can therefore clarify some obscurities of noncommutativity. (A similar point of view concerning the 
quantum Hall effect was taken in Ref. [49]) We shall explore this connection and shall demonstrate 
that the Seiberg-Witten map [50] between noncommuting and commuting gauge fields corresponds 
to the mapping between the Lagrange and Euler formulations of fluid mechanics. We shall obtain a 
simple derivation of the explicit "solution" to the Seiberg-Witten map in even dimensions [51] and 
will extend it to odd dimensions. 

The two formulations of fluid dynamics (Lagrange and Euler) can be put in the proper context 
for the noncommutative space setting. The natural Poisson (commutator) structure, present in the 
Lagrange description of a fluid, and the possibility of introducing a vector potential to describe the 
evolution of comoving coordinates, will be recognized as classical precursors of analogous noncom- 
muting entities. Within this framework, we shall show how noncommuting gauge fields respond to 
coordinate transformations, generalizing previously established results. [52] 

As explained in Section 1.1, the Lagrange description uses the coordinates of the particles 
comprising the fluid: X(t, x). These are labeled by the comoving ccordinates x, which are the coor- 
dinates of some initial reference configuration, e.g., X(0,x) = x, see (1.1.11). We may parameterize 
the evolution of X by defining 

X%x) = x l + 9 ij Aj(t,x) (6.1.24) 

which looses no generality provided 6 is nonsingular. As will be seen below, A behaves as a 
noncommuting, Abelian vector potential. 

6.2 Noncommutative Gauge theory 
(i) Commuting theory with Poisson structure 

We introduce into the Lagrange fluid description the (nonsingular) antisymmetric tensor 6. This 
allows for a natural definition of a Poisson bracket, which may be viewed as a classical precursor of 
the noncommutativity of coordinates. We define the bracket by 



so that 



(6.2.1) 
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{x\x j } = 9 ij . (6.2.2) 
It follows from the definition (6.1.24) that 

{X\ X j } = 9 ij + 9 ik e ji F M (6.2.3) 

with 

F ij = ^ i A j --£-A i + {A i ,A j }. (6.2.4) 

It is seen that the structure of the gauge field F is as in a noncommuting theory, with the Poisson 
bracket replacing the * commutator of two potentials A, compare (6.1.10). Also, in the limit that 
the deviation of X from the reference configuration x is small, that is, for small A, we recover a 
conventional Abelian gauge field. 

The above formulas are understood to hold either in even dimensions for a purely spatial Eu- 
clidean formulation (there is no time variable) or in odd-dimensional space-time for spatial compo- 
nents. (X and x are spatial vectors, without time components.) 

(ii) Coordinate transformations in the commuting theory (even dimensions) 

In even dimensions, the ^-preserving transverse diffeomorphism, which also implements the repa- 
rameterization symmetry of the Lagrange fluid, acts on X through the bracket according to (1.2.24) 
and (6.1.20) as 

5,X = f.VX = *«f§|£. (6.2.5) 
This may be presented with help of the bracket defined in (6.2.1). 

Wx)^«^ = { X(xW(x)} (6.2.6) 

Because 5X compares the transformed and untransformed X at the same argument, 5Ai = UijSX^ 
and the volume-preserving diffeomorphism (6.2.5), (6.2.6) induces a gauge transformation on A: 
[Compare (6.1.9), (6.1.10)]. 



^A(x) = V0(x) + {A(x), 0(x)} = D0 (6.2.7a) 
^4(x) = {4(x),0(x)} (6.2.7b) 

We see that the dynamically sterile relabeling diffeomorphism of the parameters in the Lagrange 
fluid leads to an equally sterile gauge transformation, under which X and F transform covariantly, 
as in (6.2.6), (6.2.7b). 
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Next we consider a diffeomorphism of the target space. 

S f X = -f(X) (6.2.8) 

In contrast to the previous relabelings, this transformation is dynamical, deforming the fluid con- 
figuration. Quantities 

C„(X) = ^e ilh ... injn {X i \X»} ■ ■ ■ {X*",X*} , (6.2.9) 

which are defined in d = 2n dimensions, respond to the transformation (6.2.8) in a noteworthy 
fashion. One verifies that 

8 f C n (X) = -V • f(X)C n (X) , (6.2.10) 

so that transverse (volume-preserving) target-space diffeomorphisms leave C n invariant. Eq. (6.2.10) 
is most easily established by recognizing that 

BX i 

CJX) = Pfaff{X i , X j } = det 1/2 {X\ X j \ = det 1/2 det — — . (6.2.11) 

ox 3 

The significance of these transformations is evident from (1.1.19), which shows that l/p(r) = 
(X)| , . when det 1 ^ 2 6 is identifiec 
target space diffeomorphisms becomes 



C„(X)| _ when det 1 ^ 2 ^ is identified with 1/po- The transformation law for p under transverse 



V(r)=f-Vp(r). (6.2.12) 

It follows that this transformation leaves invariant all terms in the Lagrangian that depend only 
on p [like the potential in (1.2.30)]. 

When we restrict the transverse, target-space diffeomorphisms to those that also leave 8 invari- 
ant, i.e., (6.1.20) (of course in two dimensions this is not a restriction), further quantities are left 
invariant. These are constructed as in (6.2.9), but with any number of brackets {X*,X ? } replaced 
by 6 l K 

It is interesting to combine the diffeomorphism of the parameter space with that of the target 
space, for a simultaneous transformation on both spaces. To this end we chose the form of the target 
space transformation to coincide with that of the reparameterization/relabeling transformation. 

f (X) = ^'^r • ( 6 - 2 - 13 ) 
As we shall show below, this results in a gauge-covariant coordinate transformation on the vector 
potential A, once a further gauge transformation is carried out. Thus we consider A = 5$ + Sf, 

AX 1 = {X\ 0(x)} - 0*^® • (6.2.14) 
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[Note that any deviation of f l (X) from 6^ <90(X) / dX^ may be attributed to 0, and can be removed 
by a further gauge transformation.] However, covariance is not preserved in (6.2.14): X on the left 
is covariant, but on the right in the Poisson bracket there occurs 0(x), which is not covariant. The 
defect may be remedied by combining AX with a further gauge transformation, 

5 gaugc X = {X,<KX)-0(x)}, (6.2.15) 
so that in A + S gauge = 6 we have a covariant transformation rule: 

^ = {X\0(X)}-^^p, (6.2.16) 
which in turn implies that A transforms as 

SAi = u l3 {X\ 0(X)} - . (6.2.17) 



To recognize this transformation more clearly, we present it as 

90(X) _ d(j)(y 
dX k dX 



M^M*',**}*^)-^), (6.2.18a) 



and use (6.2.3) to find 

SAi = OMfJ^Q = f k (X)F ki . (6.2.18b) 



Note that in the final expression (6.2.18b) the response of A is entirely covariant: it involves the 
covariant curvature F and the diffeomorphism function f evaluated on the covariant argument X. 
This expression is precisely the gauge-covariant coordinate transformation. (This was derived, by 
a somewhat different method [52].) 



(iii) Coordinate transformations in the noncommuting theory with *-products 
(even dimensions) 

The above development may be taken over directly into a noncommutative field theory by replacing 
Poisson brackets by — i times *- commutators, so that (6.2.2) goes over into (6.1.1). Eq. (6.1.24) 
remains and (6.2.3), (6.2.4) become 

[X\ X% = i0 ij + i6 ik 6 ji F ke (6.2.19) 

with F given by (6.1.10). The covariant transformation rules (6.2.16) and (6.2.17) may be used in 
the noncommutative context, provided a sensible ordering prescription is set for 0(X). This we do 
as follows. Define 

$ = / dx0(X) (6.2.20a) 
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where 4>(X) is a series of (star) powers of X: 

0(X) = c + aX* + \c ii X i * X j + IcijkX' *X j *X k + ---- (6.2.20b) 

[We are not concerned about convergence of the integral (6.2.20a), since we are interested in local 
quantities like (6.2.20b) or (6.2.23) below.] The integration over x (the argument of X) ensures 
that $ is invariant (in an operator formalism the integral becomes the trace of the operators). 
The c-coefficients in (6.2.20b) are required to be invariant agaisnt cyclic index shuffling (so that 
$ and 4> possess the same number of free parameters). Also we require <fi to be Hermitian. [This 
ensures, e.g., that Cjj is real symmetric; that Recife is entirely symmetric and that Im is entirely 
antisymmetric (which is impossible in two dimensions).] Then (6.2.16) and (6.2.17) become 

5X l = -i[X\ 0(X)]* - 0V^L (6.2.21) 

(5$ 

5 A = -iUiAX'^iX)]* - (6.2.22) 

where now the last entries employ a functional derivative: 

5$ 



d + djX 3 + c ijk Xi *X k + ---- (6.2.23) 

In two dimensions, the ordering prescription (6.2.20) and its consequence (6.2.23) preserve the in- 
variance of the [X 1 , X 7 ]* commutator against the target space diffeomorphism [last term in (6.2.21)]. 
Thereby a property of the classical Poisson bracket [c.f. (6.2.10) at n — 1] is maintained in the non- 
commuting theory. 

With 0(X) at most quadratic in X (f at most linear), one readily verifies the result inRef. [52]. 

8A, = l{r(X)^F Jl + F Jl ^r(X)} (6.2.24) 

But with more general (f containing quadratic and higher powers) there arise further reordering 
terms. 

(iv) Coordinate transformations in commuting and noncommuting theories 
(odd dimensions) 

In odd dimensions, with the ^-preserving transformation function given by (6.1.23), the relabeling 
transformation on the base space is 

fyX(t,x) = d^— 0(*,x)— X(t,x) + /(t)*X(t,x) 

dx> dx l (6.2.25) 

= {X(t,x),0(t,x)} + /(t)X(t,x) . 

The fluid coordinate X has components only in the spatial directions. Here the Poisson bracket is 
defined with the nonsingular 6^ . 
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For the target space diffeomorphism we again take the formula (6.2.13), so that the combined, 
noncovariant transformation A = 5<f, + 8f reads 

AX* = (X\<Kt,x)} + f(t)X%x) -^^1^1 . (6.2.26) 
This is modified by the gauge transformation 

5 gaugc X = (X,0(t,X) - 0(t,x)} - (X,/(t)i (t,x)} (6.2.27) 
resulting in the covariant transformation A + <5 ga ugc = S. 

SX* = {X\ <j>(t,X)} - ^'^|^ + f(t)DX* (6.2.28) 

Here DX l = X 1 + {A ,X 1 }, where A is a connection introduced to render the time derivative 
covariant against time- dependent gauge transformations, generated by 0. This is achieved when 
the gauge transformation law for A is 

^i o = 0+{i o ,0} . (6.2.29) 

The spatial components of the vector potential are introduced as before in (6.1.24), and DX l 
becomes 

DX i = e tj ^Aj - djA + {i , ij}) = ij F Oj • (6.2.30) 
The covariant transformation law of A follows from (6.2.19), (6.2.28), and (6.2.30): 

5A i = u ij {X j ,<j)(t,X)} ~ ^; X) + cju f(t)DX* 

dX (6.2.31) 

= f s (t,X)Fji + f(t)F 0i = nt,X)F^ . 
It remains to fix the transformation law of Aq. This requires specifying S^Aq. Since 



- d<j>(t,x) 

f { = dX^ (6.2.32) 



it is natural to take 



<J f i o = -0(t,X) (6.2.33) 
(The time derivative acts on the first argument only.) Thus we have from (6.2.29) and (6.2.33) 

Ai = 0(t, x) + {A , </>(t, x)} - 0(t, X) . (6.2.34) 
After adding to this a gauge transformation generated by (f)(t, X) — x) we are left with 
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8A = ?*^x* + {A Mt,x)} 

86 (6 ' 2 ' 35) 
= qY^X 1 = f%X)F i0 = F(t,X)F» . 

Eqs. (6.2.31) and (6.2.35) coincide with the formula obtained in a conventional commuting gauge 
theory [53]. 

Similar results follow within the noncommuting formalism, once the now familiar ordering 
prescription is given for 0(t, X) and $ = j dx0(t,X). In the noncommutative theory (6.2.31) 
and (6.2.35) are regained, up to reordering terms. 

6.3 Seiberg-Witten map = Euler fluid-Lagrange fluid map 

By considering various limits within string theory, Seiberg and Witten found that noncommuting 
fields can be mapped onto non-local functions of commuting fields [50]. The principle behind the 
mapping can be stated as a requirement of stability agaisnt gauge transformations in the following 
sense. Consider the noncommuting gauge potential to be a functional of the commuting gauge 
potential A^ and of 9. It is then required that a commuting gauge transformation performed on 
the commuting gauge potential with commuting gauge function A : A^ — > A A = + <9 M A, can 
be equivalently achieved by performing noncommuting gauge transformation on the noncommuting 
gauge potential with noncommuting gauge function A : A^ —> A A = (e* A ) * [A^ + id^} * (e iA ) _1 . 
Thus 

i A (A M ) = i M (A M + <9 M A) (6.3.1a) 
This in turn implies a differential equation in 9 for A^. 

BA 1 

^ = --{A*, dp A, + F^}t - (a <- 0) (6.3.1b) 

with "initial" condition A^ \e=o= A^. { , }+ denotes the *-anticommutator. One can solve (6.3.1b) 
order-by-order in 9, with the lowest order solution being 

A ll = A ll - \9 afi A a (d p A^ + Fp,j) + • • • (6.3.2) 

but in general, the equations for different a and (3 contained in (6.3.1b) are not integrable, except 
in two dimensions where 9 af3 involves a single quantity 9 a/3 = e af3 9. 

Aside from its intrinsic mathematical interest as providing a connection between commuting and 
noncommuting fields, the Seiberg-Witten map also serves a practical purpose. The noncommuting 
field strength is not gauge invariant, rather, just as in Yang-Mills theory, it is gauge covariant: 
Ffu, — > (e lA ) * F^v * (e tA )~ 1 . But unlike in Yang-Mills theory, there are no local gauge invariant 
quantities. To obtain a gauge invariant result one must integrate over x M , or (equivalently) consider 
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the trace over the Hilbert space on which the operators x M act. However, if one wishes to compare 
predictions of noncommuting electrodynamics with those of ordinary commuting electrodynamics 
(to set limits on the amount of noncommutativity in Nature) non local quantities are not useful, 
because the physical content of ordinary electromagnetism is expressed by local quantities (waves, 
energy and momentum densities, etc.). Here the Seiberg-Witten map provides a resolution: map 
the noncommuting theory onto a commuting one, from which local gauge invariant quantities may 
be extracted [54]. 

We now show that the inverse Seiberg-Witten map is equivalent to the map between Lagrange 
and Euler descriptions for fluids. The argument is first presented in two Euclidean dimensions and 
(2+l)-dimensional space-time, where (6.3.1b) is integrable (because it involves a single derivative 
variable). Then the argument is taken to higher dimensions. 

(i) Seiberg-Witten map in (2) and (2+1) dimensions 

To construct the Seiberg-Witten map in two Euclidean dimensions, we (temporarily) introduce a 
time dependence in the fluid variables (but not into the diffeomorphism functions - only spatial 
variables are transformed) and observe that (cp, vp) form a conserved 3-vector j a [also true in the 
noncommuting theory when an ordered definition for 5(X(£,x) — r) is given - this will be provided 
below]. Therefore, the dual of j a , e^ ucc j a , satisfies a Bianchi identity and can be presented as the 
curl of a potential, apart from additive and multiplicative constants. 



tfivaj 01 oc F^ v + constant (6.3.3) 
Ffiu = d„A v - d u A^ (6.3.4) 

Note j a , F^v, are ordinary functions, even in the noncommuting setting, since the noncommuting 
variables X are integrands (in an operator formalism, their trace is involved). In particular, the 
spatial tensor is determined by p. 

— A^r) - g-Mv) = F t3 (r) = - £ij (p - p ) = -e ijPo ( J dx5(X(x) - r) - l) (6.3.5) 

(The time dependence is now suppressed.) X contains A, as in (6.1.24). Since X is (noncommuting) 
gauge covariant, the integral in (6.3.5) is (noncommuting) gauge invariant. Therefore, (6.3.5) serves 
to define an (inverse) Seiberg-Witten map between the noncommuting (hatted) and commuting 
(unhatted) variables. The additive (e^po) and multiplicative (—1) constants are fixed by requiring 
agreement at small A. It still remains to give a proper ordering to the 5-function containing X. 
This we do by a Fourier transform prescription 

J dr e*- r i^(r) = -e ijPo J dx (ef x « - e* k ' x ) , (6.3.6) 
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and the ordering (Weyl ordering) is defined by the expansion of the exponential in (* product) 
powers: 4 k ' x = l + z'k-X-|k-X*k-X + -- -. 

When the exponential e^ kX is written explicitly in terms of A: exp^i(kiX l + 9kiS^ Aj), factoring 
the exponential into e* kx times another factor involves the Baker-Hausdorff lemma, (because the 
individual terms in the exponent do not * commute). This leads to an open Wilson line integral 
[55]. In that form (6.3.6) is seen to coincide with the known solution to the Seiberg-Witten map 
[51], which is now also recognized as nothing but an instance of the Lagrange — > Euler map of fluid 
mechanics. 

To construct the Seiberg-Witten map in (2+l)-dimensional space-time we consider the conserved 
current, defined in (1.1.13) and (1.1.14), except that now the time dependence is retained throughout 
and the derivative is gauged with A Q : 

j(t,r) = ydx(X + {i ,X})<5(X-r) (6.3.7) 

The operator ordering is prescribed in momentum space with the exponential (Weyl) ordering 
and (6.3.7) in the noncommuting theory becomes 

j(t,k) = J dre lkr j(i,r) = J dx - i[A , X],) . (6.3.8) 

Note that the commutator does not contribute to current conservation because it is separately 
transverse. 

ydxel kX [io,k-X^ = J dxA [k-X,ef- x ]^ = (6.3.9) 

Therefore the 3-current is conserved as before. Its dual, e IMUa j ct satisfies the Bianchi identity, so the 
Seiberg-Witten mapping reads 

/dre-d /dse- 

J dre^ r F 0l = ojij J dxe**(X* - i[A ,X%) (6-3.10) 
dx< kX F * • 

Formulas (6.3.6) and (6.3.10) may be verified by comparison with the explicit 0(9) Seiberg-Witten 
map (6.3.2), which for field strengths implies 

Fjj,v = Fjiv - e aP {F af ,F Pv - A a d p F^) . (6.3.11) 

Upon setting 9 a0 = 0, 9 ij = 9e ij and 

^{x'+evAj) = eik . x ^ + i9k . £ ij^. _ le^kik^e^AjAn^ (6.3.12) 
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it is recognized that (6.3.6) and (6.3.10) reproduce (6.3.11). 
(ii) Seiberg-Witten map in higher even dimensions 

In dimensions higher than three the correspondence between the Bianchi identity and the conserva- 
tion of particle current is lost. The derivation of the Seiberg-Witten map calls for higher conserved 
currents, whose duals are two-forms. 

The introduction of such currents can be motivated by starting again from the commutative 
particle density p as expressed in (1.1.13) and its inverse p _1 as expressed in (1.1.19). Their product 

1 = j ds*(X-r)det-^i (6.3.13) 

is independent of the fluid profile X(x) and constitues a topological invariant. The Jacobian de- 
terminant in the above can be expressed in terms of the square-root determinant (Pfafnan) of the 
antisymmetric matrix {X^, X }: 



~ 2»n! / dx K^-^)^,n,...^AX i \X^}---{X^,X^} = Po J dxS(X. - r)C B (X) (6.3.14) 

where, in analogy with the 2-dimensional case, we identified Pfaff(0) with 1/po- Removing all 
n Poisson brackets from the above recovers the full density p. The removal of a single Poisson 
bracket {X i ,X^}, then, produces a sort of residual density pij in the corresponding dimensions, 
which becomes a candidate for the Seiberg-Witten commutative field strength: 

Pij = 2"-i(w - 1)! / ( M X -'•)<,,.<,,,, ,.,.{.V-.A^}-..{.Y-.A^ } (6.3.15) 

The current dual to p^, in momentum space, 

.J j: -■>■•■■ ■ = 2W _ 1 ^° jTj J dxe^{X [jl ,X j2 } ■ ■ ■ {X j2n ~ 3 , X j2n ~ 2] } (6.3.16) 

(the indices are fully antisymmetrized) is gauge invariant and conserved, ensuring that p^- satisfies 
the Bianchi identity. 

The corresponding current in the noncommutative case can be written by turning products into 
♦-products and Poisson brackets into (— i times) ^-commutators. The ordering of the exponential 
and other factors above has to be fixed in a way which ensures that the obtained current is conserved. 
Various such orderings are possible. For definiteness, we pick the ordering corresponding to the 
choice made in: 




e^^.P] * • . . e *»n-ik-x ^ [ x ; 2n - 3)X i2n- 2 ]^ (6.3.17) 
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This corresponds to Weyl-ordering the exponential and distributing it in all possible ways between 
the different commutators. Note that the volume of the Sj-integration space reproduces the factor 
l/(n — 1)! present in (6.3.16). 

To express compactly the above and to facilitate the upcoming derivations, we introduce anti- 
symmetric tensor notation. We define the basis 1-tensors vj representing the derivative vector field 
dj, and corresponding one-forms dx 3 . We consider the fundamental 1-tensor X and the 1-form k. 

X = X j Vj , k = kjdx j (6.3.18) 
All tensor products will be understood as antisymmetric. 

v i v k = 2 ( v i A v k ~ v k A Vj) , etc. (6.3.19) 

This amounts to considering v 3 - and dx k as anticommuting quantities. Scalar products are given by 
the standard contraction . 



v . ■ dx k = 5{ (6.3.20) 

We also revert to operator notation, dispensing with *-products and writing Tr for po / da;. Finally, 
we simply write L n _ x \ for the (n — l)-dimensional Sj-integration. 



Lrm^^^-i^i'-'P') (6 ' 321) 

Overall, the current in (6.3.17) is written as the rank-(2n — 2) antisymmetric tensor J, 



J = tr/ e !Slk - x XX---e ,s "' lk - x XX, (6.3.22) 

J(n-l) 

and its conservation is expressed by the contraction k • J = 0. The contraction of k acts on each X 
in a graded fashion. Using cyclicity of trace and invariance under relabeling the Sj, this becomes 



k- J = {n- l)tr / e isik - x [k-X,X]e iS2k - x XX---e is - lk - x XX (6.3.23) 

J(n-l) 

Using the identity 

[e isk ' x ,X] = f d Sl e iSlk - x \ik-X,X]e i{s - Sl)k - x (6.3.24) 
Jo 

we can absorb the Si-integration in (6.3.23) and bring it to the form 

k ■ J = --tr / [e iS2k ' x , X]XX • • • e is - lk - x XX (6.3.25) 

2 J(n-2) 
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Finally, using once more the cyclicity of trace, we see that the above contraction vanishes. This 
proves that the tensor J is conserved and, as a consequence, its dual pj k satisfies the Bianchi identity. 
As in the 2-dimensional case, we put 

F jk {k) = p jk (k) - u jk 5{k) , (6.3.26) 

and recover the commuting Abelian field strength, which can, in turn, be expressed in terms of a 
(commutative) Abelian potential Aj. 

In the above manipulations we freely used cyclicity of trace. In general this is dangerous, since 
the commuted operators may not be trace class. Assuming, however, that X becomes asymptotically 
x for large distances, the presence of the exponentials in the integrand ensures that this operation 
is permissible. 

As mentioned previously, the fully symmetric ordering is not the only one that leads to an 
admissible pj k . As an example, in the lowest-dimensional nontrivial case d = 4 we can alter the 
ordering by splitting the commutator as 

1 r 1 

Ji k = -tre ik *[X\X k ] -> jf = -itr / dsf{s)e ls ^X 1 e i(l ^ s) ^X k . (6.3.27) 
2i Jo 

If f(s) = — f(l — s) the above will be antisymmetric in (j, k) and conserved, as can explicitly be 
verified. Further, if /(s) satisfies 

/ ds(2s-l)f(s) = l, (6.3.28) 
Jo 

then (6.3.27) will also have the correct commutative limit. We obtain an infinity of solutions 
depending on a function of one variable f(s). This arbitrariness reflects the fact that the Seiberg- 
Witten equations are not integrable and therefore the solution for = depends on the path in 
the #-space taken for integrating the equations. For d = 4 the parameter space is a plane and the 
path from a given 6 to = on the plane can be parametrized by a function of a single variable, 
just like Jj k . The various solutions are related through field redefinitions. 

(iii) Seiberg-Witten map in higher odd dimensions 

The situation in odd dimensions differs in that we need to specify separately the components of the 
conserved current in the commutative and noncommutative directions. For d = 2n + 1 the current 
is of rank 2n — 1 and it can be constructed by a procedure analogous to the even-dimensional case: 
We start from the expression for the total particle current j 74 (1.1.13) and (1.1.14) and introduce 
2n — 2 commutators, one less than the number which would fully saturate it to (c, v). The temporal 
components J -? 1 ---.? 2 ™- 2 can be expressed as a rank-(2n — 2) antisymmetric spatial tensor J°, while 
the spatial components Jioh—hn-z can De expressed as a rank-(2n — 1) antisymmetric tensor J. 
Their fully ordered expressions are 
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tr / e iSlk ' x XX • • • e is - lk - x XX, (6.3.29) 

n — 1 J( n -l) 



J = tr / e is ° k ' x DXe iSlk ' x XX ■ ■ ■ e iSn - lk ' x XX . (6.3.30) 

J(n) 

The above expressions can be unified by introducing a temporal component for the field X M , namely 
X° = t (which is obviously commutative), and extending the one-tensor X also to include X°v . 
Further, we can Fourier transform in time and define k = k^dx 11 to include also the frequency k$. 
Then the corresponding (space-time) (2n — l)-tensor J acquires the form 

J = / dt tr / e lSlk ' x DX e iS2k ' x XX . . . e iSnk ' x XX . (6.3.31) 

J J{n) 

X° is absent in XX and, since DX° = 1, only sq + s± appears in the temporal component of J; 
integrating over si reproduces the factor l/(n — 1) appearing in (6.3.29). 

The above current is obviously gauge invariant. We shall prove that it is also conserved, that 
is, it satisfies k • J = 0. The contraction is 



k • J = J dt tr |e isik ' x k • DX e iS2k - x XX ■ ■ ■ e is " k - x XX 

n 

- ^e iSlk - x L>XXX---e iSmk - x [k-X,X]e iSm+lk - x ---Xx| (6.3.32) 

m=2 

(with s n+ \ = 0). By formula (6.3.24) and a similar one for the covariant time derivative, the above 
can be rewritten as 



k • J = f dt tr f \De iSlk X e !S2k ' x XX • • • e is - lk ' x XX 

n-1 

^2 e isikX DXXX ■ ■ ■ [e iSmk - x ,X]XX---Xxj. (6.3.33) 



m=2 



Due to the cyclicity of trace, the sum above telescopes and only the first term of the m = 2 
commutator and the second term of the m = n — 1 commutator survive. Altogether we obtain 
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k-J = [dttrf (De lSlk ' x + DXX + XDX) e iS2k ' x XX ■ ■ ■ e is "- lk ' x XX 

J J(n-1) 

= I dt tr f D (e 1Slk ' x XX) • • • e !S "- lk ' x XX 

J J(n-l) 

= [dttrf — - — D (e iSlk ' x XX ■ ■ ■ e iSn ~ lk ' x XX) (6.3.34) 



'(n-1) 

dt 4-tr f -^—e iSlk - x XX ■ ■ ■ e iSn ~ lk - x XX 



dt ./( n _i) n - 1 



0. 



which proves the conservation of J. Its dual p^ v satisfies the (2n + l)-dimensional Bianchi identity 
and can be used to define the commutative Abelian field stength. 



i^-(k) = Pij (k) - cu^ik) (6.3.35) 
F 0l (k) = p 0l (k) (6.3.36) 

In the above we gave separate derivations of the Seiberg-Witten map for even and odd dimen- 
sions. The two can be unified by demonstrating that each case can be obtained as a dimensional 
reduction of the other in one more dimension. This is treated next section. 



(iv) Dimensional reduction 

It is quite straightforward to see that the even dimensional Seiberg-Witten map is obtained from 
the d = 2n + 1 map by dimensional reduction. We assume a time-independent configuration in 
which (j = 1, . . . , 2n) do not depend on t and A vanishes. In this case DX vanishes and so does 
J in (6.3.30); only the component J° in (6.3.29) survives, reproducing the 2n-dimensional solution. 

The reduction from a fully noncommutative d = 2n + 2 case to the d = 2n + l case is only slightly 
subtler. For concreteness, we shall take t = x° to be canonically conjugate to the last dimension, 
call it z = x 2n+1 , which will be reduced; that is, 

[t, z] = i6 Q (9 Q = 9°' 2n+1 ) , [t,x l ] = [z,3?] = (z = l,...,2n). (6.3.37) 

This can always be achieved with an orthogonal rotation of the x M . The reduced configuration 
consists of taking all fluid coordinates other than X 2n+l to be independent of x 2n+1 and, further, 
the gauge potential corresponding to z = x 2n+1 to vanish. Specifically, 

X i = X'fat), (6.3.38) 
X° = t, (6.3.39) 
X 2n+i = z + e Q A (-K,t). (6.3.40) 
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With this choice the corresponding field strengths become 



[X\X j ] 
[X\X°] 
[X\X 2n+l ] 



i6 ij + i6 ik 6 je F f 



(6.3.41) 
(6.3.42) 
(6.3.43) 



0, 

i9 (D X l - i[X\ Ao}) = i9 ij e 2n+lfi F 3 



with (/i, v = 0, . . . , 2n) the field strength of a noncommutative d — 2n + 1 theory. 

The corresponding d = 2n + 2 Seiberg-Witten map reduces to the d = 2n + 1 map. Indeed, 
the current J in (6.3.22), now, is a rank-2n antisymetric tensor. When all its indices are spatial 
(1, . . . ,2n) it becomes a fully saturated topological invariant, that is, a constant; this reproduces 
a constant po,2n+i- When one of its indices is and the rest are spatial it vanishes, leading to 
Pi,2n+i = 0. When one of its indices is 2n + 1 and the rest are spatial it reproduces expression 
(6.3.30). Finally, when two of its indices are 0, 2n + 1 and the rest are spatial it reproduces (6.3.29), 
recovering the full commuting (2n + l)-dimensional Abelian field strength. 

We stress that the above reductions are not the most general ones. Indeed, mere invariance of the 
fluid configuration with respect to translations in the extra dimension does not require the vanishing 
of the gauge field in the corresponding direction. This means that we could choose X° = t + H(pc, t) 
(instead of X° = t) in both d = 2n + 1 and d = 2n + 2. The corresponding reduced theory contains 
an extra Higgs scalar in the adjoint representation of the (noncomutative) U(l) gauge group. Our 
Seiberg-Witten map in this situation reproduces, with no extra effort, the space-time derivatives of 
a corresponding commuting 'Higgs' scalar. 

The above complete reduction scheme (2n + 2 — >• 2n + 1 — >• 2n — > . . . ) is reminiscent of the 
topological descent equations relevant to gauge anomalies. It is possible to consider the fluid analogs 
of noncommutative topological actions and the mapping of topologically nontrivial configurations 
[56], but we shall not consider these issues here. 

A final comment: We demonstrated that gauge theory on noncommutative spaces has an effective 
description as a classical fluid. Other nonclassical situations may also be describable in a fluid 
dynamical language, thus revealing their dynamics as pertaining to a spacial kind of fluid. In 
particular, quantum mechanical many-body states could be effectively described in this fashion. 
This will be explored in the final Section. 
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7 MISCELLANEOUS TOPICS 

7.1 Quantized fluid mechanics 

Eulerian fluid mechanics, even though it can be formulated as an independent dynamical system 
in its own right, can also be understood as a good description for the physics of a distribution of 
particles at large length scales. At these scales, large compared to the typical particle separations, 
a continuum approximation can be made and the underlying particulate nature is seldom apparent 
or needed. If this fluid mechanics is quantized, the result is a quantum field theory where the 
hydrodynamical variables of density and velocity become quantum operators. For such an analysis, 
it is mathematically irrelevant whether the fluid theory emerged as an approximation to the particle 
theory. 

For quantization, the Poisson brackets (1.2.36), (1.2.39) and (1.2.40) are replaced by i/H times 
a quantum commutator of the p and v operators. There is no ordering ambiguity in the quantal 
version of (1.2.40) because p commutes with uJij, according to (1.2.39) and (1.2.41). Ordering does 
have to be prescribed when the Hamiltonian (1.2.33) is promoted to an Hermitian operator. The 
kinetic term may be taken as ~ v 1 pv 1 or as |(pv 2 + v 2 p); the two coincide because the [p, v] commu- 
tator is a c-number. The Heisenberg equations of motion then imply that the current contributing 
to the continuity equation (1.1.16) is the ordered Hermitian quantity 

j = I(pv + vp), (7.1.1) 

which also satisfies the quantum commutator analogs of (1.2.37) and (1.2.38). In the Euler force 
equation (1.1.17) the term non-linear in v emerges upon commutation with the quantum Hamilto- 
nian as | (v k dk v l + dk v % v h ) . 

It is noteworthy that the nature of the particles which constitute the fluid, in particular their 
statistics, does not seem to be important. For this reason, attempts to relate superfluidity to 
quantized hydrodynamics have been criticized [57]. However, one can see that this is clearly not 
the whole story by the following argument. 

Consider the classical action for a (non-linear) Schrodinger field on a <i-dimensional space R d ; 
d — 2, 3 are the cases of particular interest to us. 

dtdr[ifiip* tp - — V^j* ■ Vtp - V(ip*tp)} (7.1.2) 

With the Madelung Ansatz (G.4) (m = 1), the action (7.1.2) is brought to a form appropriate for 
an irrotational fluid. 

/ = J dtdr[6p - l -p{V9f - V(p)}, (7.1.3) 
where the classical potential V acquires an addition. 

V(p) = V(p)^(V^pf (7.1.4) 
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We see that the Schrodinger field theory (7.1.2) and fluid mechanics are equivalent with a reinterpre- 
tation of certain quantities, the particular nature of the fluid (such as its equation of state) will be 
characterized by the function V(p). The Schrodinger theory (7.1.2) can be quantized as a fermion 
field or a boson field. Therefore it must be possible to quantize the fluid theory (which is after all 
the same thing, being a particular parametrization or choice of coordinates on the classical phase 
space of the theory) as either a fermion theory or a boson theory. Since the symplectic structure, 
given by the #p-term, is what is relevant to this consideration, the specific nature of V(p) should 
not matter; it should be possible to quantize any fluid, irrespective of its dynamics, with either 
statistics. 

This observation is certainly not new. The fluid theory is described in terms of the density 
p and the current j = pV0. The algebra of these observables, the current algebra, is given by 
(1.2.36)-(1.2.38). This algebra is the same whether ip, ip* obey bosonic or fermionic commutation 
rules. Various observables of interest, such as the Hamiltonian, momentum density, etc., can be 
constructed in terms of p, j. Therefore, one can take the algebra of p, j as the starting point 
and construct the quantum theory in terms of a unitary irreducible representation of this algebra of 
operators. It is known that this algebra has many inequivalent realizations, allowing for the freedom 
of choosing different statistics, bosons and fermions corresponding to different representations [58]. 
The fluid theory can indeed carry information about the statistics of the particles which compose 
it. 

More generally, even for a one-component fluid, there are additional variables needed. Consid- 
ering the case of three spatial dimensions as an example, the free action is 



The fluid velocity is given in the Clebsch form v = V# + «V/9. The Clebsch parameterization of 
the velocities can be expressed as [see Section 7.3 (i), below] 



where the group element g is an element of SU(2). 



Using (7.1.7) with (7.1.6) gives the Clebsch parametrization for v with a = cos 7. Even though we 
have the additional variables a,/3 , the current algebra is the same as before. Evidently the same 
current algebra is realized by fluids that move with or without vortices. 

Generally speaking, in two dimensions, the existence of inequivalent representations is related to 
the nontrivial connectivity of the space of fields. If the phase space is simply connected, it cannot 
support double-valued (or many-valued) wave functions and one does not have the possibility of 




(7.1.5) 




(7.1.6) 
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different statistics. The topology of the phase space is therefore relevant to the question we are 
considering and one can try to characterize the inequivalent representations for the fluid current 
algebra along these lines. To see how this arises, it is instructive to consider the configuration space 
for identical particles. For N particles in two dimensions, a point in the configuration space is 
given by (xi,X2, ...,xjv), where x« e R 2 . The identity of particles tells us that we must make the 
identification 

(xi,x 2 , ...,Xj, ...,Xj, ...,xtv) ~ (xi,x 2 , ...,Xj, ...,Xj, ...,xjv) (7.1.8) 

Further, we must impose the condition that the locations of the particles do not coincide, i.e., 
Xj 7^ Xj if % 7^ j. The resulting configuration space Cn has nontrivial connectivity; in fact, the first 
homotopy group is given by IIi(Cjv) = -Bjvj where is the braid group. The fact that IIi(Cjv) ^ 
shows that the configuration space can support many-valued wave functions. Under exchange of 
particles, the wave functions get a phase; the various phase factors for different exchanges form a 
unitary irreducible representation of IIi(Cjv)- The phase for a single exchange can be arbitrary and 
so we get the possibility of arbitrary statistics. (An analogous construction can be done in three 
dimensions; the first homotopy group is then the permutation group Sn and we get the possibility 
of fermions or bosons.) 

Returning to the case of two dimensions, X = (x,y), notice that the existence of different 
representations is also related to the existence of a flat potential, which distinguishes between the 
representations. For the iV-particle Heisenberg algebra, we can write a representation in terms of 
complex variables w = x — iy,w = x + iy. 

X; = (Wi,Wi) 

i 



Pw, — Pwi + 9 , ^2 



2k Wi — w.j 



d . . d 



dwi ' Wt dwi 

k is a parameter specifying the representation, p gives the standard Schrodinger representation, and 
the summed expressions comprise the potentials, which make the representation (7.1.9) different 
from the standard one. It is easily verified that the quantities in (7.1.9) obey the Heisenberg algebra. 
[The commutator of the momenta leads to <5(xj — Xj), but since we do not allow coincidence of 
particle locations (or the wave functions vanish at coincidences), this is zero.] We can write the 
representation (7.1.9) as 

± = U- l xU, P = U' 1 pU (7.1.10) 

where U is given by 

U = exp ^[lnK - Wj) - \n{wi - %)]J (7.1.11) 
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The representation (7.1.9) is however not equivalent to the Schrodinger representation because U 
is not single-valued and so does not give a unitary transformation. Notice that the addition to the 
expression for the momentum operators is a "flat potential", in the sense that its field strength, 
which is given by the commutators [P Wi ,P 0j ], is zero (with the condition of removing coincidences 
of particle locations). This flat potential can be written as t/ -1 [j3, U], but it should be kept in mind 
that U is not a genuine unitary transformation. 

Similar results hold in general. Let {4> a } be a set of observables, which may include the identity, 
obeying a commutation algebra of the form 

[0 a ,0 b ]=C a c b c . (7.1.12) 

Let (pi and <p^ be two representations. We write (f)^ = 4>a + A a - The fact that 0„ and (fffi 
obey the same algebra shows that A a is a flat potential, i.e., 

A b ] - [tf \ Aa] + [A a , A b ) - Cl b A c = 0. (7.1.13) 

If A a = tf[$\U\ for some unitary transformation U, the two representations are unitarily 
equivalent. However, if A a obeys the zero-curvature condition (7.1.13) but cannot be written as 
A a = W[(j)a 1 \U] for some unitary U, the two representations are inequivalent. 

Clearly, the existence of inequivalent representations and the existence of a flat potential which 
cannot be written as A a = U^[(f)a\U\ are related to the topology of the space of fields, or the 
phase space, if we are thinking in terms of canonical quantization. The IPs which we use form a 
representation of some nontrivial ITi(C); flat potentials require in general that the first cohomology 
group of the phase space should be nontrivial. We have seen this explicitly at the level of particles. 
When we generalize to a field theory one can ask whether statistics can be obtained in such terms. In 
fact, a description of the boson and fermion representations and an expression for the corresponding 
flat potential have been obtained [58]. It is thus possible to identify an operator, namely the flat 
potential A a , which can distinguish the statistics of the underlying particles. Nevertheless, the 
situation is not entirely satisfactory. The flat potential has been obtained only for a subspace with 
fixed value of the particle number N and it is in terms of the phase of the fermionic ground state 
wave function for N particles. In our opinion, this answer partially begs the question. Also, the 
connection to the topology of the phase space in unclear. A more direct method, in terms of the 
fields, would certainly be better. 

Turning to the field variables p, j, we first note that, already, the use of a group element 
to parameterize the velocities as in (7.1.7) entails making certain additional assumptions about 
the topology of the phase space. One has to determine on a physical grounds whether such a 
parameterization is justified and if so, which group should be used. For the parameterization 
(7.1.7), the phase space is given by 

V = {set of all mapsp : R 3 -> R+, g : R 3 -> G = SU(2)} (7.1.14) 
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The topology of the phase space is specific to SU (2): there is a compact U(l) direction corresponding 
to 9. Since p and 9 are canonically conjugate, the operator U = exp(i27i j dr p) shifts 9 by 2ir. 
The compactness of 9 requires that all observables be invariant under the action of U ; equivalently, 
the spectrum of J dr p must be an integer spaced. Thus the compactness of U(l) leads to an 
underlying particle description, with the integer portion of j drp being the number of particles 
N. Therefore, for hydrodynamics with an underlying particle structure, the use of a compact 9- 
variable is appropriate. Notice also that J dr p is a Casimir operator for the algebra so that one 
can also restrict to representations with a fixed value of N. If we have only (p,9), this argument 
just recaptures the Schrodinger field description. 

One can now try to see how the flat potentials mentioned above can arise, in the field description. 
So far, there is no direct construction of the potential A a in terms of p and j. In our formulation 
of the Clebsch parametrization, the phase space given by (7.1.8) does not have nontrivial first 
homotopy group. One might argue that this is because we have to impose some further condition, 
analogous to the condition of excluding coincidence of locations at the particle level. It may be 
that a similar condition must be imposed, expressed in fluid variable terms, on the space of (7.1.14) 
to obtain the necessary structure for a quantum hydrodynamics, which takes account of particle 
statistics. It is not yet clear what such a condition would be. Also, representations of the algebra 
of observables are characterized, in purely algebraic terms, by the values of the Casimir operators. 
Thus, alternatively, one can ask the question: What are the Casimirs in terms of the fluid observables 
which discriminate between the different statistics? How are these Casimirs related to the topology 
of the phase space? As yet there is no completely satisfactory answer to these questions. 

7.2 Fluids in intense magnetic fields 

Another interesting dynamical system is a charged fluid in an intense magnetic field. This exhibits 
a noncommutativity, which is the fluid mechanical version of noncommuting coordinates for a point 
particle in a magnetic field, so intense that it effects reduction to the lowest Landau level. 

(i) Particle noncommutativity in the lowest Landau level 

Before describing the motion of a charged fluid in an intense magnetic field, we review the story 
for point particles on a plane, with an external and constant magnetic field B perpendicular to the 
plane [59]. The equation of motion for the 2- vector r = (x, y) is 

mtf = -e ij v j B + f(r), (7.2.1) 
c 

where v is the velocity r, and f represents other forces, which we take to be derived from a potential 
V: f = — VV^. The limit of large B is equivalent to small m. Setting the mass to zero in (7.2.1) 
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leaves a first order equation. 

r* = (7-2.2) 
This may be obtained by taking Poisson brackets of r with the Hamiltonian 

H = V, (7.2.3) 

provided the fundamental brackets describe noncommuting coordinates, 

0V} = ^". ( 7 - 2 - 4 ) 

so that 

r l = {H ,r 1 } = {r^r'jdjV = ^'T(r). (7.2.5) 

The noncommutative algebra (7.2.4) and the associated dynamics can be derived in the following 
manner. The Lagrangian for the equation of motion (7.2.1) is 

L = imv 2 + -v • A - V. (7.2.6) 

When we choose the gauge A = (0, Bx) and set m to zero, (7.2.6) leaves 

eB 

L = —xy-V(x,y) J (7.2.7) 

which is of the form pq — h(p,q), and one sees that (^-x,y) form a canonical pair. This implies 
(7.2.4), and identifies V as the Hamiltonian. 

Additionally, we give a canonical derivation of noncommutativity in the m — > limit, starting 
with the Hamiltonian 

H = h V. (7.2.8) 

2m 

H gives (7.2.1) upon bracketing with r and v, provided the following brackets hold. 

{r\r j } = (7.2.9) 
{ 7l \r i } = 5 lj (7.2.10) 

{7r>n = -— e ij (7.2.11) 

c 

Here 7r is the kinematical (non-canonical) momentum, mr, related to the canonical momentum p 
by 7r = p — -A. 

We wish to set m to zero in (7.2.8). This can only be done provided w vanishes, and we impose 
77 = as a constraint. But according to (7.2.11), the bracket of the constraints {tt 1 ,^} = C % i = 
— — e l i is non-zero. Hence we must introduce Dirac brackets: 

c 

{O l5 2 } D = {Ox, 2 } - {Oi, n k }{C- l ) kl {n\ 2 }. (7.2.12) 
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With (7.2.12), any Dirac bracket involving 7r vanishes, so w may indeed be set to zero. But the 
Dirac bracket of two coordinates is now non-vanishing. 

{r\ ri} D = -{r\ 7r fc }^e fe V, = -^e« (7.2.13) 

en eB 

In this approach, noncommuting coordinates arise as Dirac brackets in a system constrained to lie 
in the lowest Landau level. 

A quantum mechanical perspective on this result is the following [60] . Let us label the degenerate 
quantum Landau states by (N,n), where N labels the Landau level and n the degeneracy of that 
level. Consider the computation of the matrix element in the lowest (N = 0) Landau level of the 
commutator of x as y. 

< 0, n\[x, y]\0, h >=< 0,n\xy\0,h > — < 0,n \yx\0,h > 

= J^(< 0, n\x\Nn >< Nn'\y\0, h > - < 0, n\y\Nn X Nn'\x\0, h >) (7.2.14) 

Nn' 

When the intermediate state sum is carried over all Landau levels, the two terms in the sum cancel 
each other, and we find a vanishing result; coordinates commute. Suppose however, we truncate the 
sum at the lowest level. A simple calculation gives a non- vanishing result, consistent with (7.2.4). 
This shows that noncommuting coordinates arise from a truncation of the Hilbert space. 

[An amusing generalization determines the coordinate commutator when the first N Landau 
levels are retained, both in the external states and the intermediate state sum. The result is that 
the only non- vanishing matrix element of the [x, y] commutator is in the highest Landau level. 

i fie 

< N, n\[x, y]\N, ri >= (N + 1)6 (n, n') (7.2.15) 

eB 

Here 5(n, n') is either the discrete or continuous delta function, depending whether the Landau 
degeneracy is exhibited in a discrete or continuous manner. The result reduces to the previous, 
when only the lowest Landau level is kept, and also shows that as more and more states are 
included, the non commutativity is pushed into ever higher states.] 

(ii) Field noncommutativity in the lowest Landau level 

We now turn to the equations of a charged fluid with density p and mass parameter m (introduced for 
dimensional reasons) moving on a plane with velocity v in an external magnetic field perpendicular 
to the plane, p and v are functions of t and r and give an Eulerian description of the fluid. The 
equations that are satisfied are the continuity equation (1.2.42), and the Euler equation (1.2.43), 
which acquires an additional term, beyond the force team, due to the magnetic interaction [compare 
(5.4.19) reduced to the Abelian case]. 

m^ + mv-Vu^-eVB + f (7.2.16) 
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Here f describes additional forces, e.g. — ^VP where P is pressure. We shall take the force to be 
derived from a potential in the form 

f(r) = -V-j- / d 2 rV. (7.2.17) 



6p(r) 

[For isentropic systems, the pressure is only a function of p. Here we allow more general dependence 
of V on p (e.g. nonlocality or dependence on derivatives of p).\ 

The relevant equations follow by bracketing p and v with the Hamiltonian 



provided that fundamental brackets are taken as 

{p(r),p(r')} = (7.2.19) 
{7r(r),p(r')} = V5(r-r') (7.2.20) 

{7r 4 (r), ttV)} = (rnuj{?) + ^f) 6(t - r') (7.2.21) 

where e u u;(r) is the vorticity — djv 1 , and 7r = mv. These are just the rescaled (by m) versions 
of (1.2.36), (1.2.39) and (1.2.40), except that the last is modified to include the constant magnetic 
field. 

We now consider a strong magnetic field and take the limit m — > 0, which is equivalent to large 
B. Equations (7.2.16) and (7.2.17) reduce to 

v* = — £_ e « A * f dW; (7.2.22) 
eB Or? Sp(r) J 

Combining this with the continuity equation gives the equation for the density "in the lowest Landau 
level." 

P(r) = —Pir)^^ [ d\V (7.2.23) 



eB dr % dri 8p(v 

(For the right hand side not to vanish, V must not be solely a function of p.) 

The equation of motion (7.2.23) can be obtained by bracketing with the Hamiltonian 



H = j d 2 rV, (7.2.24) 
provided the charge density bracket is non- vanishing, showing non-commutativity of the p's [61]. 

{p(r), p(r')} = ~<#aip(T)d s 6(T - r') (7.2.25) 
er> 

Hq and this bracket may be obtained from (7.2.18) and (7.2.19) - (7.2.21) with the same Dirac 
procedure presented for the particle case: We wish to set m to zero in (7.2.18); this is possible only 
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if 7r is constrained to vanish. But the bracket of the 7r's is non-vanishing, even at m = 0, because 
B 7^ 0. Thus at m = we posit the Dirac brackets 

{0 1 (r 1 ),0 2 (r 2 )} D = 

{0 1 (r 1 ),0 2 (r 2 )} - | cZMdVaOi^O.Tr^r;)}^- 1 )^^,^)!^^)^,^)}, 

(7.2.26) 

where C u is the bracket in (7.2.21), (at m = 0) so that 

(C^^n, r 2 ) = -L^n)^ - r 2 ). (7.2.27) 

Hence Dirac brackets with 7r vanish, and the Dirac bracket of densities is non-vanishing as in 
(7.2.25). 



{p(r),p(r>)} D = 

~ J d 2 r"{p(r),n\r"}p(r")e^{^(r"),p(r')} 



c 



-— e^p(r)«9^(r - r') (7.2.28) 

The p-bracket enjoys a more appealing expression in momentum space. Upon defining 

p{p) = J d 2 re ip - r p(v) (7.2.29) 

we find 

{p(p),p(q)} D = ~^WP(P + q)- (7-2.30) 

The brackets (7.2.25), (7.2.30) give the algebra of area preserving diffeomorphisms. Indeed (7.2.30) 
follows from the full diffeomorphism algebra (1.2.38) with the identification that p in (7.2.30) is 
related to f in (1.2.38) by p = ' J (),j J . 

The form of the charge density bracket (7.2.25), (7.2.28), (7.2.30) can be understood by reference 
to the particle substructure for the fluid. Take as in (1.1.9) 

p(r) = ^5(r-r„), (7.2.31) 

n 

where n labels the individual particles. The coordinates of each particle satisfy the non-vanishing 
bracket (7.2.4). Then the {p(r),p(r')} bracket takes the form describe above. 

(iii) Quantization 

Quantization before the reduction to the lowest Landau level is straightforward. For the particle 
case (7.2.9) - (7.2.11) and for the fluid case (7.2.19) - (7.2.21) we replace brackets with i/h times 
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commutators. After reduction to the lowest Landau level we do the same for the particle case 
thereby arriving at the "Peierls substitution," which states that the effect of an impurity [V in 
(7.2.6)] on the lowest Landau energy level can be evaluated to lowest order by viewing the (x,y) 
arguments of V as non-commuting variables. 

However, for the fluid case quantization presents a choice. On the one hand, we can simply 
promote the brackets (7.2.25), (7.2.28), (7.2.30) to a commutator by multiplying by i/h. 

[p(r),p(r')] = ih^dip^djSiT - r') (7.2.32) 

[p(p), p(q)] = ih^Wp(p + q) (7-2.33) 

Alternatively we can adopt the expression (7.2.31), for the operator p(r), where the r n now 
satisfy the non-commutative algebra, 

[ri,ri,] =-ih^6 m , (7.2.34) 

and calculate the p commutator as a derived quantity. 

However, once r n is a non-commuting operator, functions of r n , even S— functions, have to be 
ordered. We choose the Weyl ordering, which is equivalent to defining the Fourier transform as 

p(p) = J>*" rB - (7.2.35) 

n 

With the help of (7.2.34) and the Baker-Hausdorff lemma, we arrive at the "trigonometric algebra". 
[62] 

[p(p), p(q)] = 2z sin (J^e^q^J p(p + q) (7.2.36) 

This reduces to (7.2.33) for small ft. 

This form for the commutator, (7.2.36), is connected to a Moyal * product in the following 
fashion. For an arbitrary c- number function /(r) define 

< / >= J d 2 rp(r)f(r) = JL J d 2 pp(p)f(-p). (7.2.37) 

Multiplying (7.2.36) by /(— p)g(— q) and integrating gives 

[</>,< g >]=<h>, (7.2.38) 

with 

h(r) = (f*g)(r)- (g*f) (r) (7.2.39) 

where the * product is defined in (6.1.5). Note however that only the commutator is mapped into 
the * commutator. The product < / >< g > is not equal to < / * g >. 

The lack of consilience between (7.2.33) and (7.2.36) is an instance of the Groenwald-VanHove 
theorem which establishes the impossibility of taking over into quantum mechanics all classical 
brackets. Equations (7.2.36) - (7.2.39) explicitly exhibit the physical occurrence of the * product 
for fields in a strong magnetic background. 
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7.3 Fluids and ferromagnets 

A quantum system that has some similarity in its description to fluid mechanics is the ferromagnet 
[63]. The requirements that we put on the ferromagnetic theory is that it give rise to a local spin 
algebra, [compare (5.3.3)] 

[S a (r),S b (r')} = ihe abc S c (r)5 3 (r-r') (7.3.1) 

and that some Hamiltonian H encode the specific nature of the ferromagnet. Thus we expect the 
action to be 

I = I - fdtH. (7.3.2) 



Here Iq involves the canonical 1-form that leads to (7.3.1), where the S a are definite functions of 
the canonical variables. 

(i) Spin algebra 

In our treatment of non Abelian fluids we obtained the Poisson bracket version of (7.3.1), i.e. (5.3.3), 
by starting from (D.l) specialized to SU(2): geSU(2),K = a 3 /2i, 

S a T a = gKg-\ (7.3.3) 

However, for the ferromagnet, there is no place for a density variable p; rather we should set p in 
(D.l) to a constant and posit for the canonical portion of the Lagrange density 

C = -ttra 3 g~ 1 g. (7.3.4) 

The canonical variables are three in number, corresponding to the three parameters specifying the 
SU (2) group element. Consequently the 2-form involves a singular matrix (3x3 and anti-symmetric) 
with no inverse. Nevertheless, one can overcome this obstacle and conclude that (7.3.1) continues 
to hold with S a defined by (7.3.2). 

This is achieved with the following steps. With the notation of Sidebar D, we determine that 
the Co reads [compare (D.2 and (D.3)] 

C = <pC a b S b . (7.3.5) 

The canonical 1— form has components a a = C a b S b [compare (D.10)] and the 2-form reads 

5 _5_ 

Q>b r- h O'a — Jab) 



5(p a 5(p b 

f ab (r, r') = -5(r - r')C c a C d b S e e cde . (7.3.6) 

Consider now the left translation of g: 5g = e a a a g/2i, or equivalently 8tp a = —e b c b a , which leads 
to 5S a = e a h c e b S c . The generator of this transformation satisfies, according to (A. 2), 

|^ = -efc f b C c b C d a S e e cde = -e c C d a S e e cde . (7.3.7) 
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This equation is solved by 

G = - I drS a e a . (7.3.8) 



From the general theory, (A. 18), we knew that Poisson bracketing with G generates the above 
transformation on any function of the phase space variables, in particular on S b . 

e abc e h (r')S c (r') = {G, S b (r')} = { J dr e a (r) S a (r), S b (r')} (7.3.9) 

Stripping away the arbitrary parameter e a (r) leads to a Poisson bracket, which gives (7.3.1) upon 
quantization. 

It is instructive to obtain this result by applying the procedures in Sidebar A(b), relevant to 
singular 2-forms. Note that S a satisfies 

JLs a = c a m f mb . (7.3.10) 

Consequently it follows from (A. 10) that the admissibility criterion (A. 13) is obeyed. The bracket 
of local spins reads 

{S a (r),S b (r')} = I dr"dr>" S a (r) /°V, r")—?—S\r'). (7.3.11) 
J Oip c (r') dip d (r") 

It follows from (7.3.5), (A. 10), (A. 11) and (7.3.9) that (7.3.10) reduces to the Poisson bracket version 
of (7.3.1). 

(ii) Momentum density algebra 

The above shows that the ferromagnet in the continuum approximation, may be considered, as far 
as the canonical structure is concerned, as fluid mechanics with a constant density p. It has been 
known for some time that there is difficulty in defining a momentum density for the ferromagnet [64]. 
The momentum density must generate the coordinate transformation 5r = —d(r),5(p a = S ■ V<y9 a , 



,9a 



S-Vpa (7.3.12) 



and the generator G should solve [see (A. 17), (D.12), (D.13b)] 
SG 



J drd-V<p b C c b C d a S e e cde 
— f drS-trVgg^S + i f drV ■ Str-^- g' 1 S (7.3.13) 

8<A> J J Olf a 



The second equality follows from the first with (D.2) and (7.3.3). The last term prevents equation 
(7.3. 12) from being integrable. Therefore, a generator of local translations, i. e. a momentum density, 
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cannot be denned, except in the case that V ■ 8 = 0. Translations with transverse 8 correspond to 
volume-preserving transformations; and are canonically implemented with the generator 

G(S) = J drd -trV g g' 1 S. (7.3.14) 

They include total momentum and orbital angular momentum. 

It is instructive to see why the procedures in Section A(b) fail to a produce proper generator, 
which by Noether's theorem for the Lagrange density (7.3.4) is 

V = a a (i P )Vi Pa = 2trVgg- 1 S. (7.3.15) 

The bracket of the quantities 

{r(r),n?} = / dr '^'(^y^w) / Q Vy0(^^(?)) (7.3.16a) 
is evaluated with the help of the projected inverse to f ab , which satisfies (A. 11). The result is 

{^(r),^!-')} = V i (r')d i 8(r-r')+V j (r)d j 8(r-r') (7.3.16b) 
- A J (r') 8i S(r - r') - A^rjc^r - r), 

where 

A = a m P a m V^ a , (7.3.17) 

P™ being the projector on the zero-modes of f a b- The usual momentum density algebra should not 
contain the A modification(7.3.16), (7.3.17), [compare (1.2.38)]. This modification also prevents 
the quantities J drS (r)P(r) from realizing the full diffeomorphism algebra (1.2.62)-(1.2.64), except 
in the restricted case of volume preserving diffeomorphisms V ■ 8 = 0. Finally we remark that with 
non vanishing A the acceptability condition (A. 13) is not met, and the bracket (7.3.17) does not 
satisfy the Jacobi identity. 

A solution to all these problems with implementing arbitrary diffeomorphisms is found through 
the connection to fluid mechanics, where the ferromagnetic phase space is enlarged to include p as 
in (Dl) [63]. Then the 2-form is non-singular, the symplectic structure and conventional Poisson 
brackets exist [see (D.5)]. Either the equation (A. 17) or Noether's theorem give the generator 
arbitrary diffeomorphism, viz. as the momentum density, as 

V = 2ptrV gg" 1 S, (7.3.18) 

which satisfies the proper algebra (1.2.38), or (7.3.16c) without the A spoiler. The p-extended 
phase space can be reduced, in the Dirac sense, by imposing the first class constraint p=constant. 
This will achieve the ferromagnet's phase space. 
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7.4 Non-Abelian Clebsch parameterization 

(or, non-Abelian Chern-Simons term as a surface integral — a holographic presenta- 
tion) 

In Sidebar B we described the Clebsch parameterization (B.l) for a 3- dimensional vector potential, 
which casts the Abelian Chern-Simons density (B.ll) into a total derivative form (B.ll), so that 
its 3-dimensional volume integral obtains its value from the 2-dimensional surface that bounds the 
integration volume. One may pose an analogous question about a non-Abelian gauge potential: 
how should it be parameterized so that the non-Abelian Chern-Simons term 

CS{A) = e^\A« dj A) + 1 f abc At A) A c k ) = -£^r(~ A d j A k + ~A A, A k ) (7.4.1) 

becomes a total derivative, and its integral becomes a surface term? (In the second equality we 
use a matrix representation of the Lie algebra: A = A a T a ,trT a T b = —5 ub /2.) This question has 
been answered in different ways by physicists and by mathematicians. We shall here present both 
constructions, but first we explain the reason for the difference. 
Consider CS(A) as a 3-form on an arbitrary manifold 

CS(A) = A a dA a + l -f abc A a A b A c 

A a = A a dx i (7.4.2) 



d CS{A) =F a F a F a = - F% dx l dx j (7.4.3) 



(In our form notation, we suppress the wedge product.) It follows that dCS(A) is a 4-form. 

1 
2 

If the Chern-Simons term is a total derivative, CS(A) = dO, the dCS(A) = 0. So the possibility 
of expressing the Chern-Simons as a total derivative requires that F a F a vanish. The "physics" 
approach to the problem achieves vanishing of the 4-form dCS(A) = F a F a by working on a 3- 
dimensional manifold, which does not support 4-forms. [The Abelian Clebsch parameterization 
(B.l) is given in 3-space!] The "mathematics" approach remains with the 4-dimensional space, 
viewed as a Kahler manifold, but requires that certain components of F vanish, see below. 

(i) Total-derivative form for the Chern-Simons term in 3-space 

In Sidebar B, eqs. (B.5) - (B.10), we gave an analytic/geometric construction of the Abelian Clebsch 
parameterization for a vector potential. However, there is another approach to this Abelian problem, 
which contains clues for the non-Abelian generalization. So we re-analyze the Abelian case. 

The alternate method for constructing the Clebsch parameterization for an Abelian potential 
relies on projecting the potential from a non-Abelian one, specifically, one for SU(2). We consider 
an SU(2) group element g and a pure gauge SU(2) gauge field, whose matrix- valued 1-form is 

g- 1 dg = V a ^. (7.4.4) 
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where a a are Pauli matrices. It is known that 



tr(<T x dg) 3 



-\e abc V a V b V 



-fyW 3 



(7.4.5) 



is a total derivative; indeed its spatial integral measures the winding number of the gauge function g 
[66]. Since V a is a pure gauge, we have 



so that if we define an Abelian gauge potential A by selecting one SU(2) component of (7.4.4) (say 
the third) A = V 3 , the Abelian Chern-Simons density for A is a total derivative, as is seen from 
the chain of equations that relies on (7.4.5) and (7.4.6), 



and concludes with an expression known to be a total derivative. Of course A = V 3 is not an 
Abelian pure gauge. 

Note that g depends on three arbitrary functions, the three SU(2) local gauge functions. Hence 
V 3 enjoys sufficient generality to represent the 3-dimensional vector A. Moreover, since A's Abelian 
Chern-Simons density is given by tr(<? -1 dg) 3 , which is a total derivative, a Clebsch parameterization 
for A is easily constructed. We also observe that when the SU(2) group element g has nonvanishing 
winding number, the resultant Abelian vector possesses a nonvanishing Chern-Simons integral, that 
is, nonzero magnetic helicity. Specifically, the example of the Clebsch-parameterized gauge potential 
in (B.15), (B.18) is gotten by projecting onto the third direction of a pure gauge SU(2) potential 
constructed from the group element g = exp(a a /2i)f a a(r)2A = itr a 3 g~ x dg. An even more direct 

3 2 3 

example is given by the SU(2) group element g = e ^ 13 e^ 7 e^ e . Then A = d6 + cos'jdfi. 

Now we turn to the non Abelian problem, which we formulate in the following way : For a given 
group H, how can one construct a potential A a such that the non- Abelian Chern-Simons integrand 
CS(A) is a total derivative? 

In the solution that we present [67], the "total derivative" form for the Chern-Simons density 
of A a is achieved in two steps. The parameterization, which we find, directly leads to an Abelian 
form of the Chern-Simons density, 



for some 7. Then Darboux's theorem (or usual fluid dynamical theory) ensures that 7 can be 
presented in Clebsch form, so that 7 d7 is explicitly a total derivative. 

We begin with a pure gauge g~ x dg in some non- Abelian group G (called the Ur-group) whose 
Chern-Simons integral coincides with the winding number of g. 



dV a = -\e ahc V h V\ 



(7.4.6) 



AdA = V 3 dV 3 = -V W 3 = ijtrfo- 1 d^) 3 



(7.4.7) 



A a dA a +\f ahc A a A h A c = 1 d ll 



(7.4.8) 



w{9) = Teb / d3rCS (^ ld ^) 



24vr 2 



1 



/ 



tr{g~ l dgf 



(7.4.9) 
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We consider a normal subgroup H C G, with generators T a , and construct a non-Abelian gauge 
potential for H by projection. 

A a oc t^TV 1 ^) (7.4.10) 

Within if, this is not a pure gauge. We determine the group structure that ensures the Chern- 
Simons 3-form of A a to be proportional to tr(g _1 dg) 3 . Consequently, the constructed non-Abelian 
gauge fields, belonging to the group H, carry quantized Chern-Simons number. Moreover, we 
describe the properties of the Ur-group G that guarantee that the projected potential A a enjoys 
sufficient generality to represent an arbitrary potential in H. 

Since tr(g _1 dg) 3 is a total derivative for an arbitrary group [66] (although this fact cannot 
in general be expressed in finite terms) our construction ensures that the form of A a , which is 
achieved through the projection (7.4.10), produces a "total derivative" expression (in the limited 
sense indicated above) for its Chern-Simons density. 

Conditions on the Ur-group G, which we take to be compact and semi-simple, are the following. 
First of all G has to be so chosen that it has sufficient number of parameters to make tr(T a g~ l dg) 
a generic potential for H. Since we are in three dimensions, an if-potential A\ has 3 x d\m.H 
independent functions; so a minimal requirement will be 

dim G > 3 dim H . (7.4.11) 

Secondly we require that the if-Chern-Simons form for A a should coincide with that of g~ l dg. As 
we shall show in a moment, this is achieved if G/H is a symmetric space. In this case, if we split 
the Lie algebra of G into the if-subalgebra spanned by T a , a — 1, ... , dim H, and the orthogonal 
complement spanned by S A , A = 1, . . . , (dim G — dim H), the commutation rules are of the form 

[T a ,T b } = f abc T c , (7.4.12a) 
[T a ,S A ] = h aAB S B , (7.4.12b) 
[S A , S B ] = N h aAB T a . (7.4.12c) 

{h a ) AB form a (possibly reducible) representation of the if -generators T a . The constant N de- 
pends on normalizations. More explicitly, if the structure constants for the Ur-group G are named 
fabc, { a , b, c} = 1, . . . , dimG, then the conditions (7.4.12a-c) require that f a b c vanishes whenever an 
odd number of indices belongs to the orthogonal complement labeled by A, B, ... Moreover, f a b c are 
taken to be the conventional structure constants for H and this may render them proportional to 
(rather than equal to) f a b c - 

We define the traces of the generators by 

tr(T a T b ) = -NiS ab , tr(S A S B ) = -N 2 5 AB 

ti{I a S A ) = . (7.4.13) 
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We can evaluate the quantity tr[S A , S B ]T a = trS A [S B ,T a ] using the commutation rules. This 
immediately gives the relation TViiV = iVjj. 

Expanding g~ l dg in terms of generators, we write 

g- 1 dg = (T a A a + S A a A ), (7.4.14) 

which defines the if -potential A a . Equivalently 

A a = -+rtT{T a g- 1 dg). (7.4.15) 

iVi 

From d(g _1 dg) = —g^ 1 dg g^ 1 dg, we get the Maurer-Cartan relations 

F a = dA a +y abc A b A c = -—h aAB a A a B , 
2 2 

da A +h aBA A a a B = . (7.4.16) 



Using these results, the following chain of equations shows that the Chern-Simons 3-form for the 
if -gauge group is proportional to tr(g~ l dg) 3 . 

; {A a dA a +\f abc A a A b A c ) = -^—{A a dA a +2 A a F a ) 



16tt 2V 3J ' 48vr 2 



48vr 



(A a dA a -Nh aAB A a a A a B ) 



2 



48tt 



{A a dA a +Nda A a A ) 



2 



. — trAdA-\ trdaa 

48tt 2 \Ni N 2 

— tr(A dA +a da) 

487r 2 iVi V ; 

tr(g^dg) 3 (7.4.17) 



48yr 2 iVi 

In the above sequence of manipulations, we have used the Maurer-Cartan relations (7.4.16), which 
rely on the symmetric space structure of (7.4.12a-c), and the trace relations (7.4.13), along with 

JViJV = N 2 

We thus see that J CS(A) is indeed the winding number of the configuration g G G. Since 
tr(g~ l dg) 3 is a total derivative locally on G, the potential (7.4.15), with the symmetric space 
structure of (7.4.12a-c), does indeed fulfill the requirement of making CS(A) a total derivative. It 
is therefore appropriate to call our construction (7.4.15) a "non-Abelian Clebsch parameterization". 

In explicit realizations, given a gauge group of interest H, we need to choose a group G such 
that the conditions (7.4.11), (7.4.12a-c) hold. In general this is not possible. However, one can 
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proceed recursively. Let us suppose that the desired result has been established for a group, which 
we call i?2- Then we form H C G obeying (7.4.12a-c) as H = Hi x if 2 , where H± is the gauge 
group of interest, satisfying dimG > 3 dimi^. For this choice of H, the result (7.4.17) becomes 

CS(lfi) + CS(H 2 ) = -JL^^ dg) 3 . (7.4.18) 

But since CS(H 2 ) is already known to be a total derivative, (7.4.18) shows the desired result: CS(Hi) 
is a total derivative. 

With SU(2) as the f/r-group and H = U(l) or 50(2), we achieve Clebsch parameterization for 
an Abelian potential, as explained in (7.4.4) - (7.4.7). [T a = fj; S A = f^}] 

For an explicit working out of a non Abelian case, we consider a potential for SU(2) ~ 0(3), 
which possesses nine independent functions. We take G = 0(5), H = 0(3) x 0(2). We consider the 
4-dimensional spinorial representation of 0(5). With the generators normalized by tr(T a T b ) = —5 ab , 
the Lie algebra generators of 0(5) are given by matrices of Pauli matrices. 



2% ( 


a a 




o a 


27 f 


-I 




/ 


1 


(° 


iV2 





(7.4.19) 




T" generate 0(3), with the conventional structure constants e a b c , and T° is the generator of 0(2). 
S, S are the coset generators. 

A general group element in 0(5) can be written in the form g = Mhk where h 6 0(3), k 6 0(2), 
and 

1 / I — ~(w x w) ■ a — w ■ a \ 

M = —^^^^=^^^= (7.4.20) 
wl+w-w-i(wx w) 2 y w • a I + |(w X w) • a. J 

w a is a complex 3-dimensional vector, with the bar denoting complex conjugation, w • w = w a w a 
and (w x w) a = e a i, c w b w c . The general 0(5) group element contains ten independent real functions. 
These are collected as six from M (in the three complex functions w a ), three in h, and one in k. 
The 0(3) gauge potential given by — tr(I a g~ 1 dg) reads 

A a = R ab (h) a b + (h^dh) a 

w a dw • (w x w) + w a dw • (w x w) 



w ■ w — |(w x w) 2 I 2 



(7.4.21) 



+ e abc (dw b w c -w b dw c ) }. 
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where R ab (h) is defined by hl a h~ x = R ab h b and k does not contribute. A a is the /i-gauge transform 
of a a , which depends on six real parameters (w a ). The three gauge parameters of h G 0(3), along 
with the six, give the nine functions needed to parameterize a general 0(3)- [or SU(2)-} potential 
in three dimensions. The Chern-Simons form is 



CS(A) 



167T 2 
1 

16V 2 



(A a dA a +\e abc A a A b A c ) 



'a a da a 1 1 



: (dh h- l ) a a a ^ 



167T 



24vr 



; tr(h~ l dhf. (7.4.22) 



The second equality reflects the usual response of the Chern-Simons density to gauge transforma- 
tions. Using the explicit form of a a as given in (7.4.21), we can further reduce this. Indeed we 
find 



a a da a +\e abc a a a b a c = (-2) 



x dw) • p + (w x dw) • p 

[1 + W ■ W — ~(w X w) 2 ] 2 



(7.4.23) 



Defining an Abelian potential 



Pk = \tijk dw 1 dw 3 



w • dw — w • dw 



1+w-w— ^-(wx w) 2 ' 
we can easily check that a da reproduces (7.4.23). 



CS(A) 



16vr 2 



-a da + d 



(dh h~ l ) a a a ) 
\M 2 



48tt : 



,-tr(^ 1 dhf 



(7.4.24) 



(7.4.25) 



If desired, the Abelian potential a can now be written in the Clebsch form making a da into a total 
derivative, while the remaining two terms already are total derivatives, though in a "hidden" form 
for the last expression. This completes our construction. 

(ii) Total-derivative form for the Chern-Simons term on a Kahler manifold 

In 4-dimensional space (xi, X2, £3, £4) we can introduce complex coordinates (holomorphic and anti- 
holomorphic) appropriate to Kahler (even-dimensional) manifold. 



(z, z) — (x\ ± ix<z), (w,w) = (x 3 ±ix A ) 



(7.4.26) 



It is then required that the holomorphic and anti-holomorphic components of the curvature F^ u 
vanish. [68] 

F zw = F 2 w = (7.4.27) 

It follows that 



A z = U- x d,U, 

A- z = -(d-M)(u-y, 



A w = U' 1 d w V, 

A^ = -(d^)(v-y. 



(7A.21 
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and after a further complex gauge transformation we may replace (7.4.24) by 

A g = g~ 1 d x g A w = g~ 1 d w g 

A- z = = A w (7.4.29) 

with Hermitian g = WU . When we define d± as the holomorphic and anti holomorphic derivatives 
d+ = dz-^ + dw-J^, oL = dz^= + dw£=, then A = A z dz + A w dw + A 2 dz + A^dw = g -1 d + g = a [with 
requirement (7.4.27) and in the selected gauge (7.4.29). This has the property of Abelianizing the 
gauge expressions. 

F = d_a, trF 2 = tr(d_ada) = d_tr(ad_a) (7.4.30) 

The last formula is also consistent with the Chern-Simons formula, since the cubic contribution in 
(7.3.1) vanishes, leaving the Abelianized quantity 

CS = tr(ad_a). (7.4.31) 

Additionally one can show that 

CS = tr(ad_a) = d + tt + cL$. (7.4.32) 

[In fact the above is established by varying tr(adja)] In this way one arrives at the final result [68] 

trF 2 = d_d + tt (7.4.33) 

It is to be emphazized the that (7.4.32) holds only with the holomorphic restriction (7.4.27), 
and gauge choice (7.4.29). No explicit formula for Q is available; it is determined by "intergrating" 
a variation: 

SSI = Itrad.v v = g~ l 5g (7.4.34) 

For $, which however does not contribute to trF 2 , but is needed for a reconstruction of the Chern- 
Simons term (7.4.32), we find 

5$ = tr(ad + v) (7.4.35) 

Evidently this construction does not give a parameterization for an arbitrary vector potential, only 
one which can be gauged to (7.3.28). 

It is instructive to see some explicit expressions. In the Abelian case, with A = a = d + 8, we 
have, F = d_d + 6, F 2 = (d_d + 6)(d_d + 6) = d-(d + 6d_d + 6),CS = AF = d + 6d_d + 9 = d + (6d_d + 6), 
Q = 9d-d + 6, $ = 0. Note that the Abelian connection is parameterized in terms of one function 
9, or two functions if the gauge freedom is included. But a general connection in 4-space, requires 
four functions. 

In the non- Abelian SU(2) case, we take g = e~ 0a , then 

A = a = g^d+g = (d + (9 a sinh0) + (cosh 9 - l){ie abc d + 9 b 9 c - 9 a d + 9))^- (7.4.36) 
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with 9 = y/ 9 a 9 a , 9 a = 9 a /9. The Chern-Simons term is reconstructed from 




(7.4.37b) 



(7.4.37a) 



Note that three functions are involved (9 a ) in the parameterization of the connection; six if the 
gauge freedom is included. However a SU(2) connection on a 4-space requires twelve functions. 
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